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Calibration of Precision Airplane Mapping Cameras 
By Francis E. Washer and Frank A. Case 


An instrument is described that permits the registration of all the information necessary 


to calibrate a precision airplane camera on 


single negative. Twenty-five collimators 


arranged at 7.5° intervals along two meridians at right angles provide optically distant tar- 


gets. These targets are photographed by means of the camera to be tested. 


Measurements 


made on the negative yield information on the equivalent focal length, distortion, resolving 


power, prism effect, orientation of the lines joining opposite pairs of collimation index 


markers, and location of the principal point. 


This instrument was designed and built to 


make the calibration of precision cameras required for all precision cameras used in Govern- 


ment mapping projects. 


A brief account of the calibration of a typical camera and a dis- 


cussion of the physical significance of calibrated focal length are given. 


I. Introduction 


mapping 
Bureau of 


The number of precision airplane 
cameras submitted to the National 

Standards for calibration has increased steadily 
since the first formal specification covering this 
type of camera was promulgated by the United 
States Department of Agriculture in March 1940. 
1] The demand 
mapping cameras has increased so rapidly since 


for calibration of precision 


1945, that new equipment [2] for this work has 
been developed and built to supplement the pre- 
cision lens testing camera that has been in use 
since 1935 3]. 
was primarily developed for use in photographi- 
length, 


The precision lens testing camera 
cally determining the equivalent focal 
distortion, and resolving power of lenses mounted 
in barrels or shutters. The instrument was later 
adapted to measuring these optical quantities for 
lenses mounted in cameras and further to locate 
ihe principal point in airplane mapping cameras 
: 

With increased volume of work, it was evident 
that the method of calibration involving use of 
he precision lens testing camera possesses several 
lcficiencies, First, four negatives are required, 


he making of which is a time consuming opera- 


on. Second, different sizes and shapes of the 


Ss in brackets indicate the literature references at the end of this 
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various makes of cameras cause difficulties in 
mounting them for test. Third, the entire field 
of the camera cannot be covered with a single 
photograph. The new camera calibrator was 
therefore designed to remedy these difficulties. 
Inasmuch as the optical characteristics of lenses 
intended for use in airplane mapping cameras can 
be measured on the old precision-testing camera 
prior to their installation in cameras, the new 
calibrator was designed primarily as a pointing 
instrument. It provides 25 beams of parallel 
light from known directions. The information 
derived from it is therefore mainly of a directional 
nature and includes data on the uniformity or 
lack of uniformity in the distortion characteristics 
of lens-cone combinations. It permits more 
rapid and accurate location of the principal point, 
together with a determination of the equivalent 
focal length of the lens as mounted in the camera. 
It provides quantitative information on the 
magnitude of the prism effect and tangential dis- 
tortion. Moreover it provides all this information 
on a single negative, thus avoiding possible changes 
that result from small movements between suc- 
provides a more stable 


cessive exposures. It 


support for the camera during test and also 
greater versatility in the types of cameras that 
can be tested. 

The scope of the present paper is limited to a 
description of the instrument and the instruments 





Fiat RE l. Collimator bank of camera caltbrator. 


Collimator bank before installation in the instrument. It consists of 25 c« 


limators arranged in the form of When mounted in the camera 


calibrator, the central collimator points vertically upward with the remaining 


- «0 ® 


24 collimators spaced at 7.5° intervals along the two meridians 


developed for use in its calibration. In addition, 
the calibration of a precision airplane mapping 
camera as conducted on the new calibrator is 
described. The manner of interpreting the nega- 


tive is given for completeness 


II. Description of the Camera Calibrator 


The camera calibrator and the details of its 
construction are shown in figures 1 and 2. Figure 1 
shows the heart of the new instrument, which is 
a bank of 25 collimators arranged in the form of a 
cross. In use, it is mounted beneath a table as 
shown in figure 2. The central collimator points 
vertically upward, and the remaining 24 collim- 
intervals along the four 
to 45°. Departure from 


interval used in the Bureau’s lens testing 


ators are spaced at 7.5° 
arms of the cross from 0 
the 5° 
camera [3] is necessitated by space limitations. 
The collimators are of the fixed-focus type with 
the tubes cut to proper length so that each reticle 
lies in the focal plane of the corresponding objec- 
tive. The collimators are mounted in a special 
casting, whose inner and outer envelopes are 
spherical. The collimator axes are normal to the 
correspondingly point 
Accuracy of 


spherical surfaces and 
toward the center of curvature. 
pointing of the collimators depends to a large 
extent upon the accuracy of counterbored holes 
and flanges by which the collimators are mounted, 


2 


but the reticles may be moved about a little ir 
the focal planes of the collimators so as to correct 
small deviations in the pointings. These adjust 
ments are made with capstan screws, so arrange: 
as to afford complete translational freedom of the 
reticle. After adjustments, clamping by appropri- 
ate screws guards against any subsequent move- 
ment. The clear apertures of the collimators ar¢ 
18 mm in diameter, which is somewhat smaller 
than those of the lenses commonly used on pre- 
cision mapping cameras. This, however, does not 
impair the instrument so far as its main purpose is 
concerned. ‘To guard the collimators against 
disturbance the lamp housings, mounted inde- 
pendently on the arms of the cross, completely 
surround the target end of the collimators but 
do not touch them at any point. A 6-v. frosted 
flashlight bulb serves to illuminate the reticle. 
Provision is made for the insertion of a filter and 
glass diffusing screen between the lamp and the 
target. 

The complete instrument is shown in figure 
The table top is a flat steel plate having a large 


») 


Figure 2. Camera calibrator. 


Instrument as arranged for use. The camera under test points vertically 
downward. Its proper orientation is effected by the autocollimating 
telescope and the plane-parallel piece of glass placed on the focal plane of 
the camera. Images on a finished negative appear every 7.5° along each 
diagonal, with the collimation index markers registered on each side of the 


negative 
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ircular opening in its center. The collimator 
ank is mounted beneath the table and centered 
ith respect to this opening in the table top. 
(he camera holder consists of a tripod placed 
wer the central opening in the table. A circular 
pening in the center of the tripod permits the 
echt from the 25 collimators to reach the lens 
any camera placed on the holder. The focal 
plane of the camera can be adjusted to normality 
with the axis of the central collimator with the 
iid of the adjusting screws on the tripod and the 
sutocollimating telescope, whose axis is bent 
“0 by the prism mounted in front of its objective. 
he adjusting screws also provide a limited amount 
vertical adjustment to compensate differences 

, distance between the front surface of the camera 
lens and the front end of the camera of the various 
types of camera. Four flashlight bulbs mounted on 
adjustable arms reaching out into the circular 
opening of the table top serve to illuminate the 
collimation index markers of the camera under test. 


III. Test Chart 


The test charts of the new instrument include 
resolving power patterns, although the limited size 
of the collimator objectives precludes the evalua- 
tion of the resolving power of the camera lens at 
is maximum aperture. However, the inclusion 
of the chart does permit study of the resolving 
power at reduced aperture. Since a new test 
chart for use in the precision lens testing camera 
had been under consideration for some time, it 
was decided to make a new chart at this time and 
to use a reduced version of it as the target reticle 
of the camera calibrator. This new test chart is 
shown in figure 3. The chart was first made on a 
large seale and then photographed to the desired 
size, using Eastman high resolution plates, type 
(49GH. The size of the target in the 0° collimator 
:Smm square. Large scale charts are also made 

nd subsequently reduced for each of the various 

ngular separations from the axis, so that the 
osine correction in the vertical dimensions and 
osine-squared correction in the horizontal dimen- 
These 


orrections are necessary in order that the corre- 


on of each target pattern can be made. 


nding patterns will all be imaged at the same 
ce regardless of angular separation from the 
“is for a given lens-camera combination. The 
hree-line patterns reduce in size in a geometrical 
ogression, the ratio being ¥/2, or approximately 
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Closer steps in the series are admittedly 


1.189. 
desirable, for example ¥/2, but this smaller ratio 
could only have been achieved at the expense of 


compressing the range of permissible values be- 


cause of the small size of the reticle. The spacings 
present in the patterns in the reticle range from 
2.3 to 161 lines per millimeter in 26 steps. These 
patterns, when imaged by the lens under test, are 
changed by the ratio of the focal length of collima- 
tor lens and lens under test. Table 1 shows the 
range of resolution in the image plane for the three 
focal lengths most frequently encountered. The 
nominal maximum f-number of the lenses are 
given. Since the diameter of the collimator lens 
reduces the effective aperture of the lens under 
test, the effective operating f-number is given for 
each lens. The line showing limit of resolution 
shows the maximum theoretical resolving power 
for each lens at its effective f-number. With this 
chart therefore, the upper limit of resolution is 
determined by the performance of the lens and 
limitations of the emulsion of the photographic 
plate and not by insufficient range of the resolving 
power patterns. A considerable amount of work 
has been done in other laboratories, wherein a 
circular target is used instead of a line target. 
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FIGURE 3 Test chart. 

rhe test chart that is used in preparing the reticle for use in ¢ 
consists of patterns of parallel lines in two orientations with sp 
in geometric progression by steps equal to w/3. The circle targets 


for comparative study of line and circle targets 








TABLI l Limits place | on effect ‘ f-n embe r by the diam- 
eter of the collimator objectives of the camera calibrator for 
three lenses of different focal lengths; the maximum 
theoretical resolution of these lenses operating at these 
f-n umbers, and the range of the re solving powers obtain- 
able for these different focal lengths provided by the test 


charts of the camera calibrator 


Equivalent focal length, mm 130 1590 210 
Usual nominal f-number 6.3 6.3 6.8 
Effective f-number 7.2 8.3 11.7 


Maximum theoretical resolution in lines 

per mm, as determined by effective 

f-number 1s 172 122 
Minimum resolution in lines per mm 


provided by test chart 1.7 4.1 2.9 
Maximum resolution in lines per mm 
provided by test chart 32 289 205 


A series of circular patterns is therefore included 
in the test chart for comparison of results obtained 
with it with the results obtained with the 3-line 
patterns. 

The cross in the center of the chart is the fiducial 
mark with reference to which the angles separating 
When imaged on 
a camera test negative, the measured separations 


the collimators are measured. 


of the crosses for various collimators serve for the 


determination of equivalent focal length and 
distortion. 

The National Bureau of Standards is at present 
using a high-contrast target and high-contrast 
fine-grained plates in the testing of photographic 
This 


is being done because existing Government speci- 


objectives and airplane mapping cameras. 


fications require that lenses for use in aerial map- 
ping projects be tested under these conditions. 
It is recognized that there may be merit in the 
procedures followed in other laboratories that 
require the use of low-contrast targets and fast 
emulsions comparable to those used in actual 
aerial photography. Investigations are now in 
progress at this Bureau on the effect of contrast 
When this 
investigation is concluded, it is probable that a 
different type of target will be recommended for 
use in the testing of photographic objectives and 


of target on performance of lenses. 


cameras with corresponding changes in the re- 
quired values of the resolving power for certifica- 
tion of lenses for Government use. 


IV. Calibration of the Instrument 


Adjustment of the angular relationships and 
measurement of the angles separating the colli- 


4 





mators in each meridian is necessary before the 
instrument can be used in the calibration of pre- 
cision cameras. 


1. Adjustment of Angular Relationships 


It is required that the angles between corre- 
sponding collimators on opposite sides of the center 
be as nearly equal as possible, and moreover that 
the axes of the beams of parallel light emerging 
from the collimator in a given meridian be parallel 
to and in near coincidence with a common plane. 
This adjustment is made before the collimator 
bank is mounted under the table. The collimator 
bank is placed face downward over the opening 
The lamp 
housings are removed to expose the target reticles. 


in the table top, as shown in figure 4. 


A front surface mirror is placed at the center of 
curvature of the spherical casting. The mirror, 
shown at A, is mounted in a special device that 
permits easy leveling and rotation about a hori- 
zontal axis. 

A microscope equipped with a vertical illumina- 
tor is mounted over the reticle of the central 
collimator, and the mirror surface is brought into 
a position normal to the axis of the central 
collimator by means of the leveling screws, with 


the collimator, mirror, and microscope serving as 
If now the —7.5° 
one side is illuminated with the desk lamp, the 


an autocollimator. reticle on 





Figure 4. Arrangement of apparatus for adjusting the 


angular relationships. 
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second microscope may be used to view the reticle 
of the +7.5° collimator on the opposite side while 
idjusting the reticle to bring it into coincidence 
with the reflected image of the —7.5° reticle. 
The two 7.5° angles are then equal, and the axes 
° collimators fall in 
The microscope with the vertical 


-O 


f the +7.5°, 0°, and —7.5 
one plane. 


lluminator is then brought over the +7.5 


reticle and the mirror rotated and made normal 
to the axis of this collimator. The +15° and 


22.5° collimators are then adjusted. The proc- 
ess is repeated by using the + 22.5° collimator as 
l 30° 97 FO 


the autocollimator to adjust the +-37.5°, 


and +-45° collimators. On returning and using 
the 0° collimator as autocollimator, all of the 


remaining collimators on the same meridian can 
be adjusted. Following this adjustment it can 
be stated that the axes of all collimators in a 
given meridian lie in the same plane, opposite 
angles are equal, and moreover angles between 
adjacent collimators are equal. The entire process 
is repeated to adjust the collimators on the meri- 
dian at right angles to the first meridian. It 
must be mentioned that, although the angles 
separating collimators along a given meridian 
are equal, the corresponding angles between 
collimators in different meridians are not exactly 
equal. The accuracy of the initial boring of the 
casting is attested by the fact that no difficulty 
was encountered in equalizing the angles and 
making the collimators coplanar, although the 
amount of lateral adjustment available was less 


than 1 mm. 


2. Measurement of the Angles 


The angles between adjacent collimators are 


determined by a comparison method using a 


reflecting biprism and telescope, as shown in 
figure 5. The standard of comparison is a biprism 
having a known angle between the two surfaces 
forming the obtuse angle. The two surfaces are 
aluminized and form two mirrors that maintain a 
onstant angle with respect to one another. The 
angle formed between one surface and the exten- 
sion of the other is 3.7600°. 


hould be one-half of the 7.5° 


Ideally this angle 
angle between the 
ollimators, but its present value is sufficiently 
lose, and moreover the amount of departure from 


75° is known. The biprism is mounted in a 
pecial device that permits rotation of the biprism 


bout a horizontal axis coincident with the edge 
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separating the two aluminized surfaces. As shown 
in figure 5, this device is mounted on the camera 
holder in such a manner that the edge of the 
biprism coincides with the center of curvature 
and is normal to the plane in which the collimator 
axes for a given meridian lie. The error that 
arises from unsymmetrical use of the collimator 
objectives is made negligibly small by the use of 
collimator objectives that have very small longi- 
tudinal spherical aberration and making certain 
that the reticles lie in the focal plane of the 
objectives. 

To measure an angle, the telescope, equipped 
with a micrometer for varying its pointing, is 
sighted at the center line of the biprism. The 
biprism is rotated about its edge until the image 





Reflecti ng bipr ism for measurement of angle S 


Figure 5. 


reflected by one face of the biprism of the reticle 
in one collimator coincides with the cross-hairs 
of the 
micrometer reading is recorded. The reticle of 


viewing telescope. The lower telescope 
the first collimator is darkened and the second one 
illuminated. The image reflected by the second 
face of the biprism then appears in the field of 
view of the telescope. If the angle 
collimators is exactly twice the biprism angle, no 


between 
shift is evident. If the second image is shifted, 
the lower telescope micrometer is used to correct 
the pointing and the difference in micrometer 
reading noted. The 
reading for the two pointings when translated 


difference in micrometer 


into angular shift gives the difference in angle 
between the telescope and twice the biprism angle 


By rotating the biprism, the angles separating 
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TABLE 2. Values of the individual angles separating adjacent collimators and their deviation from 7.5° for each collimator bank 


Collimator bank 


Nominal angle Deviation Angle 


7. 4046 
7. 4978 
7. 4960 
7. 4955 
7. 4930 
7. 4992 


determined for 
The 


biprism holder and viewing telescope are then 


each successive collimator can be 
all the collimators lving in the meridian. 


rotated 90° about a vertical axis, and the above 
process is repeated to obtain the angular separa- 
tions of the collimators in the second meridian. 
The measured values of the angles are given in 
table 2 to show the agreement between the various 


angles, together with their departure from 7.5° 
The tolerance on angle for boring the casting was 
1 min. (+0.0167°); that the 


results obtained indicate the accuracy of boring 


the table shows 


was amply close. 


V. Operation of the Camera Calibrator 


To calibrate a camera on the instrument, one 
must be sure that the focal plane of the camera is 
To 


do this the operator adjusts the micrometers on 


normal to the axis of the central collimator. 
the viewing telescope until the image of the 
collimator coincides with 
The 


camera cone is then placed on the camera holder, 


center cross on the 0 
the cross-hairs of the viewing telescope. 


adjusted laterally until the lens of the camera is 
coaxial with the 0 The 
then adjusted vertically until the images formed 
of 


vignetting, the camera cone is rotated until the 


collimator. camera is 


by the outer collimators show a minimum 
two diameters along which the two perpendicuiar 
rows of images fall coincide with the diagonals of 
the focal plane. The camera is then clamped to 
further An optically flat 
plane-parallel glass plate aluminized on its upper 


prevent its motion. 
surface, is placed on the focal plane of the camera. 
Then using the viewing telescope as an auto- 
collimating telescope, the focal plane of the cam- 
era is adjusted to normality with the axis of the 
autocollimating telescope by small adjustments of 
the screws supporting the camera holder. Flash- 
light bulbs attached to the ends of small rods 
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supported by the table top are then moved about 
until the light from the bulbs illuminate each of 
the four collimation index markers mounted in the 
The rods are then clamped. The 
plane parallel is removed; the room is darkened, 


camera cone. 


and a photographic plate is placed on the focal 
plane of the camera, and an exposure is made. 
It is customary to place a heavy piece of flat glass 
over the photographic plate during exposure to 
hold the plate in place and help to flatten it. 

The finished negative shows the images of the 
targets every 7.5° from the center to the corner 
A schematic drawing of a 
finished The 
collimation markers are registered on the sides of 
the negative, so that there is little likelihood of the 
that 


along each diagonal. 


negative is shown in figure 6. 


target images being fogged by the light 


registers the markers. 


Fiaure 6. Schematic drawing of test negative obtained with 


camera calibrator. 
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Eastman V-G spectroscopic plates are generally 
ised. In those instances, where a_ precision 
amera is to be used for infrared photography, 
Eastman IV—R spectroscopic plates are used. 


VI. Calibration of a Precision Camera 


the location of the principal point with respect 
to the collimation index markers, determination 
of the angle formed by the intersection of lines 
oining opposite pairs of collimation index markers, 
determination of the equivalent focal length of 
the lens as mounted in the camera, determination 
of the distortion, and resolving power from the 
center to the corner of the image plane. Detail 
requirements relating to all these quantities are 
contained in specifications set up by various 
governmental mapping agencies concerning cam- 
eras that are to be used in Government mapping 
projects [5]. It is further stipulated that precision 
type cameras intended for use in these projects 
be tested for compliance with specifications by 
this Bureau. 

It frequently happens that a camera on final 
examination does not comply with the require- 
ments, but that it can be brought into compliance 
with the requirements by small lateral movements 
of the lens with respect to the collimation index 
markers or by small movements of the markers 
themselves. Following these adjustments, a re- 
check of the camera shows that it complies. To 
ensure that nothing can then happen to change 
this adjustment, holes are drilled and dowels 
placed so that the movable members remain fixed 


with respect to one another. 


VII. Interpretation of the Negative 


With the new camera calibrator, all of the 
nformation necessary to evaluate performance 
of the lens-cone combination is contained in a 
single negative. However, it is customary to 
nake two negatives and combine the results of 
measurements from these to minimize any pos- 
sible error. A schematic drawing of a typical 
negative is shown in figure 6. The markers are 
designated A, B, C, and D for convenience. The 
etter A is reserved for the markers bearing the 
urow indicating direction of flight. The Roman 
umerals J, JJ, IJJ, and IV serve to designate 


he rows of images formed by the four collimator 
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Calibration of a precision camera consists of 
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FIGURE 7. Equi pme nt used in the determination of the 90 


condition. 


banks radiating out from the center. The central 
image formed by the 0° collimator is designated 
the center cross (C. C.). The intersection of 
lines joining opposite pairs of collimation index 
markers is designated the center of collimation, 
(C of C). 


1. Determination of the 90° Condition 


The first measurements on a given negative are 
made to determine whether or not the lines AB and 
CD are perpendicular within +1 min. This deter- 
mination is made on a special device shown in 
figure 7.2. The central feature of this device con- 
sists of a piece of plate glass bearing four short 
radial diamond lines cut in the surface near the end 
of a 4.5-in. radius. The error in the positioning 
of the lines does not exceed +:5 sec. This plate is 
embedded with the lines uppermosé in the metal 
frame that bears the plate clamps, so that the sur- 
face of the metal frame and the upper surface of 
the glass lie in the same plane. The negative is 
placed, emulsion side down, over the glass plate in 
such a manner that the index lines of the collima- 
tion markers and the four diamond lines are in near 
coincidence. The negative is then clamped down, 
so that no relative movement can occur during 
measurement. The metal frame is mounted on a 
spindle, so that a given marker can be brought 
under the viewing microscope. The separation of 


2 The authors wish to express appreciation to William P. Tayman, who 
developed this device for holding the 90° standard plate and the negative in 
t 


contact during measurement, which represents a definite advance over the 


earlier equipment used in making this determination 














diamond line and index line on the negative are 
then measured. This process is repeated for the 
remaining three markers. From these four meas- 
urements, the departure and the direction of de- 
parture of the collimation index markers from the 


90° condition can be determined 


2. Location of the Center Cross with Respect to 
the Center of Collimation 


The center cross, as used herein, is the point 
where an infinitely distant object point lying on a 
line perpendicular to the focal plane is imaged by 
the camera lens. The coordinates of the center 
cross are determined with respect to the rectan- 
gular coordinate system formed by the lines 
joining opposite pairs of collimation index markers 
These 
measurements are made with the aid of the device 
figure 8. The 


having the center of collimation as origin. 


shown in negative is placed, 


emulsion side uppermost, upon the flat metal 
plate and so oriented that the line AB is ap- 
parallel to the bench ways upon 
rests. The then 


clamped to prevent any movement with respect to 


proximately 


which the slide negative is 


the metal plate. By moving the slide along the 
bench ways, the measuring microscope can be set 
in turn on marker A, center cross, and marker B 
and readings taken. From this data the lateral 
displacement of the center cross from the line AB 


ean be obtained. The negative is then rotated 


Ficure 8. Equipment used in locating the center cross with 


respect to the center of ( ollimation. 

















| 
4. 





Ficure 9. Schematic drawing showing the center cross with 


respect lo the center of the collimation. 


90° and similar settings made to determine the 
lateral from the line CD. The 
results of measurement on a typical negative are 
The magnitudes 


displacement 


shown schematically in figure 9. 
of the displacement are exaggerated for purposes 
of clarity in the figure. 


3. Location of the Principal Point with Respect to 
the Center Cross 

The principal point of precision-type aerial 
mapping camera is defined as that point where a 
perpendicular dropped from the near nodal point 
of the lens meets the focal plane [6]. For an ideal 
lens, the principal point and the center cross 
coincide. In practice, however, the center cross 
is shifted away from the principal point either 
because of prism effect in the lens or because of 


nonparallelism of the surfaces of the filter on the 
front of the lens, which forms a part of the optical 


This is illustrated in 


figure 10, which is a schematic drawing showing 


system of the camera. 


the displacement of the axial ray and rays inclined 
at angle 8 with the axis before and after placing a 
prism in front of an ideal lens. 

In figure 10, O is the principal point and is the 
point where light from an infinitely distant object 
point lying on the axis of an ideal lens would be 
imaged if no prism effect is present. On inter- 
posing a prism having a refractive index of 1.5 
and angle @ the axial ray is bent away from the 
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ormal by an amount ¢=a/2 and now cuts the 
wal plane at O’, which is now referred to as the 
nter cross. To locate the position of the 
rincipal point with respect to the center cross the 
stance O’O must be determined. 
The determination of O’O is made possible be- 
use the rays inclined at angle 8 are also displaced 
y amounts ¢«, and ©, and cut the focal plane at 
points .Y, and .Y, instead of Y, and X,. By 
ssuming that the ray passing through the prism 
ind meeting the focal plane at 0’ is near the region 
of minimum deviation, it follows that we may 
egard «, as equal to & and can compute the rela- 
tive magnitude of € in terms of «. The magni- 
tude of € is the average value of « and & and al- 
than « and X, X,=X, X,>0’0. 


None of these quantit ies can be measured directly, 


vavs greater 


it their magnitudes may be obtained as follows: 


From the figure, 


OO’ X,=f{tan (8+-«)—tan «| —=ONX,—O'O 

OX, =Jl{tan (8B—e) +tan «|=ON,+0'0 
or O'X}=OX,4+ X; X,—0'O 

O'XN’=OX TY, A,+-00. 





O X,—O'X,=2(X, X,—O0’0)=AD 


X, A2.=X, NX, and 0X,=0OX, 


1 
In the case of an actual lens having prism effect, 
Y,, 0’, and YY, are the observed images and the 
separations O’ Y, and O’'X, can be 
Consequently, the value of AD) can be determined 


measured. 


for each pair of equal known angular separations 
from the axis from measurements on the negative. 
This information is used most conveniently by 
comparing the measured values of AD) with values 
of AD computed for a definite set of values of 
angle,a, focal length /, and separation 


angle ~B. 
formed and the observed values of AD can be 


prism 
The computations can be readily per- 


nterpreted in items of q@ sufficiently accurately 
by this process, because AD is linear in a@ for a 
fixed angle 8 for small value of prism angle a. 
lable N3X2, 


VO, and AD for the case of a lens having a focal 


3 gives the computed values of & 


ngth of 150 mm combined with a prism of index 
> having an angle @ of 0.05°. It may be noted 
om the table that AJ) increases rapidly with 
creasing 8. For the case shown, AD is negli- 
ble at 7.5° but increases to 0.3 mm at 45°. It 
evident from the table that the value of O’O 
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Ficure 10 


prod iced by inte rposing a thin prism in 


Schematic drawing showing the image shift 
ont of an ideal 


lens. 





deduced at the wider angles is more sensitive than 
the value determined at the narrow angles. It is 
also clear that the presence of appreciably large 
prism effect may cause marked asymmetric dis- 


tortion. 
Tassie 3. Image dis place ments of various points in the 
focal plane for a prism of index 145 anda ngle a 0.05° 


combined with a lens whose focal le ngth is 150 mm 


) 
8 © |X8X=| oo | ox 
oN 
Degrees Dearees mm mm mm 
0 0. 0250 0. 06, 0. 065 0. 000 
7 0254 OOS Of 006 
1 on o74 * OS 
22 0285 OS6 nt 042 
i) OS16 110 (nt O90 
7 O364 151 Ot 172 
4 043 228 Ot 2t 


Although it is feasible to determine the prism 
effect 
calibrator by this process, it is customary to 


from a single negative from the camera 


combine the results from two negatives with the 
camera rotated 180° 
symmetry between the marking of each negative. 


being about its axis of 


This procedure increases the accuracy and simpli- 
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fies some of the computations. The coordinates 
of the principal point P. P., with respect to the 
center cross, are measured along the diagonals 
of the focal plane, as indicated in figure 11, 
because the calibrator is so constructed that the 
images from the various collimations lie along 


the diagonals. 


4. Location of Principal Point with Respect to 
Center of Collimation 


To locate the principal point with respect to 
the center of collimation, it is only necessary to 
combine the results obtained in section VII, 
2 and 3. Care must be taken to avoid errors 
in sign and to make allowance for the fact that 
the measurements are made in two coordinate 
systems rotated 45° with respect to one another. 
Figure 12 shows the resultant obtained by com- 
bining the measurements illustrated in figures 
9 and 11. 

The usual requirement in precision mapping 
cameras is that the principal point must not be 
distant from the center of collimation by more 
than 0.03 mm. In the event that this require- 
ment is not met, it is customary to adjust the 
markers with respect to the lens until the separa- 
tion of principal point and center of collimation 
is less than 0.03 mm. When this has been done, 
dowels are placed so that no relative movement 
of principal point and center of collimation can 


occur 














w 0 
B 


Figure 11. Schematic drawing showing the principal 


point with respect to the center cross. 
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Figure 12. Schematic drawing showing the principal 


point with respect to the center of collimation. 


5. Determination of Equivalent Focal Length and 
Distortion 


The same negative from which the location of 
principal point is made is used to determine the 
equivalent focal length and distortion. In fact 
the same measurements, described in section VII, 
3, for use in locating the principal point with 
respect to the center cross are used. 

The equivalent focal length [6] is 
theoretically by the equation 


defined 


f=lim 


y' ' 
{ 
so tan B 


where y’ is the distance from the principal focus 
to the center of the image in the image space 
focal plane of an infinitely distant object point 
that lies in a direction making an angle 8 with the 
axis of the objective. If a photographic objective 
were free from distortion, the quotient would be 
invariant with respect to 8. For many photo- 
graphic objectives, the distortion is negligible for 
points distant from the center of the useful field 
by not more than one-fifth of its radius. Con- 
sequently, it is often possible to obtain a satis- 
factorily accurate value of f by a single deter- 
mination of 8 and y’ for a point lying near the axis. 

The negative made for an airplane camera cone 
on the camera calibrator has images at the center 


of the field and at six known angles proceeding in 


7.5° steps from the center to the edge of the useful 
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field. Consequently, on substituting the meas- 


red separation from the 0° to 7.5° image and the 
weasured value of tan 8 into eq 1., the value of the 
quivalent focal length can be determined. The 
to 7.5° separation occurs four times on the nega- 
ve along radii separated by 90° of azimuth, so 
hat four independent determinations of f are pos- 
ble from a single negative. The accuracy of this 
etermination is further increased by the fact that 
he measurements made on opposite sides along 
he same diameter automatically compensate for 
ny error of plate tipping. It must be mentioned 
that practically the term equivalent focal length 
s used for a lens cone combination is actually the 
scale factor for use in interpreting distances meas- 
ired on our aerial photograph in the central area. 
Inasmuch as photographic objectives are mounted 
in the aerial camera in such manner as to yield 
wst over-all definition, it is clear that the equiva- 
lent foeal length as above determined is the scale 
factor for the plane of best average definition. 
Having established the value of the equivalent 
focal length for the lens camera combination, the 
distortion referred to the equivalent focal length 


can be readily found. To evaluate the distortion 


Y., Yo, Yo, . . ~ be the separation of the 
images on the negative from the central image 
corresponding to the angles im the object space 
B,, B, Bs. Let products f tan 8,, f tan B, f tan Bs, 
be computed. Then the values of the distortion 


D),, Dy, Ds 


, are given by the relation 
D,=y,—f tan 8, 
D,=y,—f tan p, 


D,=y,—f tan Bs 


Positive values of the distortion indicates a 


displacement of the image away from the center 


of the negative. 


6. Determination of Calibrated Focal Length 


The lenses used in aerial mapping are not free 
from distortion. Consequently the value of the 
quivalent focal length determined for the vicinity 
of the axis is not necessarily the best scale factor for 
ise in interpreting measurements made on the 
negative at points well removed from the axial 

ca If one holds fast to the value of equivalent 
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focal length determined for the axial region, the 
measurements made on the negative in the axial 
region are free from distortion, whereas in the 
extra-axial regions the errors from distortion may 
become serious. 

The customary procedure for alleviating this 
condition is to use a new scale factor that prov ides 
a better over-all accuracy in interpreting distances 
measured on the negative, although it does so at 
the expense of introducing negligibly small errors 
in the axial region This new scale factor is called 
the calibrated focal length [6]. The term eali- 
brated focal length, its significance, and the man- 
ner of its evaluation have been sources of con- 
siderable misunderstanding in photogrammetric 
circles since the inception of the term. It is 
proposed in the following paragraphs to present 
a study of the variation of equivalent focal length 
and distortion that occur for a lens affected with 
distortion when the equivalent focal length is 
based on the expression 


, 


 . 
tan 6 
finite 


and 8 is allowed to take any value at 


intervals between 7.5° and 45°. It is hoped that 
careful consideration of this analysis will lead to 
better understanding of the physical significance of 
the calibrated focal length. 

In table 4, the row marked y’ lists the measured 
separations from the axial image of the images 
These 


measurements are from an actual negative made 


occurring at 7.5° intervals out to 45°. 


on a lens-cone combination. In the second row, 
values of f, the equivalent focal length, are given 
WwW here each \ alue of fis obtained from the equat ion 
y’ 
tan 3 


It is clear that a different value of f is obtained for 
each value of 8. This variation in f/ appears be- 
cause of the distortion in the lens. It is possible 
to determine the values of the distortion for each 
of the values of f at each angular separation from 
the axis. This has been done, and the results are 
shown in the table. One of the natural conse- 
quences of this procedure is that, for whatever 
value of the equivalent focal length is determined, 
the distortion is necessarily zero for that value of 


8 when the distortion is evaluated with respect 
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to that particular focal length. The location and 
the magnitude of the positions of maximum 
positive and negative values also change for each 
These effects are shown graphically in 
It may be noted in figure 13 that the 


value of f. 
figure 13. 
curves of distortion are very similar, and that the 



























































7° 5° 225° 30° 375° 45 

Figure 13. 

with equivalent focal length computed from the relation 
J=y'/tan B for varying angles of B. 


Variation of distortion in the same image plane 


The lowest curve showing the distortion referred to the calibrated focal 
length indicates that the calibrated focal length is equal to the equivalent focal 
length for values of 8 of approximately 16° and 42.5°. 
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Taste 4. Variation of the values of equivalent focal length 
and distortion with angle upon which the computations are 


based for a typical lens 


V5 


mm mm mm mim mm 
20. 064 40. 487 62. 182 88.112 


fi he fs 


y 


/=tan B’| 152.400 | 152.443 | 152. 535 


Distortion referred to the equivalent focal length f=y’/tan 8 


0. 000 0.012 0. 056 
006 . 000 038 
O18 025 000 
O28 O46 0383 

~, 025 039 “22 


007 026 o79 
Distortion referred to the calibrated focal length f-=152.451 mm 


0. 006 0. 002 35 0. 094 


variation of the curves is very like that which 
would be produced by rotating the entire curve 
above the 0° point. 

The map mi«er is primarily interested in the 
scale factor that introduces the least error wherever 
used throughout the entire area of the photograph. 
Thus, while use of the scale factor equal to the 
equivalent focal length at 7.5° gives zero error 


° 


within the 7.5° zone, it is at the expense of large 


P Inasmuch as there is 


error in the 37.5° zone. 
usually no preferred area in a given photograph, 
it is obvious that it would be preferable to use a 
scale factor that would reduce error in the region 
of large error even if a small error were introduced 
in the region of zero error. Considering the 
curves in figure 13, it is clear that use of J, instead 
of f; as a scale factor will give reduced distortion 
in the 37.5° zone with the introduction of a small 
amount of distortion at 7.5°. Selection of one of 
the other values of f changes the error pattern 
either for the better or for the worse. 

It is, at present, standard practice to consider 
the preferred error pattern as that one for which 
the absolute values of maximum positive and maxi- 
mum negative distortion are equal. The pre- 
ferred value of equivalent focal length for which 
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this condition prevails is called the Calibrated 
Focal Length. It seldom happens that the cali- 
brated focal length coincides with the equivalent 
focal length at one of the fixed angles (i. e. one of 


° intervals) so that it is necessary to com- 


he 7.5 
pute the increment Af that must be added to the 
equivalent focal length determined at one of the 


fixed angles (usualiy the 7.5° angle) to yield the 
calibrated focal length. This may be done with 
the aid of the following formula: 


pam T dD, 
tan 8,,+tan Bn 


where D,, and D, are the values of distortion 
existing at the angles 8, and 8, referred to the 
equivalent focal length f. The calibrated focal 
length f, is then obtained from the relation 


f=f+af. 


The above formula is a general one, and the new 
values of distortion referred to the calibrated 
focal length at angles 8, and 8, are equal in mag- 
nitude and opposite in sign. Usually 6, and 8, 
are the angles at which maximum positive and 
naximum negative distortion occur. These equa- 
tions are essentially the same as those reported in 
a paper by Sewell [8]. 

The values of the distortion referred to the 
calibrated focal length are shown in table 4 and 
figure 13. The computation is made for8,,=37.5 
and 8,=45°. It may be noted that the same 
value of f, will be obtained regardless of which 
value of f is selected as a base of computation, so 
long as the corresponding values of distortion are 
ised. 

It must be emphasized that the determination 
of the calibrated focal length involves no shift 
whatever of the position of the focal plane of the 
The cali- 


brated focal length is simply that value of the 


camera cone with respect to the lens. 


quivalent focal length that serves as the pre- 


ferred seale factor in interpreting distances 


measured on the photograph. 
7. Determination of the Resolving Power 


The camera calibrator provides an excellent 
neans of checking whether or not the camera lens 
; mounted in the camera is capable of producing 
isable definition from the center to each corner 
f the field. If, as sometimes happens, the resolu- 
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tion along the four radii from the center to the 
corners of the field is not uniformly good this will 
be detected on the negative and thus minimize the 
possibility of a camera being certified as satis- 
factory when the resolution satisfies minimum 
definition requirements along one diameter of the 
field but does not do so at some points along the 
Table 5 lists 


the observed values of the resolving power for a 


diameter at right angles to the first. 


typical 6-in. lens in the focal plane of the camera 
in which it is mounted. These resolvings powers 
are measured at effective aperture //8.3 and may 
be slightly higher than the values that would be 
found at maximum aperature. 

It is clear from this table that some variation 
exists in the values of the resolving power for the 
same angular separation from the axis on different 
meridians. For the lens shown, the variations 
are not such as to bring the resolving power down 
at any point to values lower than the usually 
specified minimum of 15 lines per millimeter for 
this type ot lens. 


TABLE 5. Observed values of resolving power for a typical 
6-in. lens as determined from a negative made with the 


camera calibrator 


Resolving power for tangential lines in 
per millimeter with angular separation 
the axis « 


Resolving power for radial lines in lines 
millimeter with angular separation 
the axis of 


8. Flatness of Platen 


The specifications for precision airplane mapping 
cameras [5] usually contain the requirement that 
the surface of the platen shall not depart from a 
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plane by more than +0.0005 in. The platen is 
the flat surface against which the film is pressed to 
ensure its planeness during the instant of ex- 
posure. This requirement is included in the 
specifications, because small departures from flat- 
ness in this locating surface would be imparted to 
the film and would in turn affect the location of 
images on the final negative. 

It is therefore customary in the course of 
calibration of a precision camera to check the 
platen for flatness. It has not been considered 
expedient at this laboratory to make a detailed 
contour map of the entire platen surface. Instead 
the practice here has been to check the platen for 
flatness along selected lines or diameters of the 
platen and to conclude that the platen is satis- 
factorily plane if no departure from planes in 
excess of the specified tolerance occurs along these 
lines. 

A simple device, shown in figure 14, has been 
constructed that permits a rapid measurement of 
the departure from flatness along a selected line 
+0.0001 in. The device 
consists of a metal beam supported by two legs 
Midway between these two legs, 


with an accuracy of 


set 8 in. apart. 
a sensitive dial indicator gage is mounted with the 
lower end of the plunger, which motivates the 
gage lying in the line common to the two main 


outer supports. The two outer supports and the 





Dial indicator device used for checking flatness 


Fiaure 14. 
of platen. 


The instrument is graduated in microns and permits the rapid meas- 
urement of departures from flatness with an accuracy of +0.0001 in. 
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central plunger have convex spherical surfaces at 
the lower portions, which touch the surface to be 
tested. The two smaller supports serve only to 
hold the instrument in a vertical position and to 
prevent the instrument from tipping over. The 
instrument may be rocked on the two outer legs. 
In use the small front leg is in contact with the 
surface being tested, and the movement of the 
plunger is in a direction closely perpendicular to 
the surface. 

The instrument is graduated in microns. 
Before each use, it is placed on an optical flat to 
ensure that the gage reads zero for a truly flat 
It is then placed on the platen under 
The instrument is 


surface. 
test and the reading noted. 
moved across the platen, rotated 90°, and moved 
It is also moved about the platen 
If the pointer on the gage 


across again, 
in other orientations. 
does not depart from zero by more than +134 
during the course of this operation, the platen is 
adjudged satisfactory. 


VIII. Discussion 


The method used at this Bureau for determining 
the location of the principal point of airplane 
cameras has been in use for some years. This type 
of work was first performed on the precision lens 
testing camera and it is now being performed on 
the new camera calibrator. In recent years, the 
question has been raised whether the principal 
point as herein defined and determined is the 
proper point to use in the photogrammetric sense 
in the interpretation of aerial photographs. In 
the absence of lens distortion, particularly of the 
asymmetric type induced by prism effect, there 
can be no question of the adequacy of the method 
and its interpretation. 

However, in the presence of unbalanced distor- 
tion, there may be some merit in the criticisms 
that have been directed at this concept of the 
principal point. It is not the aim of the present 
authors either to refute or concur in the various 
other suggested concepts of what constitutes the 
best definition of principal point in cameras in the 
presence of asymmetric distortion. Rather, these 
matters are included in the present paper to show 
that cognizance of these alternative procedures is 
being taken; and it is hoped that ultimately some 
common basis may be found wherein these differ- 
ences can be reconciled. 
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The most recent publication dealing with one 
of the alternative definitions of principal point is 
that of P. D. Carman of the National Research 
Council of Canada [9]. In that paper a proof is 
presented that indicates that a point having 
properties identical with the center cross as de- 
fined in section VII, 2 of the present paper, ought 
to be used as the principal point. This definition 
of principal point is accepted in Canada as evi- 
denced by a paper by R. H. Field [10], who de- 
scribes the method of locating the principal point 
ised at the National Research Laboratories at 
Ottawa. 

In the September 1948 issue of Photogrammetric 
Engineering, E. D. Sewell [8] presents the concept 
of a “Point of symmetry,” which is a point about 
which all radial distortions are to be symmetrical. 
This point does not coincide with either the center 
cross or principal point as defined in the present 
paper, although all three points coincide in the 
absence of asymmetrical distortion. 

To illustrate the differmg locations of these 
pomts, consider the following case. A camera 
whose lens showed negligible prism effect was cali- 
brated at this laboratory and the principal point 
and center of collimation brought into coincidence. 
Following this, a thin prism was placed in front 
of the lens with its base facing the direction A, 
and a test negative was made. The center cross 
was found to be displaced in the direction A by 
0.239 mm. The principal point, as determined by 
the method described in this paper, was found to 
coincide still with the center of collimation within 

0.002 mm. The point of symmetry was found 
to be displaced in the direction A by 0.496 mm 
the computation being based on targets separated 
37.5° from the axis and the 0° target following the 
method described by Sewell). The values herein 
given are far greater than one likely to be found in 
practice, as the angle of the prism used was 
approximately 12 min. 

It is probable that in lens-cone combinations 
where little asymmetric distortion exists there is 
little to be said in favor of any one of these con- 
cepts over the others. When appreciable asym- 
metric distortion is present, the authors are not in 
a position to state which of these points ought to 
be used in photo-interpretatioa. This must be 
said in favor of the principal point as located by 
this Bureau: It is an invariant point for a given 
lens-cone combination; its location is not appre- 
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ciably affected by placing filters of differing pris- 
matic power in front of the lens; it is possible that 
our present procedure of reporting, however, 
should be amended to include magnitude of the 
prism deviation and its direction for a given lens- 
cone-filter combination so that the center cross 
can be located by those users who prefer it in 
their interpretation processes. 

In the opinion of the authors, there are three 
courses that may be followed to eliminate these 
difficulties. The first course is to set an upper 
limit to the amount of prism deviation permissible 
on a lens-cone combination, possibly 0.015 mm 
for the axial ray. A lens-cone combination that 
has excessive prism deviation would then be un- 
acceptable unless it were reworked to reduce the 
effect. 
denced by the discussion contained in a paper [11] 
by J. V. Sharp and H. H. Hayes. 


The second course of action is to make accurate 


That such a course is practicable is evi- 


measurements of the prism effect, either on the 
lens as mounted in its camera cone or on the lens 
alone. Havingaccurate knowledge of the magnitude 
and direction of the prism deviation, it should then 
be possible to neutralize the effect by using a filter 
on the front of the lens, which is itself a thin prism 
instead of the usual plane parallel. Other re- 
searches in progress at this laboratory indicate that 
such a course is feasible. It is possible that small 
residual effects may still remain, but it is believed 
that such residual effects will be negligible com- 
pared to the asymmetric distortions known to be 
produced by the prism effect. 

The third course of action resembles the first, ex- 
cept that an upper limit on “tangential distortion” 
rather than prism deviation is set. This is pre- 
ferred by J. V. Sharp, because it allows for the 
possibility ° of the existence of a form of “tan- 
gential distortion” that deviates from the pattern 
predicted by the hypothesis that tangential dis- 
tortion is produced by prism effect. Even in this 
case it is possible that a compromise involving 
partial neutralization of the prism effect may be 


advantageous. 


The authors express their appreciation to other 
members of the staff of the National Bureau of 
Standards for assistance rendered during the vari- 


This possibility was brought to the senior author's attention in the courss 


of an informal discussion. 
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ous phases of development of this equipment. 
Robert E. Ward of the Instrument Shop did a 
major portion of the machine work and in par- 
ticular performed the highly accurate boring oper- 
ation on the casting that holds the collimators in 
alinement. George L. Buzas constructed the 
target holders and the devices used in calibrating 
the instrument. Edgar C. Watts assisted in the 
design of testing devices and made the drawings 
used in this paper. William P. Tayman made 
the reticles used as targets in the collimators and 
is the principal operator of the completed instru- 
ment. Arthur A. Magill, Roland V. Shack, Irving 
Malitsky, and William P. Tayman assisted in the 
Wilbur W. Bran- 


calibration of the instrument. 


non made the biprism used for calibration of the 


instrument. 
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A Permeameter for Magnetic Testing at Magnetizing 
Forces up to 300 Oersteds 


By Raymond L. Sanford and Philip H. Winter 


A permeameter is described for magnetic testing at values of magnetizing force up to 


300 oersteds. 


and 1 em thick, with a preferred length of 28 em. 


It is designed to test specimens of rectangular cross section up to 3 em wide 


The instrument requires only a single 


specimen and is more simple and rapid to operate than the Burrows permeameter, which for 


many years has been generally accepted as the standard instrument for magnetic testing in 


this range 


dimensions so that it does not require calibration by any other permeameter. 


It is an absolute instrument, in that its constants are derived from its own 


It is estimated 


that the values of magnetizing force and induction, which are obtained, are accurate within 


| percent, except that for very low values of magnetizing force the precision is the deter- 


mining factor, and the uncertainty may be of the order of 0.05 oersted. 


I. Introduction 


For many years the Burrows Compensated 
Double-Yoke permeameter, either in its original 
form,' or as modified by Gokhale,? has been 


generally accepted as a standard instrument for 


the magnetic testing of straight bars. As it is an 
absolute instrument in the sense that its constants 
are derived from its own dimensions, it does not 
require calibration by reference to any other 
permeameter. For this reason, it has usually 
been employed for the calibration of test bars to 
be used as reference standards for calibrating and 
checking the performance of other permeameters. 
However, the Burrows permeameter has certain 
disadvantages. It is somewhat complicated and 
tedious to operate and requires duplicate speci- 
mens, one of which is used as an auxiliary. Fur- 
thermore, it is extremely sensitive to the effect of 
magnetic nonuniformity along the length of the 
specimen which leads to errors which are indeter- 
minate and may be large. Consequently, its use 
in recent years has been limited in general to the 
testing of magnetically uniform bars to be used 
as standards for calibration purposes. 
a need for an 
C. W. Burrows, The determination of the magnetic induction in straight 


rs, Bul. BS 6, 31 (1909) S117 


S. Patent 1,559,085 (1925). 
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absolute permeameter that would have an accu- 
racy comparable with that of the Burrows instru- 
but which would be less complicated to 
For this reason the development of such 
The result is the 


ment 
operate. 
an instrument was undertaken. 
MH permeameter,’ designed to operate at mag- 
netizing forces up to 300 oersteds. It is capable 
of testing specimens of rectangular cross section 
up to 3 cm wide and 1 em thick, with a preferred 
length of 28 em, although specimens 24 em in 
length can be tested satisfactorily. It is absolute 
in that its constants can be derived from its own 
dimensions, so that it does not require calibration 
by any other permeameter. It is estimated that 
values of either magnetizing force or 
percent, except 


induction 
can be determined to within 1 
that for very low values of magnetizing force, the 
uncertainty due to lack of precision may be of the 
order of 0.05 oersted. It is simple and rapid to 
operate and requires only a single specimen. 


II. Description of the Permeameter 


The permeameter is shown in figure 1. The 
specimen, of rectangular cross section, is sur- 
rounded by the main magnetizing coil. The 
magnetic circuit is completed by symmetrical 


? The designation MH stands for medium H (magnetizing force). 
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U-shaped yokes that are supported by end pieces, 
which in turn are supported by four Bakelite 
blocks. These end pieces also support the speci- 
men and provide for close contact between the 
yoke and the specimen. The symmetrical yoke 
structure is of advantage because it promotes a 
uniform cross-sectional distribution of magnetic 
flux in the specimen. This is especially important 
in determining points on the hysteresis loop where 
the magnetizing force and the induction are in 
opposite directions. Auxiliary magnetizing coils 
are located around the ends of the yokes and sup- 
porting end pieces that make contact with the 
specimen. These auxiliary coils are of advantage 
because they promote a uniform longitudinal dis- 
tribution of flux in the specimen by overcoming the 
contact reluctance between the specimen and the 
voke The flux distribution in the region adjacent 
to the center section of the bar is thereby improved, 
and since the value of the magnetizing force is 
obtained by measurements in this region, the 
accuracy of the measurements depends upon the 
degree of uniformity that can be obtained. 
Values of flux-density, B, are obtained by means 
of a test coil surrounding the bar at the middle. 
An auxiliary test coil whose area-turn value is 
adjusted to be equal to that of the B-coil is located 
adjacent to the B-coil, but not surrounding the 


specimen, and connected in series opposing the 


B-coil. Thus the galvanometer deflections ob- 
served are a measure of the intrinsic flux density, 
B,, rather than the total flux density. This 
obviates the necessity of calculating the air-flux 
correction for bars of different cross-sectional area, 

Values of magnetizing force, HT, are obtained by 
means of //-coils so mounted that they can be 


made to rotate 180° in the region adjacent to the 


Ficure 1. The MH permeameter 


center section of the bar. There are two //-coils 
mounted one above the other on a rotating brass 
table and adjusted to equal area-turn values. 
When these coils are connected in series, opposing, 
with their axes parallel to the specimen, and ro- 
tated through 180°, the galvanometer deflection 
observed is proportional to the gradient of the 
field with respect to the distance from the surface 
of the specimen, which will be referred to as the 
radial gradient. The value of // as measured by 
the deflection when only the upper coil is connected 
to the galvanometer is extrapolated to the surface 
of the specimen by using the observed value of the 
radial gradient. 

The use of rotatable //-coils, or 
has been found to have certain definite advantages 
Identity in 


“flip-coils,”’ 


over the use of stationary //-coils. 
the observed values of magnetizing force for both 
directions of the magnetizing current is a good 
criterion of a satisfactorily demagnetized condi- 
tion. Also, values of magnetizing force for points 
on the hysteresis loop can be obtained by direct 
observation instead of by calculation in terms 
of the difference between two large quantities. 
This improves both the precision and the accuracy 
of the measurements. 

The magnetic circuit is shown in detail in figure 
2. The laminated yokes are of No. 16 gage silicon 
sheet of armature grade and are 4 by 5 cm in cross 
section. The specimen rests on end pieces that 
extend between the supporting Bakelite blocks. 
This provides a good magnetic contact between 
the yokes and the flat surface of the specimen 
The ends of the U-shaped sections 
over and slide on the end 


at the ends. 
are recessed to fit 
pieces so that contact is made at the edges of the 
specimen. This construction provides good con- 
tact between the yokes and the specimen without 
the necessity of clamping. This is an important 
feature, especially when strain-sensitive materials 
are to be tested. A supporting platform asso- 
ciated with the test-coil assembly prevents strains 
that might result from sagging of the specimen. 
The over-all length of the yokes is 28 em, and the 
specimen is preferably of this length. However, 
bars as short as 24 cm can be tested if necessary. 

The main magnetizing coil is 13.8 em long, 
wound on a Bakelite form whose inside diameter 
The winding is of No. 18 AWG Formex- 
insulated magnet wire and has 1,540 turns uni- 
The four compensat- 


is 7 em. 
formly wound in 12 layers. 
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cuit of the MH permeameter. 


FIGURE 2 Vagnet 


ing coils on the yokes are 5.8 em long, and each 
has 430 turns of No 
lavers. As it was desired to 


18 Formex-insulated wire 
wound in eight 
avoid the complication of having to adjust the 
compensating current for each point to be deter- 
mined, the effect of the compensating coils on 
the degree of uniformity of the field adjacent 
to the specimen was determined for various ratios 
of the number of turns in the main and compen- 


Although 


could not be obtained over the whole range for a 


sating coils perfect compensation 
single ratio of turns, that ratio which appeared 
to give the best results on the average was adopted 
and the coils were connected in series with the 
main magnetizing coil. Observations of — the 
radial gradient were then made for different types 
of specimen as related to the value of magnetizing 
foree. It was found that, within the experimental 
error, the percentage correction to the magnetizing 
force varies linearly with the permeability of the 
specimen, as shown in figure 3. It was also 
found that, although under some conditions the 


gradient was not quite linear, the variation from 
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Figure 3 Relation of H-correction to permeability 
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linearity was so small that no significant error is 
introduced from this cause. 

The test-coil assembly is shown in figure 4. It 
is held in position within the main magnetizing 
coil form by Bakelite rings. The platform extend- 
ing between these rings supports the specimen and 
is recessed at the middle so that the upper /7/ coil 
can be located close to the surface of the specimen. 
The two H coils are mounted below the specimen 
on a small table that is rotated by a horizontal rod 
through a rack and pinion gear by an iron-clad 
solenoid and plunger located far enough away 
from the magnetic circuit to avoid interference 
Adjustable stops limit the rotation of these “‘flip- 
coils” to 180 
value of approximately 


The upper coil has an area-turn 
8,250 The 


lower coil has a slightly greater area-turn value and 


em*-turns. 


is shunted with a resistance so as to make its 
effective value equal to that of the upper coil 
Each coil is approximately 2.9 cm long, 4.44 em 
wide, and 1.11 em thick. The axes of the coils 
are 1.11 em apart, and the axis of the upper coil is 
0.87 cm from the surface of the specimen. The 
windings are of No. 34 Formex-insulated magnet 
wire wound on Micarta forms. 

The B-coil, of the same length as the //-coil, is 
wound on a Micarta 
It has 100 turns of No. 31 single silk 


form and surrounds the 
specimen. 
and enamel-insulated wire. The inside dimen- 
sions are such as to accommodate a specimen 3 em 
wide and 1 em thick. The compensating coil of 
the same area-turn value as the B-coil is mounted 
above it and is connected in series, opposing. 

The usual switches and rheostats are employed 
to control the magnetizing current supplied by a 
storage battery. The measurements are made 
with a heavily overdamped ballistic galvanometer 
calibrated by means of a standard mutual inductor 


in the regular way.‘ 
Ill. Performance 


The use of a symmetrical yoke structure 
promotes uniformity of flux distribution, particu- 
larly across the section of the specimen. However, 
magnetic leakage due to the presence of joints in 
the magnetic circuit causes a certain lack of 
uniformity of the field near the surface of the 
specimen in the region occupied by the //-coils 
Many tests were made to determine whether or 








* NBS Circular C456, p. 23 (1946 





FiGuRB 4 


not properly located auxiliary magnetizing coils 
would be effective in improving the distribution. 
From a survey of the leakage flux around the 
magnetic circuit it was determined that auxiliary 
magnetizing coils placed on the end sections of 
the vokes would greatly improve the flux distri- 
bution in the region where the measurements were 
to be made. Distribution tests were then made 
with various specimens with and without the 
auxiliary coils. It 
auxiliary coils the longitudinal distribution of the 


was found that with the 
flux in the region adjacent to the center section of 
the specimen was uniform and, with the use of a 
double //-coil to determine the radial gradient, the 
value of magnetizing force at the surface of the 
specimen could be determined with satisfactory 
accuracy. 

It is generally recognized that the degree of 
magnetic uniformity of the specimen affects the 
results of a nature and 
magnitude of the effect depends not only upon the 


magnetic test. The 


nature and magnitude of the nonuniformity but 
also upon the characteristics of the testing appara- 
tus. In the M// permeameter, the test is limited 
to the middle 3 cm of the bar. It is conceivable 
of course that extreme nonuniformity in the part 
of the bar near the middle might affect the ob- 
served value of magnetizing force unduly, but such 
a condition is relatively rare. However, valid 
comparisons between permeameters that include 
different 
cannot be made unless the test specimens are 


lengths of the specimen in the test 
uniform. 
A source of uncertainty often ignored, but which 


should not be disregarded, is the effect of mechani- 


20 


T+ st-coil assembl yf. 


It is well known that the high-per- 
form, are 
How- 


cal strain. 
meability alloys, especially in sheet 
likely to be exceedingly sensitive to strain. 


ever, the effect of strain in solid materials of other 


types is often overlooked. For instance, in a 
recent test of a solid bar that was not quite 
straight, the induction corresponding to a given 
value of magnetizing force could be changed by 
as much as 2 kilogausses by varying the clamping 
To avoid this uncertainty, the MH 
permeameter was so designed that good contact 


pressure. 


between the yokes and the specimen could be 
obtained without the necessity of clamping. 

In estimating the accuracy of any permeameter, 
it is necessary to take into account not only the 
sources of error peculiar to the individual appara- 
tus, but also those inherent in ballistic methods 
in general. The factors inherent in_ ballistic 
methods include the value of the mutual inductor, 
the area-turns of the test coils, the adjustment of 
the galvanometer sensitivity, ability of the gal- 
vanometer to integrate the value of a delayed 
impulse correctly, temperature effects, degree of 
uniformity of the specimen, and observation errors. 
The principal points peculiar to the M// per- 
meameter are the distribution of the magnetic 
field in the space adjacent to the specimen occu- 
pied by the H-coils and the longivudinal and 
transverse distribution of the magnetic flux within 
the specimen. 

The value of the standard mutual inductor has 
been determined to within +0.05 percent by 
comparison with an absolute standard of high 
accuracy. 

The area-turns of the //-coil have been deter- 
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nined by a method previously described,’ with an 
stimated accuracy of +0.2 percent. 

The current in the mutual inductor for calibrat- 
w the galvanometer is easily adjusted to within 

0.05 percent by using a potentiometer and a 
tandard shunt. 

Recent theoretical and experimental studies of 
he behavior of the ballistic galvanometer in use 
vith this permeameter have shown that under the 
lost severe conditions likely to be encountered, 
he galvanometer integrates correctly to within 
0.1 pereent 

In view of the fact that the specimen is not 
closely surrounded by the magnetizing coil and 
that the currents are relatively low, the heating is 
negligible under normal conditions. 

With the understanding that the test results 
relate only to the middle 3 em of the specimen, 

rors due to nonuniformity along the length of 
the specimen are usually of negligible magnitude. 
The observational error of a single reading may 


be as much as 0.2 mm. 


However, by taking a 
sufficient number of readings, the probable error 


With suitable 
adjustment of the galvanometer sensitivity, this 


can be reduced to about 0.1 mm. 


will correspond to not over 0.2 percent for values 
of magnetizing force of 5 oersteds or over, but for 
values of less than 5 oersteds the precision is the 
determining factor and the uncertainty is of the 
order of 0.05 oersted The observational error 
in the measurement of the induction can easily be 
kept within 0.2 percent by suitable adjustment of 
the galvanometer sensitivity. 

The accuracy of the observed values of mag- 
netizing force is limited mainly by the degree of 
uniformity of the field in the space occupied by 
the //-coils and adjacent to the specimen. In 
order to obtain a reasonable degree of precision, 
it Was necessary to choose dimensions for the test 
coils that represent the best compromise between 
the conditions for sensitivity and for accuracy. 
From the results of the experimental exploration 
of the field, it is estimated that errors due to non- 
iniform distribution do not exceed 0.3 percent. 

Experimental evidence indicates that the longi- 
tudinal and transverse distributions of flux within 
he part of the specimen covered by the B-coil 
ure sufficiently uniform to make errors from this 

urce negligible. 

For purposes of comparison, normal induction 
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FIGURI 5 Vorma induction and h ysteresis curves 10 
tungsten-steel ba 
Curves represent data obtained with solenoid. Check points represent 
data obtained with MH permeameter 


and hysteresis tests were made on a tungsten steel 
bar about 94 em long with a rectangular cross- 
sectional area of about 1.8 em*. After confirming 
the magnetic uniformity of the bar, it was first 
tested in a long solenoid and then in the M// 
permeameter sefore testing the bar in the 
solenoid, the distribution of flux density in the 
center region of the solenoid was checked, and it 
was found that the radial gradient of /7 was very 
small, so that values of 7/7 could be obtained with 
a single flat //-coil with an error due to this gra- 
dient of less than 0.2 percent. 

The curves in figure 5 are from values obtained 
with the solenoid, and the check points are plotted 
from values determined with the M// perme- 
ameter 
points shows that the discrepancy between the 


values determined with the solenoid and MH 


A comparison of the curves and check 


permeameter is less than the observational error 
and well within the normal tolerance. 
Considering all the probable sources of error, it 
is estimated that by the exercise of care in taking 
the readings and averaging several readings for 
each point, values of induction can be obtained 
that will be accurate to within about 0.5 percent, 
and values of magnetizing force can be obtained 
that will be accurate to within 0.5 percent, or 0.05 
oersted, whichever is larger. For routine tests, 
values of either induction or magnetizing force 
accurate to within | percent should be obtained 


without difficulty. 


WasHINGTON, March 9, 1950. 
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Heat Capacity of Sodium Between 0° and 900° C, the Triple Point and Heat of Fusion 


By Defoe C. Ginnings, Thomas B. Douglas, and Anne F. Ball 


Using an improved ice calorimeter and furnace, the enthalpy changes of two samples 


of pure sodium have been accurately measured by a drop method at a number of tempera- 


tures between 0° and 900° C. 
and entropy, 


both the solid and liquid. 


based on zero values at 0° C, 


Equations are derived to fit the data, and values of enthalpy 
as well as the heat capacity are tabulated for 


Sources of significant experimental error are examined critically, 


and some theoretical implications of the results are discussed qualitatively. 


I. Introduction 


liquid sodium for heat 
transfer at high temperatures have been known 
for some time. Not 


values of the heat capacity of sodium been limited 


The ady antages of 


only have experimental 
to the temperature range below 300° C, but these 
values have been widely discordant at 200° to 
00°. It was the purpose of this investigation to 
measure the heat capacity of pure liquid sodium 
with higher accuracy and at higher temperatures 
than had been done previously. 


II. Samples 


The measurements were made on two samples 
of sodium enclosed in stainless steel containers of 
the type shown in figure 1. The stainless steel 
used was chosen for its resistance to attack by 
sodium at high temperatures. The sodium had 
been distilled once in glass in vacuum at or below 
00° C. The weights of the samples were 5.3404 
and 6.0297 g, respectively, and each was ad- 
mitted in vacuum to a stainless steel cylinder 
weighing 16.840 ¢ and having a capacity of 9.7 
m’. After admitting helium to a pressure of one- 
thirtieth atm, each cylinder was sealed by welding 
had subsequently tested for leakage by noting 
vhether exposure to 850° C for 15 min in an 
vacuated vessel caused any change in weight. 
lhe welding process utilized pulsed inductive 
eating {1].! 


ossible to seal the container without changing 


With this welding technique, it was 


ts weight by more than a milligram. 





res in brackets indicate the rature reference it the end of this 


Specific Heat of Sodium 


The sodium sample used for most of the measure- 
ments was analyzed spectrochemically at this 
Bureau after the heat-capacity experiments. Of 
the 47 tested for, 
common metals, only three were detected. 


including all the 
_Cal- 


elements 








_—_— 


ad 5/8 oa | 


Figure 1. Scale diagram of sample container. 
cium and lithium were found in amounts estimated 
0.0001 and 0.001 
whereas potassium was estimated to be between 
0.001 and 0.01 percent. 
pletely soluble in water. 
The spectrochemical analysis did not include a 


at between weight percent, 


The sample was com- 
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Ficure 2. Schematic diagram of furnace and ice calorimeter. 


A, Calorimeter well; B, 


container; FE, ice bath; F, 


24 


beaker of mercury; C, 
copper vanes; @, gate; 


glass capillary; D, sample 
HH, platinum heater; J, ice 


search for the common nonmetals. However 
other samples from the same original source anc 
distilled by the same procedure were reported t 
contain as much as 0.01 percent by weight o 
oxygen (as Na,O). The typical content of hydro 
gen before distillation amounted to 0.006 percent 
by weight. A melting curve observed by the 
authors on the principal sample is described in this 
paper. This indicated that the impurity in 
solution in the liquid but not in the solid amounted 
to about 0.0016 mole percent at the triple point. 


III. Method and Apparatus 


The method as applied to an earlier apparatus 
has already been described [2, 3, 4]. In brief, the 
method is as follows. The sample, sealed in its 
container, is suspended in a furnace until it comes 
to constant temperature. It is then dropped into 
a Bunsen ice calorimeter, which measures the heat 
evolved by the sample plus container in cooling to 
0° C. A similar experiment is made with the 
empty container at the same temperature. The 
difference of the two values of heat is a measure of 
the change in enthalpy of the sample between 0° C 
and the temperature in the furnace. From 
enthalpy values of the sample so determined for a 
series of temperatures, the heat capacity can be 
derived. 

1. The Ice Calorimeter 

The construction, calibration, and use of an 
earlier Bunsen ice calorimeter in this Bureau have 
been described [2]. The present model, shown 
together with the furnace in figure 2, incorporates 
several improvements. In the first place, the 
mercury-water interface is now located in the 
bottom of the inner glass calorimeter vessel. This 
has two advantages. First, the area of the inter- 
face is much larger than before, so that for a given 
influx of heat, the level of mercury in the calorim- 
eter changes very little. The calorimeter and its 
contents are slightly compressible, so that a change 
in pressure in the calorimeter results in a change in 
volume that must be distinguished from the change 
in volume due to heat input. With the earlier 
calorimeter, the change in volume produced by the 
change in head of mercury affected the calibration 
factor by about 0.01 percent, but with the present 
mantel; J, K, L, silver cylinders; M, mercury; N, Inconel tubes; P, Pyrex 
containers; R, mercury reservoir; S, platinum shields; 7, mercury “‘temper- 


ing” coil; V, needle valve; W, water; Y, porcelain spacers; 1, 2, 3, 4, thermo 


couple junctions. 
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ilorimeter, the effect of the change in head is only 
bout 0.004 percent. The electrical calibration 
xperiments automatically compensate for this 
ffect in heat measurements. However, reduction 
if the magnitude of the effect permits determina- 
on of the calibration factor of an “ideal” ice 
alorimeter (one that operates at constant head) 
vith greater accuracy. A second reason for locat- 
ng the mercury-water interface in the bottom of 
the calorimeter is to avoid danger of breaking the 
inner glass vessel. The inlet tube, now being 
filled with mercury rather than water, is not 
subject to blocking by freezing, an event that 
would burst the glass vessel. The present inlet 
tube includes a coil (7'), for bringing the mercury 
to the temperature of 0° C before it enters the 
calorimeter. 

The present calorimeter has been equipped 
with a tin-plated system of copper vanes, as 
pictured in figure 2. This results in more rapid 
distribution of heat from the well to the ice, so 
that the lag in reaching thermal equilibrium is 
now caused almost wholly by the thermal resis- 
tance between the heat source (sample in con- 
The latter 
lag has been considerably reduced by using a slow 


tainer) and the inner wall of the well. 


stream of dry helium gas up the well instead of 
carbon dioxide gas (as shown in fig. 1 

Unlike the earlier calorimeter, the new one has 
Cali- 
brations were made with a removable heater pro- 


no built-in heater for electrical calibration. 


vided with two zones of thermal contact with the 
calorimeter well for trapping the heat that would 
otherwise escape through the leads. This heater 
has a resistance of about 100 ohms instead of the 
10 ohms in the earlier calorimeter. This reduces 
the uncertainty of accounting for the heat de- 
veloped in the leads between the calorimeter and 
the jacket to a few thousandths of a percent of 
the total heat input to the calorimeter. 

About one hundred electrical calibration experi- 
ments were made with the new calorimeter, with 
These 


experiments gave a calibration factor of 270.46+ 


nergies varying from 12,000 to 25,000 j. 
1.03 abs j/g of mercury, which agrees with the 
value of 270.41+0.06 obtained for the earlier 
alorimeter within the uncertainty assigned to the 
atter figure. About a hundred additional cali- 
ration experiments were made with low energies 
about 100 j) in a search for small absolute errors 
that might not be detectable in high-energy 


Specific Heat of Sodium 


experiments. These experiments indicated that 
the calorimeter failed to register about 0.4 j, an 
amount of heat that is negligible in most measure- 
ments. 


2. The Furnace 


The furnace is shown in position over the ice 
calorimeter in figure 2. It was mounted on a 
ball bearing so that it could be swung aside for 
access to the calorimeter. This furnace incorpor- 
ated certain improvements that reduced temper- 
ature gradients in its central region and permitted 
more accurate temperature measurements over 
most of its temperature range. 

In order to minimize end effects, the furnace 
was made 24 in. high, instead of 18 in. as with the 
earlier one [3]. The furnace was electrically heat- 
ed by means of No. 26 platinum wire (/7) wound 
on a grooved alundum tube having an inside 
diameter of 2 in. This heater was made in 
three separate sections corresponding in elevation 
to the three silver evlinders, which were located 
inside the alundum at J, A, and LZ. Silver was 
used because of its high thermal conductivity. 
By maintaining the silver cylinders J and L at 
the same temperature as the evlinder K, the tem- 
perature gradient in A could be made negligible. 
The silver cylinders were supported by porcelain 
spacers ()). The silver and porcelain cylinders 
were lined with an Inconel tube (14% Cr, 80% 
Ni, 6% Fe) having an inside diameter of 1 in. in 
the central and lower parts of the furnace but with 
a diameter of one-eighth in. in the upper part. 
This tube served to enclose the sample container 
and its suspension wire (No. 32 Nichrome \V), 
so that an inert atmosphere could be maintained 
at high temperatures where the stainless steel 
container would oxidize in air. Helium was used 
for the atmosphere in order to minimize the time 
necessary to hold the capsule in the furnace. 

Figure 2 shows some of the vertical holes drilled 
through the silver and porcelain. These were 
placed 90 deg apart azimuthally. Three of them 
(one-fourth-in. in diameter) extended through the 
length of the furnace while two (three-sixteenth-in. 
in diameter) extended only to the top and bottom 
respectively of the center silver cylinder as shown 
in figure 2. These holes were fitted with thin- 
walled Inconel tubes, one of which served to hold 
the platinum-rhodium thermocouple sometimes 
used for temperature measurement of the central 
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silver section (K). In the two %¢-in tubes were 
placed differential thermocouples used to keep the 
end silver sections (J, LZ) at the same temperature 
All thermocouples were 
insulated in four-hole porcelain tubes, which fit 
A platinum resist- 
ance thermometer was encased in a_platinum- 
rhodium tube, which also fit snugly in one of the 
one-fourth-in. Inconel tubes. 


as the central section. 


snugly in the Inconel tubes. 


In another one-fourth-in. tube were placed 
three small auxiliary heaters, located at the eleva- 
tions of the three silver cvlinders. It was possible 
with these heaters to add heat at a low rate to any 
silver cylinder for fine-temperature regulation. 
It was possible to maintain the central silver 
cylinder constant to 0.01 deg and the end silver 
evlinders within a few tenths of a degree of the 
central silver. These auxiliary heaters were very 
useful in fine-temperature regulation since their 
heat was mostly delivered to the silver in 30 sec, 
thus avoiding the larger and troublesome lag that 
would have been encountered had the fine regula- 
tion been dependent on adjusting the currents in 
the main heaters. 

In addition to the thermocouple mentioned, 
several more were installed just outside the main 
heaters (/7). These 
measure the temperature of this region and were 


thermocouples served to 
used in a safety device to turn off automatically all 
current in the furnace if its temperature ap- 
proached 961° C (the melting point of silver). 
The suspension of the container in the furnace 
and its drop into the calorimeter is similar to that 
described earlier [3]. The braking started after 


the container entered the calorimeter. The 
weight of the falling system was kept constant in 
all experiments. Two platinum shields were used 
above the container in the same manner as for- 
merly except that the diameters of the shields were 
increased from fi-e-erghths in. to three-fourths in. 


since the diameter of the calorimeter well was 


increased from three-fourths in. to seven-eighths 


In. 


3. Measurement of the Temperature of the 
Capsule in the Furnace 


Up to 600° C, a strain-free platinum resistance 
thermometer was used to measure the temperature 
of the central silver cylinder that surrounded the 
This thermometer was calibrated at 
Temperature 


capsule. 
this Bureau on the International 
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Seale [5] at the ice, steam, and sulfur points, and 
its resistance was always checked at the ice point 
after experiments at 600° C. Between 600° and 
900°, the resistance thermometer was removed 
from the furnace, and the platinum-rhodium 
thermocouple was used. This thermocouple also 
was calibrated at this Bureau and was checked 
against a standard thermocouple before and after 
the measurements. 
ance thermometer with the thermocouple showed 
agreement within 0.2 deg C in the temperature 
Both the resistance thermom- 


A comparison of the resist- 


range up to 600°C. 
eter and thermocouple were made to have good 
thermal contact with the silver, so that their 
indications would not be affected by heat conduc- 
tion along their leads. Experimental verification 
of this was found in tests with varying depths of 
immersion. The resistance thermometer winding 
and the principal junction of the thermocouple 
were both located at approximately the height of 
the sample container in the furnace. 

It has been assumed that the sample container 
attains the temperature of the central silver, 
which surrounds it except at the bottom.  Al- 
though calculations showed this assumption to be 
justified, experiments were made to test this. 
These experiments consisted of first making regu- 
lar drop experiments at 700° C with the end silver 
pieces (J and ZL) regulated in the normal manner 
within a few tenths of a degree of the central 
Then another series of drops was 


silver (A). 
made, regulating both end silver pieces about 50 
deg colder than the center. In other words, the 
temperature difference between the ends and the 
central silver was exaggerated by a factor of over 
one hundred. Even under these extreme condi- 
tions, the experiments gave a result only 0.1 
percent lower. This indicated two main facts. 
First, with the exaggerated temperature difference, 
either the temperature gradient in the central 
silver was small, or the thermocouple was so 
placed that it represented the temperature of the 
Actu- 
ally calculations show that the gradient is negligi- 


Second, this also indicated 


sample container in spite of the gradient. 


ble in normal use. 
that any lowering of the temperature of the sample 
container below that of the main silver cylinder K 
due to radiation loss from the sample container 
out the bottom of the furnace is negligible. 
solid angle had 


Calculations considering the 


already shown that this should be true. 
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Since the temperature of the sample container 
is not directly measured, it was necessary to 
low sufficient time for the container to reach 
1c temperature of the silver. Special tests have 

en made to determine this time, similar to those 
escribed in detail in an earlier publication [3]. 
n the present experiments using helium in the 

nace, the container with sodium acquired the 
emperature of the furnace (within 0.01 deg 
vithin 15 min, whereas the empty container took 
ess time in proportion to its heat capacity. In 
the experiments with sodium just above its melting 
point, extra time was taken to allow for absorp- 
tion of the heat of fusion. 

It was also necessary to consider the possibility 
of azimuthal temperature gradients in the silver 
due to the periodic addition of small amounts of 
heat by means of the auxiliary heaters that are 
located in one of the vertical holes in the silver. 
Calculations showed that no significant tempera- 
ture gradients resulted with the relatively small 
amounts of heat supplied in this way. 


4. Errors Due to Cooling of the Sample Container 
During Its Drop into the Calorimeter 


Some users of the “drop’’ method in the past 
have simply subtracted the calculated enthalpy of 
the sample container from the measured enthalpy 
of the container plus sample. Of course, this 
procedure does not make allowance for heat lost 
during the drop. By making drop experiments 
with the empty container, most of this heat lost 
is accounted for, since it is approximately the same 
in both types of experiments. However, for the 
most accurate measurements, it is desirable to 
test the validity of this assumption that the heat 
loss is essentially the same for both empty and 
filled containers. The heat loss during the drop 
Is dependent on the temperature difference be- 
tween the container surface and its surroundings. 
This temperature difference will be slightly greater 
for the filled container, since its total heat capacity 
s higher than that of the empty container. This 
difference between the amounts of heat lost by the 
mpty and the filled containers has been estimated 
by the authors for a sample assumed to have 
infinite thermal conductivity and perfect contact 
with the container. This gives the maximum 
possible difference The radiation loss was esti- 
mated on the basis of a container with an area of 
‘5 em? with an emissivity of approximately unity. 
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The heat lost by radiation from a container at 
900° C 
temperature for one-quarter sec is about 80 j. 


dropping through a region near room 


The convective heat loss is much more difficult to 
estimate. Using the formula given in the Inter- 
national Critical Tables [6] for horizontal cylinders 
as being somewhat analogous, the convective loss 
in helium was calculated as about 80 j. Using 
these figures with the known heat capacities of 
the empty and filled containers, it is estimated 
that there would be a maximum error of about 4 j 
in the value for the enthalpy of the sample at 900° 
C. This is equivalent to about 0.06 percent on 
enthalpy, and a somewhat greater error on heat 
capacity. 

Due to the approximate nature of the calceula- 
tions of this error, some experiments were made 
to test if the heat losses calculated in this w av were 
reasonable. Eight experiments were made at 
300° C, dropping another container having about 
the same emissivity. In these experiments, the 
time of fall of the container was varied by factors 
of two and four greater than the normal time of 
fall (almost free fall). Assuming that the heat loss 
is directly proportional to the time of fall, it was 
estimated that the total heat lost in a normal fall 
at 300° C was about 6 j. This figure is consider- 
ably less than the value of about 16 j calculated in 
the manner described above (12 j by convection 
and 4 j by radiation). This would indicate that 
the heat loss is less than the calculations indicate 
and that at 900° C 
be 0.02 to 0.03 percent. 


, the error due to this loss may 
Since this error is small 
and is not known accurately, no attempt was made 


to apply a correction. 


IV. Results of Enthalpy Measurements 


A total of 129 measurements of enthalpy was 
made, 67 with the two empty containers and 62 
with the two containers containing sodium. The 
average deviation from the mean was +0.02 per- 
cent with a given sample, and the difference 
between the values of enthalpy for the two samples 
did not exceed 0.3 j/g or 0.03 percent, whichever is 
larger. The detailed results of individual runs 
are given in table 1. Measurements were often 
repeated at lower temperatures after experiments 
at higher temperatures without disclosing any 
evidence that significant permanent changes had 
occurred in the samples meanwhile. 
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sample lated 
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s, 44y 
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&, 452 } 
® International Temperature Seale of 1948 [5] 
b On alternate container of same mass 
© Weighted less because of unsteady heat leak. 
{On alternate container, containing 6.0297 g of sodium. he other con 


tainer had the same mass but contained 5.3404 g of sodium, 


No corrections for impurities were made, since 
these were believed to be much smaller than the 
accidental error. All weights were reduced to 
an in-vacuum basis. Small corrections were made 
to bring the results to a common basis at the listed 
temperatures. The values of enthalpy (/7,—H,) 
given in table 1 are at the vapor pressure (satura- 
tion pressure). Small corrections were made for 
the thin oxide coatings that the containers acquired 
in the furnace, using published heat capacities of 
Fe and Fe,0, [7]. In addition, in order to evaluate 
the heat capacity of the liquid alone, it was neces- 
sary to make small corrections at the highest 
temperatures for the heat given up in condensing 
the sodium vapor. These corrections, based on 
published vapor pressure data and the volumes of 
the vapor, amounted to about 0.03 percent of the 
enthalpy at the highest temperature. The total 
corrections of all kinds did not amount to more 
than 0.1 percent of the enthalpy. 

The experimental values found for the enthalpy 
of sodium (at saturation pressure) in absolute 
joules per gram, are represented by the following 
empirical equations. //, is the enthalpy of the 
solid at 0°C, and t represents degrees C. 
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solid) —H, (solid) = 1.19926t+ (3.247) (10-4) 





(3.510)(10-*®)t (0° to 97.80° C) 


liquid H, (solid 98.973 +-1.436744t 

2 90244 107-*)# (1.54097) (107-")é 
Len 

24000e ‘** (97.80° to 900° C) 


[he last term of the latter equation (liquid) 
ounts for the difference between A, H, and 
net heat measured. As shown by Osborne 
this difference is equal to PV/m, where P is 
vapor pressure of the sodium, V is the volume 
the container, and m is the mass of sodium. In 
temperature investigated, this term 
aches a maximum value of 0.02 percent at 900° C 
These 


juations give 113.2 abs j/g as the heat of fusion 


range 


nd 1s entirely negligible for the solid. 
the triple point, 97.80° C. 
\fter dropping the last term in the equation for 
enthalpy of the liquid, the two preceding 
quations then give the measured heat differences 
curve. Therefore their 
the heat 


ilong the saturation 


capacities 


direct differentiation gives 

(’,) under saturation temperature and pressure, 
s follows: 

( solid) = 1.19926-+- (6.494) (10~*)t 


1.0531) (10-5) # (0° to 97.80° C 


C, (liquid) =1.43674— (5.8049) (10~‘)t 


1.6229) (10-*)t (97.80° to 900° C 


From these equations directly result the following 


quations for the entropy (also at saturation 


, absolute joules g~' deg K™', in excess 


T represents degrees Kelvin 


essure 
that at 0° C. 
aking 0° C=273.16° K 


Sy (solid Sor sex (solid 4.16241 logy)7 


‘a 


5.1036) (10 T+ (5.2656) (10~°) 7?—9.14016 


(0° to 97.80° C 


»pecific Heat of Sodium 


Sr (liquid) - Sor73 16°K (solid) 3.75276 logy T 


(8.3303) (10~*) 7+- (2.3112) (107*) 7°—8.67398 


(97.80° to 900° C) 
Values of enthalpy, heat capacity, and entropy 


calculated from these equations for rounded tem- 


peratures are listed in table 2. 


TABLE 2. Enthalpy, 8 pe cific heat, and entropy of sodium 


remy u H,— Hox ( Sr-—S » 
( 16s deg 
0 0 1. 1991 0 
25 2 1. 2221 10590 
) 61.2 1.2 259 
93. Qt l l W11E 
97.80 (solid 123. 68 l i SSSI] 
[ Melting 113.2 
97.80 (liquid 236. 86 1. 3845 69090 
100 239. WO ] s. OWU0E 
200 375. OF ] 43 1. 0222 
20) 508. 03 ] 42 1. 2757 
400 637. 09 1. 2786 1. 4832 
mM) 764. OF 1. 619 1. 6590 
600 SSO. 8&2 1. 2548 1. 8121 
700 1,01 } 1. 2569 1. 9481 
S00 1, 141. ¢ 1, 268 071 
G00 1, 269. 4 1. 2887 2. 1853 


V. Determination of the Triple Point and 
the Temperature Change During Melting 


In order to measure the triple point and the 
amount of impurities (liquid-soluble, solid-insolu- 
ble) in the sodium, a resistance thermometer (60 
ohms of No. 40 enameled copper wire) was wound 
directly around the lower half of the sample 
container used in most of the enthalpy measure- 
ments. This container was raised into the fur- 
nace, and the copper resistance thermometer was 
then calibrated at 94° 100° C 
platinum resistance thermometer located in the 


and against th: 


silver cylinder (A). Starting with this container 
and the furnace at equilibrium at 94° C, the 
container was heated adiabatically using the cop- 
per resistance thermometer as an electric heater 
and raising the temperature of the furnace at 
about the same rate. After an appropriate time 
interval the electric power was turned off, and 
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the container was allowed to come to a constant 
temperature, which was observed by the copper 
resistance thermometer. This cycle was re- 
peated about 20 times, using electric power of 
one-fourth watt over a total heating period of 90 
min. It was therefore possible to observe the 
relation between the temperature of the sodium 
and the fraction that had melted. 

The temperature and fraction of the sample 
melted were consistent with Raoult’s law within 
the precision of the measurements, and indicated 
the triple point of the sample to be 97.792° C and 
that of pure sodium to be 97.80°. This triple- 
point lowering corresponds to a liquid-soluble, 
solid-insoluble impurity of 0.000016 mole/gram- 
atom of sodium. The freezing point of pure 
sodium at l-atm pressure is calculated to be 
97.81°C. 


VI. Effect of Mechanical State of the 
Solid 


E. Griffiths [9] claimed that he found reproduc- 
ible changes in the heat capacity of crystalline 
sodium amounting to as much as 1.5 percent, 
depending on whether in its previous heat treat- 
ment, the sample had been “quenched” or 
“annealed’’. Similar tests were made by the 
authors on the two sodium samples reported here, 
and with conditions of quenching and annealing 
somewhat more extreme than employed by 
Griffiths. The furnace temperatures used were 
chosen to cover the temperature range where he 
found the greatest difference. The results are 
summarized in chronological order in table 3. 

The first set of experiments (Set 1) was made 
after freezing slowly over an interval of 1 hr and 
then cooling at the rate of about 5 deg/hr. Pre- 
ceding the experiments in Set 2 and on the same 
day, the molten sample was quenched by putting 
its container into a tube immersed in ice. The 
treatments before Sets 3 and 4 were respectively 
similar, except that the annealing carried out on 
the day preceding Set 3 involved freezing the 


sample over a period of 3 hr and then cooling 


the solid to room temperature at a rate varying 
from 3 to 5 deg/hr. Before observing Set 5, the 
sample was remelted and requenched by immersing 
its container directly into a dry ice-acetone bath 
(—78° C). It was then kept at this temperature 
overnight until just before the experiments. 
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TaBLe 3. Experiments to determine the reproducibility 


the mechanical state of solid sodium 


Measured heat 
(Na and Con- 
tainer 
Tem- Mass of 
perature Nain Prior heat treatment 
t sample Indi- 
vidual Mean 
experi- set 
ment 


ahe j 


SOO. 7 


S61. 7 
Annealed 
863.4 


S62. 2 


Quenched 


6. 0297 Quenched 


6. 0297 Requenched rapidly by 


dry ice 


The results tabulated in table 3 indicate that the 
quenching and annealing treatments did not 
change the heat capacity of the sodium by more 
than 0.2 percent and that these changes seem to 
be accidental and within the precision of the 
measurements. 


VII. Reliability 


An index to the reproducibility, or “‘precision”’, 
of the enthalpy measurements is afforded by the 
deviations from the mean, as shown by the results 
of the individual measurements showr in table | 
Another index is provided by the deviations (also 
listed in table 1) from the smoothed values given 
by the empirical equations adopted. All the 
results for liquid sodium lead to a probable error 
of the mean of from 0.01 to 0.03 percent on 
enthalpy (depending on the temperature), and as 
much as 0.3 percent on the derived heat capacity 
values, with a somewhat increased uncertainty be- 
low 200° and above 800° C because of the difficulty 
of determining the derivative of an empirical 
function accurately near its ends. For solid 
sodium, the relative uncertainties become con- 
siderably greater at the low temperatures, owing 
to the inability to determine heats by the ict 
calorimeter with a precision of better than 0.4 j 
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In addition, all factors suspected of contributing 
preciable systematic errors were analyzed to 
termine their likely contributions. Direct com- 
risons between the platinum resistance ther- 
ometer and the platinum-rhodium thermocouple 
» to 600° indicated differences in temperature of 
om 0.0 to 0.2 deg. 


rors due to the very small temperature gradients 


Calculations indicated that 


nside the silver cylinder in the furnace were 
ntirely negligible. As far as uncertainties in 
easuring the temperatures on the International 
lemperature Seale are involved, the resulting 
wobable error? in the enthalpy is thought to be 
mly about 0.02 percent between 100° and 600°, 
ind perhaps three times this much at 900°, where 
was necessary to resort to measurement using a 
hermocouple. As a consequence, the probable 
error in the heat capacity due to uncertainty in 
emperature measurement may be considered to 
‘within 0.03 percent from 100° to 600°, but may 
reach 0.2 percent at 900° C 
In case the sample contained 0.01 percent by 
weight of oxygen as sodium oxide, as discussed 
earlier in this report, the error in enthalpy or heat 
capacity due to its gradual deposition from solu- 
tion as the temperature fell might be as much as 
0.05 percent. A similar error would be present 
f the liquid sodium dissolved appreciable amounts 
of the steel container; but recent measurements on 
the solubility of iron in liquid sodium by L. F. 
Epstein [8], of the Knolls Atomic Power Labora- 
tory, indicate a solubility at 900° C of 22 parts per 
million, which would lead to negligible error in the 
The analysis of the dif- 
ference between the heat lost, during the drop 


present: Measurements 
nto the calorimeter, by the sample container and 
by the empty container, showed that this effect, 
greatest at the highest temperatures, might lead 
to errors of 0.03 percent on enthalpy and 0.1 per- 
cent on heat capacity at 900°. Taking due con- 
sideration of the sources of error mentioned and 
the 
that the probable error of the enthalpy values of 


other less significant ones, authors believe 


sodium reported here may be considered to be 
tween 0.1 and 0.2 percent (except below 100°, 


s mentioned above), whereas the probable error 


robable error’ as used in the remainder of this paper includes the 


estimate of certain systematic errors ro each factor subject to 
there was assigned a figure such that it was believed that the actual 
was as likely greater as smaller than this figure These figures were 
mbined statistically to yield an over-all probable error in the given 


imental quantity 


pecific Heat of Sodium 


SSOH354 oO 


of the heat capacity was similarly estimated to be 
between 0.3 and 0.4 percent. 

One experimental check on the over-all accuracy 
of the apparatus in measuring enthalpy was made 
by measuring the heat delivered to the ice calo- 
rimeter by a Monel capsule containing water and 
dropping from 250° C. By thus determining in 
several measurements the difference in heats for 
two amounts of water differing by about 6 g, a 


250 


value of 1042.05 abs j/g for a)’ of water, de- 


fined elsewhere [10], was obtained. This figure 
differs by only 0.02 percent from the value of 
1041.85 published in the latest report [11] on the 
thermal properties of water as accurately measured 
earlier in this laboratory with an adiabatic 
calorimeter. 

The value found for the triple point is considered 
to have a probable error of 0.03 deg, dependent 
mostly on the uncertainty in the calibration of the 
the 


A probabie error of 0.3. percent 


copper resistance thermometer used for 
measurements. 
is assigned to the heat of fusion This figure is 
based on the estimated reliability of the enthalpy 
values in the neighborhood of the triple point 
the 


melting curve were actually in liquid solution at 


If all the impurity indicated by measured 
94°, then 0.25 percent of the sample was liquid at 
In that case, the value of the 
enthalpy of the solid given in table 1 would be 0.2 


this temperature. 


j/g too high due to this effect, making the cal- 
culated heat of fusion 0.3 j/g too low. 


VIII. Comparison with Results of Other 
Investigations 


A number of other investigators [9, 12, 13, 14, 
15, 16, 


capacity of solid and liquid sodium above 0° C, 


17] have made measurements of the heat 


and most of these are shown for comparison in 
figure 3, where the full-line curves represent the 
values of the present investigation as given by 
the equations presented earlier in this paper. 
(NBS) given for the 
liquid in figure 3 was calculated from two adjacent 


Each experimental point 
observed enthalpies given in table 1. Rengade’s 
[12] measurements, which made use of an ice 
from 15° to 100°, and 
although he stated that his sample was not very 


calorimeter, extended 
pure, his precision amounted to about 0.2 per- 
Griffiths [9] employed a dynamic method 
and covered the range 0° to 140°, claiming a 


cent. 
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Figure 3. Comparison of heat-capacity values 

reproducibility of better than 0.1 percent, except 
for the variation of about 1.5 percent on the solid 
caused by varying the heat treatment, as men- 
tioned Eastman and Rodebush’s [13] 
measurements were mostly at low temperatures, 
litaka [14] 
made determinations of the heat capacity of the 


above. 
but afford a comparison at 20° C. 


liquid between 125° and 290°, the latter being the 
highest temperature attained prior to the present 
investigation. He dropped his samples from a 
furnace into his calorimeter and took due pre- 
caution to avoid certain systematic errors, but 
his precision was not better than several percent, 
and it will be noted from the graph that his 
values for heat capacity are considerably higher 


TABLI } Values for the me lting point and heat of fusion 


ot sodi um 


Heat of fusion 


rhis ins estigation 

Ladenburg and Thiele [1S 

Edmondson and Egerton 

119] 

rammann [18 7 

Griffiths [9 v7 115.35 (113.4, less 
curate 

Bridgman [20] 07. 63 125. 5 

Rengade [12 

Bernini [11 07. 63 74.5 

litaka [14 


Joannis [16 


than those of the other intestigators. Dixo: 
and Rodebush [15] made measurements with 
precision of about 1 percent, from 121° to 178 
Their method consisted in measuring the adiabati 
temperature-pressure coefficient, and seems t 
have been subject to experimental difficulties du 
to the high thermal conductivity of the sodium 

The values reported for the melting point an 
heat of fusion of sodium are given in table 4 
These melting point values in themselves are not 
of course, a reliable criterion of the relative purity 
of the various samples used. 


IX. Discussion of Results 


The values of the heat capacity of sodium 
reported herein have practical importance, inas 
much as they cover the temperature range between 
0° and 300° C with greater accuracy and precision 
than has been the case heretofore, and in addi- 
tion extend the measurements to much higher 
temperatures, where experimental values for 
sodium have been entirely lacking until now. At 
the same time these values should be of consider 
able theoretical interest, since the extensive 
interest in the structure and properties of the 
liquid state in recent years has led to the theoret- 
ical calculation of the heat capacity of certain 
liquids, including sodium, whereas accurate experi- 
mental measurements on liquid metals have been 
relatively very meager. 

The accuracy of the results given in this report 
is sufficient to show clearly the trends of the 
heat capacity curves of sodium with temperature, 
both below and above the triple point. The fall 
of the liquid curve to a minimum, a phenomenon 
in contrast to the steady rise with temperature 
for the great majority of liquids, is by no means 
unique in the case of metals, being paralleled 
capacity of mercury [21]. The 
l 


by the heat 


entropy of fusion of sodium, 7.03 j deg 
I 


gram 
atom is one of the lowest values for metals 
According to current theories, the metal does not 
acquire complete randomness on melting, and 
the increased heat capacity of the liquid over what 
it displays at much higher temperatures is to b 
attributed to a further absorption of ‘communal 
entropy above the melting point. Taking th 
values of the present investigation and extra 
polating linearly to the melting point that pa: 
of the curve (fig. 3) above 700°, it was calculate: 
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t the excess heat capacity between the triple 
nt and 600° is equivalent to an entropy increase 
3 j deg™' gram-atom™', a figure of the right 
er of magnitude to be explained in this way. 
[t has been pointed out further that the increase 
heat capacity of the solid alkali metals with in- 
asing temperature is enhanced in the regions 
temperature just below the respective melting 
nts, and this has been attributed to an acqui- 
tion of a part of the randomness, most of which 
urs in fusion. The equation which Rengade 
2| gave for the heat capacity of solid sodium 


to the melting point, being linear with 


from 15° 
temperature, does not of course show this feature. 


lhe equation found for the solid in this investiga- 

mn. however, does have the sign of curvature 
just referred to, and so in this respect ts at vari- 
ance with an unmodified Debye function above 
the characteristic temperature. 

\side from any theoretical explanation of its 
approximate values, the heat capacity of liquid 
sodium is, as is well known, one of the determinants 
of the curvature of the curve. 
Ladenburg and Thiele [18], in a critical evalua- 


vapor-pressure 


tion of numerous vapor-pressure data of sodium, 
were forced to assume an average heat capacity of 
Although their calculations are based 
function for the 


the liquid. 
on an approximate partition 
vapor which is not accurate at these temperatures, 
the effect of this on the vapor pressures is probably 
relatively small, and they were able to show that 
to 500° C 


for the range 100 an average heat ca- 


pacity of liquid sodium of 1.326 j g~', deg™', based on 
the measurements of Dixonand Rodebush [15], have 
with the data 
theoretically caleulated ‘chemical constant” 


the latter being litaka’s 


better agreement vapor-pressure 
na value of 1.452, 
i} mean value. The results of the present in- 
stigation lead to an average heat capacity of 


10 in this temperature interval. 


X. Summary 


he enthalpy of sodium has been measured in 
to 900° C and with an accuracy of 
to 0.2 percent over most of the range. Values 
heat capacity are derived from the enthalpy 


range 0 


ies With an accuracy of 0.3 to 0.4 percent over 
The triple point of sodium 
+ (2).03° C 


of fusion is 113.2+0.4 abs j/g. 


of the range 


shown to be 97.80 whereas the 


Specific Heat of Sodium 


This investigation was made possible by the 
splendid cooperation of the Knolls Atomic Power 
Laboratory, Schenectady, N. Y., in preparing 
the samples and developing techniques of handling 
and sealing the samples without changing the 
mass of the containers. For this assistance, the 
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Mannal and William D. Davis, who planned and 
supervised the work; and to William A. Ruggles, 
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High-Temperature X-Ray Study of the System Fe,O-Mn.O, 
By H. F. McMurdie, Barbara M. Sullivan, and Floyd A. Mauer 








A series of compositions in the system Fe;Q.-Mn;0, was prepared and examined at 


high temperatures by X-ray diffraction methods. Although Fe;O, and Mn,;0, form a con- 





tinuous series of substitutional solid solutions stable at room temperatures, those composi- 





tions containing less than 60 percent of Mn;Q, are cubic, of the spinel type, and the rest are 






tetragonal. Those compositions that are tegragonal at room temperature pass through a 





two-phase region when heated and acquire the cubic spinel-type structure at high tempera- 






tures. The temperature range of the two-phase (exsolution) region varies with composition. 






The relation to this system of the minerals vredenburgite and hausmannite is discussed. 


Mason prepared a series of compounds by co- 


I. Introduction 


A few years ago, while the oxides of manganese 





precipitating mixtures of the hydroxides of iron 






and manganese and heating the precipitate at 
approximately 1,200° C to obtain the proper 





were being studied in this laboratory [1],' it was 





found that Mn,O,, on being heated, inverts from 






state of oxidation. In this way, homogeneous 






the tetragonal structure to a cubie form of the . 
ce . oo os » fo 9 
spinel type at 1,170° C. preparations covering the range from 20 to 100 


The low temperature form of Mn,O, (hausman- 






mole percent of Mn,O, were obtained. 





nite) is tetragonal, with a=5.75 A and c=9.44 






























{ (2). A larger unit cell, which is more readily — ee eee 

compared in volume and atomic arrangement to 1100 

that of the spinel-type [3] Mn,O, obtained at 

higher temperatures, may be constructed by con- — out emer 

sidering a equal to the diagonal of the cross sec- _ 

tion of the unit cell normal to ¢. When the two 

forms are compared in this way, it is found that 800 + 

the inversion involves only a slight distortion of 

the structure. As the change from tetragonal _—. : 

to cubie at 1,170° C. requires little change in S cos | 

structure, it is rapid and reversible as would be 5 ' 

expected. ¢ 500 b ; 
Since Fe,O, and Mn,O, are known to combine’ # 

400 } Two ] 
to give cubic spinel-type structures (jacobsite), PHASES 
the inversion of pure Mn,QO, described above 300 | 
suggests that a study of the system Fe,0,- 
\In,O, would be of interest. ' | 

II. Previous Studies 

The main previous work on the system Fe,O,- ° i ws Fy i 1 i } 
\IngO, is that of Mason [4], from whose disserta- FexQ, 10 20 30 40 50 60 70 80 90 Mn,0, 





on the diagram in figure 1 is reproduced. Fy Rs 





Figure 1. Diagram of the system FesOg-MnigO, after 
Vason [4] 





gures in brackets indicate the literature references at the end of this 






The System Fe,O,-Mn.O, - 


These preparations were studied by means of 
powder photographs. The sampies containing 
less than 60 percent of Mn,O, gave cubic patterns, 
whereas those containing more than 60 percent of 
Mn,O, were tetragonal. The tetragonal form 
showed a gradual increase in axial ratio with 
increase in manganese content. The axial ratio 
varied from unity (at 60 percent of Mn,O,) to 1.16 
(for pure Mn,O,). Mason assumed that inasmuch 
as these compounds were formed near 1,200° C, 
they represented the conditions at that temper- 
ature and were not in equilibrium at room tem- 
perature. This 
confirmed by the existence in nature of the mineral 


interpretation seemed to be 
vredenburgite. This mineral has been found by 
X-ray diffraction and microscopy [5,6] to be com- 
posed of two phases, one cubic and the other 
tetragonal. The presence of these two phases 
indicates that the mineral is formed by exsolution. 
Mason, this mineral 
represents equilibrium conditions at room tem- 


According to two-phase 
perature, and on the basis of an analysis of the 
two phases present, he indicated an area of im- 
miscibility in the portion of the system between 
54 and 91 percent of MnO, His approximation 
of the region in which exsolution would occur is 
outlined in figure 1. He did not indicate the posi- 
tion of the boundary between the cubic and the 
tetragonal phases in the one phase region. 

In Mason’s survey of the minerals comprising 
the system Fe,O,-Mn,O,, he fixed limits of com- 
position for four minerals: magnetite, jacobsite, 
vredenburgite, and hausmannite. According to 
him, 
whose composition lies within the area of im- 


“Vredenburgite signifies those substances 


miscibility, i. e., between 54 and 91 percent 
Mn,QO,.”’ 
substances in this system containing from 91 to 
100 percent of Mn,O,. The remainder of the 


system, from pure Fe,O, to 54 percent of Mn,Q,, 


The name hausmannite was given to 


was arbitrarily divided at a composition corre- 
sponding to 50 percent of Fe,O, and 50 percent of 
MnFe,O,, that is, at 16.7 of Mn,Q,. 
Substances containing less than 16.7 percent of 


percent 


Mn,O, were called magnetite, and those containing 
from 16.7 to 54 percent of Mn.O, were designated 
jacobsite. The two-phase mineral vredenburgite 
was believed by Mason to represent equilibrium 
conditions at room temperature and was therefore 
called 8-vredenburgite. The 
burgite was reserved for a high-temperature form 


name a-vreden- 
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that Mason expected to find. Later [7], he dix 
find minerals in this range of composition tha 
consisted of a single tetragonal phase. 

Von Eckermann [8] suggested that the chang: 
from the cubic to the tetragonal structure might 
take place at a different composition if the tempe1 
aturé were increased. Basing his hypothesis on 
Mason’s paper, and on earlier work by Montoro 
|9], he suggested that compositions high in Mn,O 
might remain tetragonal up to the melting point 
whereas those richer in iron would change on 
heating from cubic to tetragonal at approximately 
1,000° C. As Mason pointed out in a later paper 
[10], such a change to lower symmetry on heating 
seems unlikely. 

When work in this laboratory disclosed that 
heating pure Mn,O, brought about a change from 
the tetragonal to the cubie structure at 1,170° ¢ 
it became evident that further work was needed 
to clarify the picture of the high-temperature 
Because of the rapid 


portion of the system 
inversions that occur in this system, it is not 
practical to investigate the phase relations by 
examination of quenched samples. High-tem- 
perature X-ray diffraction studies seemed to be 
the only practical method of attacking the prob- 
lem, and, since suitable equipment was available 
[11], work 
of the system. 


was undertaken to revise the diagram 


III. Method and Materials 


The high-temperature X-ray diffraction appara- 
tus described by Van Valkenburg and McMurdie 
[11] was employed to hold the specimens at various 
temperatures while the diffraction patterns were 
being made. Unfiltered FeK radiation was used, 
and the patterns were automatically recorded on 
a Geiger-counter spectrometer. The equipment 
used can record only those lines for which the 
Bragg angle is less than 45°. For this reason, it 
was not possible to determine the cell dimensions 
with the accuracy attained by Mason. 

A series of mixtures at intervais of 5 percent 
by weight ? from 60 to 100 percent of Mn,O, was 
prepared by a method similar to that of Mason 
Manganous and ferrous sulfates were weighted 
out in the proper proportions and dissolved 


The hydroxides were coprecipitated with NaOH 


2 Although Mason's compositions are expressed in mole percent and th 
in this paper in weight percent, the difference in this system is very slig! 
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nd washed by centrifuging until free from sul- 
es. The 
200° C for approximately 18 hours to convert 
hem to the desired state of oxidation. 


precipitates were then heated at 
Analysis 
several of the samples for iron and manganese 
howed that the composition was as originally 
janned in each case 
X-ray diffraction patterns of each of the samples 
vere recorded at room temperature and at various 
levated When a 


ecorded at a high temperature, the sample was 


temperatures. pattern was 
isually held at that temperature for about half 
an hour. Recording the pattern required about 
1) minutes, after which the sample was raised to 
he next temperature without intermediate cooling. 
In a few cases samples were held at one tempera- 
ture for periods up to 24 hours. The change in 
the relative amounts of the two phases during 


this period was so slight that it seemed reasonable 


TABLE | 


The System Fe.O,-Mn.O, 


to assume that the samples were at equilibrium 
when held at one temperature for half an hour. 

In order to prevent oxidation of the sample at 
temperatures between 200° and 1,000° C, it was 
necessary to pass a stream of nitrogen through 
the furnace. In no case where nitrogen was used 
did lines of the higher oxide appear. 


IV. Results 


Table 1 shows the phases present and the unit- 
cell dimensions * for each of the samples at various 
that, at room 


temperatures. These data show 


temperature, each sample consists of a single 


phase, the one containing 60 percent of Mn,O, 
When 


the samples containing 65 to 95 percent of Mn,O, 


being cubic while all others are tetragonal. 


As mentioned previously modified unit cell s been used to describe 
tetragonal phase l 


tures 


Unit-cell dimensions of com positions in the sustem Fe ( ».-Mn ), at various te m pe ratures 





are heated, each passes from the tetragonal phase 
through a region in which both the tetragonal 
and cubic phases are in equilibrium. At high 
temperatures the cubic phase alone is present. 
The pure Mn,O,, of course, passes directly from 
the tetragonal to the cubic form with no two- 
phase region. As any one sample is heated to a 
temperature near the inversion point, the axial 
ratio does not approach unity as one might expect. 
The only change in unit-cell dimensions is a 
slight over-all increase due to thermal expansion. 
For this reason, the cubic and tetragonal patterns 
do not become indistinguishable at the inversion 
point. 

Although the 60-percent Mn,O, sample gave a 
cubic pattern at room temperature, this pattern 
was characterized by poorly defined peaks and 
considerable line broadening. The fact that the 
pattern improved when the temperature of the 
sample was raised to 600° C indicates that the 
cubic phase is probably not stable at room tem- 
perature. The cubic pattern obtained can be 
attributed to the fact that the samples were held 
at 1,200° C during preparation, at which tem- 
perature all would exist in the cubic phase. The 
60-percent-Mn,O, sample appears to have been 
quenched when cooled to room temperature, and 


it was impossible to observe any transition on re- 
heating. Extrapolation from data obtained for 
samples of higher Mn,O, content indicates that, 
in the case of the 60-percent-Mn,O, sample, the 
cubie and tetragonal phases should be in equilibri- 
um at about 400° C. The specimen was held at 
this temperature for a period of 24 hours, but 
there was no significant change in the pattern 
obtained at the end of this time. The furnace 
temperature was then manually lowered in small 
steps over a period of 2 days. When the temper- 
ature of the specimen was lowered to 250° C the 
lines characteristic of the tetragonal phase were 
discernible. The specimen was eventually cooled 
to room temperature, but the tetragonal pattern 
was still not distinct, and the cubic pattern showed 
some line broadening, indicating that the lattice 
was under strain. Thus, while it seems reason- 
able to assume that the cubic form containing 60 
percent of Mn,O, is not at equilibrium at room 
temperature, the tendency to supercool makes it 


impossible to obtain a specimen that is tetragonal. 


Data on samples containing 50 and 55 percent of 
Mn,O, would be desirable, but because of the 
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slower rate of inversion of compositions in thi 
range, it does not seem likely that one could b 
reasonably sure that the patterns obtained repre- 
sented the conditions at the temperature use: 
during the test. 

In figure 2 the results obtained from the pat- 
terns made at elevated temperatures are plotted 
Points have been plotted for each of the samples 
to show the highest temperature at which th 
tetragonal phase alone was discernible, the 
lowest and highest temperature at which both the 
cubic and tetragonal phases could be detected, and 
the lowest temperature at which the cubic phase 
appeared alone. It should be borne in mind that 
a second phase may not be detected if it constitutes 
less than 10 percent of the specimen. Figure 2 
shows the temperature range over which two 
phases could be distinguished, but it must be 
assumed that traces of the second phase too small! 
to detect may be present outside this temperature 
range. As explained above, no satisfactory data 
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Figure 2. Partial diagram of the system FexsOq-Mng0. 


, Highest temperature at which the tetragonal-phase alone was discernibl 
A, limits of the temperature range over which both the cubic and tetre 
gonal phases could be detected; ©, lowest temperature at which the cub 
phase appeared alone. 
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CELL DIMENSIONS. A 


DIMENSIONS, A 


CELL 





ere obtained for the 60-percent sample. The 
otted portion of the curve indicates the probable 
mits of the two-phase region at this end of the 
irve. As shown in the figure, the pure Mn,O, 
nverts at 1,170° C and has no two-phase region. 
The cell dimensions at room temperature are 
hown in figure 3, and the dimensions for both 
shases occurring at 600° C appear in figure 4. 
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Unit-cell dimensions of the system FesOy- MgO, 
at 600° C. 


The System Fe,.O,-Mn.O, 
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V. Discussion 


Whereas the end compounds of the series 
Fe,O,-Mn,O, are of the cubic and tetragonal 
structure, respectively, at room temperature, the 
unit cells of the two compounds have been shown 
to be very similar. As iron and manganese have 
similar atomic radii, it is not surprising to find 
that the two oxides are completely miscible at 
room temperature and form a series of substitu- 
tional solid solutions. If the unit-cell dimensions 
at room temperature are plotted for the series of 
compositions between Fe,O, and Mn,O, [4], it is 
found that the size of the cubic cell increases 
regularly until a composition of about 60 percent 
of Mn,O, is reached, at which point the cubic 
lattice is no longer stable. When the composition 
exceeds 60 percent of Mn,Q,, the cubic lattice 
becomes deformed and a tetragonal lattice results. 
The axial ratio increases continuously from unity 
at about 60 percent of Mn,QO, to 1.15 at 100 per- 
cent of Mn,O,. The volume of the unit cell in- 
creases almost linearly with the substitution of 
manganese throughout the whole range of compo- 
sitions. This behavior is typical of substitu- 
tional solid solutions. 

Although Fe,O, and Mn,O, appear to be com- 
pletely miscible in the solid state at room tem- 
perature, a two-phase region has been found at 
higher temperatures. For example, at 600° C 
manganese may be substituted for iron to obtain 
a single cubic phase until a composition of about 
69 percent of Mn,O, (by weight) is reached. 
Likewise, iron may replace manganese in the 
tetragonal structure of Mn,QO,, giving a single 
tetragonal phase. When a composition of 81 
percent of Mn,O, and 19 percent of Fe,Q, is 
reached, however, no more iron will go into the 
structure of Mn,QO,. For this reason, there is a 
range of compositions at this particular tempera- 
ture within which each sample is a mixture of two 
phases. 

The vredenburgite occurring in nature usually 
shows two phases, cubic and tetragonal. Mason 
believed that this two-phase mineral, rather than 
the prepared samples, represented equilibrium 
conditions at room temperature, and that the 
compositions of the two phases making up the 
mineral therefore represented the limiting values 
of the solubility of Mn,O, in Fe,O, and of Fe,O, in 
Mn,Q, at normal temperatures and pressures. — In 
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his survey of the minerals comprising the system 
Fe,Q,-Mn,O,, he 


those substances whose compositions lie within 


designated as  vredenburgite 
the area of immiscibility, that is, between 54 and 
9] percent of Mn (), 


by examining samples of vredenburgite and com- 


He arrived at these limits 


paring the lattice dimensions of each of the two 
phases present with those of a series of prepara- 
tions that covered all compositions from Fe,O, to 
Mn,O, at intervals of 10 mole percent. 

The present investigation shows that there is no 
region of immiscibility at room temperature and 
that Mason’s criterion for establishing the limits 
of composition for vredenburgite is not based on 
equilibrium conditions at room temperature 

Mason also found specimens in the range from 
54 to 91 percent of Mn,O,, which consisted of a 
single tetragonal phase. Believing that the single- 
phase type was a high-temperature form not in 
equilibrium at room temperature, he proposed 
ealling it ea-vredenburgite and reserved the name 
8-vredenburgite for the two-phase material. In 
choosing these designations, Mason followed the 
convention that the high-temperature form should 
be indicate!’ by the prefix a. However, the 
present study indicates that the two-phase form 
(8-vredenburgite) is the one that is stable at 
higher temperatures. To prevent confusion, it 1s 
suggested that Mason’s designations be retained, 
the two-phase material being called 8-vreden- 
burgite 

Figure 5 is a photomicrograph of §8-vreden- 
burgite from India. The cubic phase appears as a 
black background, upon which the bright lines 


As this 


material is unstable at room temperature, we must 


due to the tetragonal phase may be seen 


assume that, in cooling, the specimen remained in 
the two-phase region long enough for fairly large 
phase to form by 


ervstals of the tetragonal 


exsolution. Subsequent cooling must have been 
rapid enough so that reaction between the two 
phases could not proceed and a degree of quench- 
ing took place 

Figure 6 is a photomicrograph of natural 
hausmannite from Sweden that shows polysyn- 
thetic twinning, as has been noted previously, 15] 
typical of material that has undergone a simple 
inversion [12]. In this case, the material must 
have had the cubic spinel structure before cooling 
through the inversion temperature 


Because of the difficulty in distinguishing iron 
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Figure 5 Photomicrograph of B-vrede nburgite from Indi 


by reflected polarized light 


rhe tetragonal phase appears as bright lines an 


dark background x 1S0 


from manganese by their difference in scattering 
power, no information could be obtained about thi 
distribution of cations among the positions avail 
able in the spinel structure. Barth and Posnjak 
[13] have shown that even the valence of a par 


i 


ticular atom does not determine which metal 


position it will assume 


VI. Summary 


At room temperature, FesQ. and Mn,Q, appear 


to be completely miscibi a the solid state 


Manganese may be substi ited for iron in the 


Ficure 6 Photomicrogr tph by reflected polarized lig 


showing polysynthetic twinning of hausmannite fro 


Swede n 180 
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ibie spinel-twpe structure of Fe,O, until a com- 
osition of about 60 percent by weight of Mn,Q, is 
ached. Compositions containing more manga- 
ese take on the tetragonal structure of pure Mn,Q,, 
nd the axial ratio increases continuously from 
nity at 60 percent of Mn,O, to 1.15 at 100 percent 
The volume of the unit cell increases regularly 
th the addition of hanganese. 
In the temperature range from 380° to 1,170 
compositions above 60 pereent of Mn,0, are 
iot completely miscible. When a composition 
with the tetragonal structure is heated, a tem- 
perature is reached at which a cubic phase of the 
spinel type starts to form. This cubic phase is in 
equilibrium with the tetragonal phase over a 
range of temperatures that depends on the compo- 
sition The percentage of tetragonal material 
decreases with further heating until a temperature 
is reached at which all the material has the cubic 


spinel-tvpe structure. From this temperature up 


to the melting point, only the cubic phase is 
present. 

The natural two-phase mineral 6-vredenburgite 
s not at equilibrium at room temperature. The 
one-phase mineral a-vredenburgite is the low- 


temperature form 
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A Quantitative Study of the Carbon Monoxide 
Formed During the Absorption of Oxygen 
by Alkaline Pyrogallol 


By Marthada V. Kilday 


The amounts of carbon monoxide evolved during analytical absorptions of oxygen by 
alkaline pyrogallol were studied with respect to the composition and methods of preparing 
solutions, the temperature of reaction, the rate of flow of the sample through the solution, 
the capacity of the reagent for the absorption of oxygen, the partial pressure of reacting 
oxygen, and the manner of bringing the gas into contact with the reagent Less than 0.02 
percent of carbon monoxide was formed during analyses of gases containing as much as 80 
percent of oxygen when the recommended solutions were used in bubbler-type absorption 
pipettes, and when the samples were passed into the solutions at a rate greater than 20 
milliliters per minute at about 25° C. This was true until the solutions had absorbed 
approximately twelve times their volume of oxygen, after which the formation of carbon 
monoxide increased significantly. The NBS colorimetric carbon monoxide indicating gel 
was used for the quantitative determinations of carbon monoxide. A method was developed 
for preparing pyrogallol solutions that vielded reproducible amounts of carbon monoxide 


so small that no significant error was caused in the usual analysis 


unfounded, as no significant amount of carbon 


I. “atroduction aed, a 
monoxide is formed under the conditions of the 


For many vears it has been known that under usual analvsis when the recommended solutions 


some conditions the analytical absorption of — are used. 


oxygen by alkaline pyrogallol is subject to an error 





Many investigators have tried to determine the 
introduced by the formation of carbon monoxide — amount of carbon monoxide formed under the 


This carbon monoxide is 


during the reaction. 
measured as part of the residue after absorption, 
so that the oxygen absorbed appears to be less 
than was actually present in the original sample. 
lf a combustion analysis follows, the carbon 
monoxide may contribute further error. 
Pvrogallol is so superior to other reagents some 
times employed for the removal of oxygen, with 
spect to the exchange of gases inert to it, that 
rrors of physical 


solubility are very greatly 


duced by its use [1].'. However, the evolution 
{ carbon monoxide has long been invoked as an 
irgument against the use of pyrogallol. It will 


shown in this paper that 


this argument is 


Figures in brackets indicate the literature references at the end of this 


Carbon Monoxide From Alkaline Pyrogallol 


conditions of the conventional analytical proced- 
ures, but usually they have found that the con- 
centration of carbon monoxide was too low for the 
methods of analysis employed. Therefore they 
have resorted to measuring the relatively large 
amounts of carbon monoxide produced under 
extreme conditions, such as high partial pressures 
of oxygen, weakly alkaline pyrogallol solutions, 
and high temperatures of reaction. Even with 
these high concentrations of carbon monoxide, 
which came within the range of their analytical 
methods, there was a great divergence in the 
results obtained by various investigators. 

During World War II, this Bureau developed a 
colorimetric indicator for the quantitative deter- 


mination of carbon monoxide in low concentrations 
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up to 0.1 percent. Higher concentrations were 


determined by making simple dilutions of the 
samples |2]. 

The NBS 
satisfactory in the direct analvis of the carbon 
the residues of oxygen 
During the prelim- 


indicator was found to be very 


monoxide present in 
absorptions by pyrogallol. 
inary work, it was found that, although a given 
solution \ ielded reproducible amounts of carbon 
monoxide during a series of like analyses, there 
was significant variation in the amount of carbon 
different solutions pre- 
This led to 
an investigation of the methods of preparing the 
dissipating the heat 


monoxide generated by 
pared according to the same formula. 


pyrogallol solutions. By 
produced while dissolving the pyrogallol in a 
saturated solution of potassium hydroxide, solu- 
tions were prepared that not only gave repro- 
ducible results with respect to carbon monoxide 
generated, but which consistently evolved less 
carbon monoxide during an absorption. 

A satisfactory method of preparing the pyro- 
gallol solutions having been established, it was 
possible to proceed with the determination of 
carbon monoxide evolved under various condi- 
tions. The data obtained in this experimental 
study are reported so that the gas analyst may use 
them in estimating the error to be expected as a 
result of carbon monoxide formation under the 
conditions of the analysis. 


II. Description of Apparatus 


A modification of the Shepherd gas analysis 


apparatus [3] with spherical, interchangeable 
ground-glass joints instead of rubber connections 
was used in this study, and Apiezon grease “L’”’ 
was the stopcock lubricant employed. Figure | 
is a diagram of the distributor as modified for this 
work. Below stopcock | is the absorption pipette. 
Two mercury reservoirs, connected under stop- 
cocks 2 and 3, were used for diluting and mixing 
the residues from the absorptions of oxygen, which 
were too rich in carbon monoxide to be analyzed 
directly with the NBS colorimetric 
Below stopcock 4 is a tube containing gold sponge 


indicator. 


to remove mercury vapor, which interferes with 
the test for carbon monoxide. The colorimetric 
carbon monoxide indicator tube was attached with 
a pure gum rubber connection below the gold 
sponge Between stopeocks 5 and 6 was a eali- 
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— COLORIME TRIC 
INDICATOR 


Figure 1. Distributor showing modifications provided fo 


the carbon monoxide determination and for measuring the 


rate of flow of the gas. 


brated flowmeter that was provided to measure the 
rate of flow of the residue after absorption through 
the carbon monoxide indicator tube. The position 
of the stopecocks could be manipulated so that the 
flowmeter containing mercury was either included 
from the distributor train. The 
previously de 


or excluded 

burette and 

scribed [4]. 
The volumetric analyses were made in a lab- 


compensator were 


oratory where the temperature was maintained at 
about 25° C. 
for making physical measurements and was there- 
fore free of the atmosphereic contaminations and 


ry’ -- > 
The room was reserved primarily 


variations present in many chemical laboratories 
The nine types of absorption pipettes used in 
this investigation are shown in figure 2. 


III. Method for Determination of Carbon 
Monoxide 


The NBS colorimetric indicator for carbon 
monoxide and the “Laboratory Method” deseribed 
by Shepherd [2] were used in the analyses for car- 
A detailed 


description of the analytical procedure as adapted 


bon monoxide reported in this paper. 


to this problem is given in section V. 

The chemicals of the indicator tube are inclosed 
in a 15-em length of 7-mm Pyrex glass tubing 
The 2-em length of vellow indicating gel, which 
turns green upon exposure to carbon monoxid 
is protected at both ends by a 5-cm length o! 
purified silica gel, which removes water vapor and 
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organic vapors that might interfere with the deter- 
mination of carbon monoxide. The guard gel is 
held in place by absorbent cotton, and a cork is 
placed in each end when the tube is not actually 
being used. 

To determine the percentage of carbon monoxide 
present in an unknown sample, a measured volume 
of the sample is passed through the indicator tube 
at a controlled rate of flow. The color developed 
by the indicating gel is compared with standards 
prepared by exposing tubes to known concentra- 
tions of carbon monoxide. The percentage of 
carbon monoxide expressed in this paper is based 
on the volume of the original sample and refers to 
the amount of carbon monoxide that would have 
been generated by 100 ml of original sample. In 
other words, if it »s stated that a solution generated 
0.02 percent of carbon monoxide, it means that 
0.02 ml of carbon monoxide would have been 
present in the residue of absorption of 100 ml of 
original sample. 

In this work there was a linear relationship 
between the concentration of carbon monoxide 
and the probable error of the carbon monoxide 
determination. From 0.0020 to 0.039 percent of 
carbon monoxide, the error was about 17 percent 
of the carbon monoxide present. These errors 
could have been reduced to less than half of these 
values (1) by reducing the interval of time in 
exposing standards to the known mixtures of 
carbon monoxide and (2) by having the concen- 
tration of carbon monoxide in the known mixtures 
nearly the same as the partial pressure of carbon 
monoxide that exists in the unknown. These 
refinements were not employed in this work be- 
cause most volumetric gas analysis apparatus are 
If the analyst is 
concerned with extreme accuracy, a determination 


not capable of such accuracy. 


of carbon monoxide should be made for each indi- 
vidual analysis, as there is some variation in the 
amount of carbon monoxide produced even when 
the conditions of the analyses are similar. 


IV. General Methods for Preparing 
Solutions 
Prior to the present study, the gas chemistry 
laboratory had used a method of preparing pyro- 
gallol solutions that was previously described [3] 
and is quoted below (page 145 of this reference): 


The reagent used for the removal of oxygen in the 


absorption pipette (platinum-disk bubbler) is 
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an alkaline pyrogallol solution. This is made by dis- 
solving erystalline pyrogallol in a minimum amount 
of hot water and adding this solution to the concen- 
trated potassium hydroxide solution ... The re- 
sulting solution closely approximates one recom- 
mended by Anderson," who found a solution con- 
taining 15 g of pyrogallol per 100 ml of potassium 
hydroxide solution of 1.55 specific gravity to be 
satisfactory. 

In order to prevent air from reaching the reagents 
a laver of light mineral oil, about 10 to 15 mm in depth, 
is poured over the reagent in the reservoir arm of the 
pipette. This also permits the preparation of the 
pyrogallol solution within the pipette itself and out 
of contact with air. 
tion is first introduced, the oil is placed over this 


The potassium hydroxide solu- 


solution and the pyrogallic acid solution in water is 


introduced beneath the layer of oil from a pipette. 


® Anderson, R. P., Ind. Eng. Chem. 7, p. 588 (1915 


It was necessary to heat the pyrogallol crystals 
with water in order to prepare the solution within 
a reasonable length of time, and this solution, 
while still hot, was poured through a long-stemmed 
funnel to the bottom of the saturated potassium 
hydroxide solution in the pipette. This was fre- 
quently accompanied by considerable bubbling 
and splattering as a result of the additional heat 
produced during solution. 

Preliminary work of this study indicated that 
this method of preparation produced solutions 
vielding inconsistent results with respect to the 
Con- 
sequently, an improved method was developed 
that produced solutions that vielded reproducible 


concentration of carbon monoxide evolved. 


amounts of carbon monoxide when used under 
similar conditions. This solution, quite happily, 
vielded significantly less carbon monoxide than 
previous solutions. A description of the im- 
proved method follows. 

The saturated potassium hydroxide solution is 
prepared by saturating distilled water with solid 
potassium hydroxide at 25° +1° C. This solu 
tion contains about 55 percent of potassium hy- 
droxide, or 120 g of potassium hydroxide per 100 
g of water, and has a specific gravity of about 1.55 

The flask or pipette in which the pyrogallo! 
is to be stored is placed in an ice wath, and th 
pyrogallol crystals are placed in the bottom of th: 
container. The saturated solution of potassium 
hydroxide is added rapidly, and the flask is sealed 
with an air-tight cover, or a layer of light minera 
oil is added if the solution is prepared in a pipett: 
The new solution should remain in the bath unt 
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he ice melts and the bath reaches room tempera- 
ire, at which time it may be used immediately or 
stored for future use. If any crystals remain 
indissolved, agitating the solution in the pipette 
V passing nitrogen into it will complete solution. 
Data will be presented that show no significant 
lifference in the amount of carbon monoxide 
ormed by solutions that were freshly prepared 
ind those that were stored out of direct sunlight 


for as long as 21 months. 


V. Procedure for Analyses 


In the course of this study, five partial pressures 
of oxygen were used: approximately 21 (air), 40, 
60, 80 percent of oxygen in nitrogen, and 100 per- 
cent of oxygen. Air was sampled directly from 
the laboratory in which the analyses were made, 
and the other mixtures and the approximately pure 
oxygen were stored under pressure in gas cylinders 
from which samples were taken in the following 
manner. A needle valve was used to control the 
sampling rate, and the sample was passed through 
and connected by 
The sample 


copper tubing to A, figure 1, 
means of nitrometer pressure tubing. 
was thoroughly flushed from A through B, and 
stopcock ( was turned to admit the sample under 
pressure to the mercury-filled burette. 

The sample was saturated with water vapor from 
the burette walls, and the volume was measured in 
the usual way. Unless otherwise stated, the rate 
of flow of the gas to the absorption pipette was 
controlled only by the maximum rate at which the 
mercury could flow into the burette and by the 
back pressure exerted in the pipette; this was 
usually about 190 to 250 ml/min. After enough 
passes had been made into the pyrogallol to assure 
complete absorption of the oxygen, the residue was 
returned to the burette, saturated, and measured 
Ls before 

From this point the procedure depended upon 
the total volume of the residue as well as the 
concentration of carbon monoxide. Dilution was 
if the amount of carbon monoxide 
xceeded 0.04 percent, or if..the residual volume 


hecessary 


vas less than 70 ml, which was the minimum set for 
If no dilution 
vas needed, the determination of carbon monoxide 


ih determination in this work. 


vas made immediately. 
When the volume of the residue was too small, it 
vas measured at atmospheric pressure and air was 
awn into the burette. The volume of the result- 
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ing mixture was also measured at atmospheric 
pressure, and the concentration of carbon monoxide 
present in the residue was computed. 

If the concentration of carbon monoxide in the 
residue after absorption was expected to exceed the 
range of the colorimetric carbon monoxide indi- 
cator, a fraction of the residue was stored in one 
of the mercury reservoirs, and the remaining 
portion was diluted and treated as in the procedure 
described above. If further dilution was neces- 
sary, the second mercury reservoir was used to 
store the residue from the first dilution, and a 
measured volume of air was taken into the burette 
and mixed with the first diluted residue. 

The amount of carbon monoxide in the residue 
or in the diluted 
follows. Six indicator tubes were used in the 


residue was determined as 


analysis, one for the residue and five for color 
standards. An indicator tube was attached to the 
volumetric gas analysis apparatus as shown in 
figure 1. The distributor was flushed, and 45 ml 
of the residue was passed through the indicator 
tube at 90 ml/min, which is equal to a 30-sec 
exposure. The standards were prepared by pass- 
ing a known concentration of carbon monoxide 
through indicator tubes at 90 ml/min for varying 


periods of time—usually at 5-sec intervals. The 
five standard mixtures of carbon monoxide in air 
were: 0.0025-, 0.0052, 0.0097, 0.020,, and 0.038, 
percent. The percentage of carbon monoxide in 
the unknown (residue) was determined by com- 
paring visually the color developed in this tube 
with the colors of the tubes that had been exposed 


to known concentrations of carbon monoxide. 


VI. Data and Discussion 
1. Effect of Variations in Pyrogallol Solutions 
(a) Experimental Methods of Preparation 


An inspection of the grap)is in figure 3 reveals the 
inconsistencies found in the amount of carbon 
monoxide formed by solutions prepared by the 
“old” method previously described. These graphs 
show the amount of carbon monoxide evolved 
during absorptions of oxygen from air by solutions 
of various compositions individually prepared in 
several platinum-disk bubbler pipettes (see table 
1). It may be seen that the amounts of carbon 
monoxide evolved during absorptions of oxygen 
from air under similar conditions by any one 
solution (those in any one small block) were 
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£20 4.024 028 6.032 0036 0 004 008 O12 O16 020 024 .028 .032 
co PRODUCED,% 


Percentage of carbon monoxide produced during absorption of oxygen from air by solutions prepared according t 


the old method. 


m of the solution The number in the upper left-hand corner of each small block indicates which platinum-disk bubbk 
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same composition 
ock) the amount of carbon monoxide produced 
metimes varied by as much as a factor of 10 


TARLE 1. Composition of solutions studied 


Solution 


Pyro 
gallol per 
“ le { 
Compe I Potassium ml o 
n Ne " hydroxide saturated | 
: . H,0 KOH so- 
iturated . 
mo lution * « 
at 2 
nl ni 
+ 235 0 
s-2 ay 2 
s 17 “w) 
4 1™) 75 
. 12 100 
=) ys 
s.7 v1) 
S-X 4 17 
s-Z4 4 
s-10 
S-B7 
S-BZ 
e crystalline pyrogallol was obtained commercially under the name 
Acid pyrogallic, analytical reagent." When the old method of preparing 
ns was used, an additional 10 ml of water was added to dissolve the 


line pryogallol 








he amount 0 rogallol was equal to about 15.5 g of pyrogallol per 100 ml 
ed solution of saturated potassium hydroxide Thus the solutions of 
ympositions would logically be placed between solutions S-X and S-Z, 
for the dilutions of the saturated potassium hydroxide solution. 

oncentration of potassium hydroxide was the same as that used in 
1 8-1, which was four woduce the least carbon monoxide in the 
howing the effect of changing the concentration of potassium hydroxide 
olution 

he only significance attached to S-X and S-Z is that the former was of a 
ition in general use by this laboratory before this work, and the latter 

uind to be the most desirable composition when the improved method of 

ring solutions Was u 

solution S-7 stored in an Erlenmeyer flask and the cap sealed with 
ffir 
s of solution S-Z stored in an Erlenmeyer flask and the cap sealed 
A zon grea I 


The fact that the variations were not caused 
by the use of different platinum-disk bubbler 
pipettes is shown in solutions S—3, S-4, S—5, and 
S-7 where inconsistencies were found when the 
sume pipette was used for different preparations 

solutions of the same composition, 

\ thorough investigation was made of the 

ssible er cses for the variations shown in figure 

First, the accuracy of the determination of 
bon monoxide was verified by taking into the 
irette samples of a known concentration of 
rbon monoxide and determining the carbon 
moxide present by the same method used in 


analysis of the residues of absorption. In the 
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ssonably consistent, but in various solutions of 
those in any one large 


ten analyses made, the maximum deviation from 
the known value was 10 percent of the carbon 
monoxide present, with an average deviation of 
6 percent. The accuracy of the determination of 
carbon monoxide was therefore discounted as a 
Other 


which will be discussed in detail later, were: 


significant variable. possible variables, 


1. Concentration of the solution with re- 
spect to pyrogallol. 

2. Concentration of the solution with 
respect to potassium hydroxide. 

3. Rate of flow of sample into pyro- 

gallol solution. 

4. Age of solution (actual age and also 
the degree expended ; 

5. Temperature of reaction during ab- 
sorption of oxygen. 

6. Number of passes into pyrogallol solu- 
tion. 


These factors were controlled sufficiently to 
eliminate them as significant variables in the runs 
represented in figure 3. This left the actual 
method of preparing the pyrogallol solutions as 
the most probable cause of the extreme deviations. 

It was believed that variations in the temperature 
of solutions might have caused differences in the 
solutions themselves, which in turn caused 
differences in the amount of carbon monoxide 
evolved during the absorption of oxygen in differ- 
ent runs. As the measurement of the maximum 
temperature developed at several points would 
have been rather difficult in this instance, the 
effect was shown by preparing solutions in which 
it was known that the temperature of the solu- 
tion had been varied by large amounts. 

Figure 4 shows the amount of carbon monoxide 
formed in absorptions from gaseous mixtures 
containing 21 percent and 80 percent of oxygen 
by solutions that were prepared in five different 
ways, each varying with respect to the tempera- 
ture allowed to develop during preparation. 

The greatest amount of carbon monoxide per 
unit volume of oxygen was evolved by the solution 
prepared by method 2 (see fig. 4) in which the 
pyrogallol was added all at once to the saturated 
potassium hydroxide solution at 53° C in an 
atmosphere of nitrogen. The yield of carbon 


monoxide was 0.023 to 0.059 percent during 


absorptions of oxygen from air 
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Percentage of carbon monoaride produced during absorption of} o1~ygen by solutions prepared according to vario 


methods. 


The analysis number on the | shows the order in which the analyses were performed 
i}, method 1; 2, 8-X, 8G, method 2; 3, 8-X, 8G, method 3; 4, 8-X, PD-4, method 4 


ire indicated as follows: 1, 5-4, 
air); @, | Oy: in N 


In method 1 the pyrogallol was dissolved, by 
heating it with a little water (to approximately 
ie” © 
saturated potassium hydroxide solution at room 


and added in small increments to the 
temperature, with vigorous stirring after each 
addition. In this case, the dissipation of heat 
was more rapid, and the carbon monoxide evolved 
during the absorption from a mixture containing 
80 percent of oxygen was about 0.05 to 0.06 per- 
cent, which is of the same order as that produced 
by the previously described solution when only 21 
percent of oxygen was absorbed. Since the higher 
partial pressure of oxygen will produce more carbon 
monoxide in the same solution, it was assumed 
that the carbon monoxide evolved by the solutions 
during absorptions was probably dependent to 
some extent upon the rapidity of the dissipation 
of heat during the preparation of the pyrogallol 
solutions. 

In preparing the solution according to method 3, 
it was found unnecessary to dissolve the pyrogallol 
crystals in water, as previously described, before 
adding them to the saturated potassium hydroxide 
solution. The solution was made by simply add- 
ing the crystalline pyrogallol to the saturated solu- 
tion of potassium hydroxide at room temperature 
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The composition number of the solutions and the pipettes uss 


5, 8-X, PD-7, method 5 . 21% O, 


and placing the mixture in a refrigerator overnight 
to accelerate the dissipation of the heat of reac- 
tion. Here the carbon monoxide evolved was re- 
duced by more than half of that found in method | 
in the absorption from a mixture containing 80 
percent of oxygen, and the reduction in carbon 
monoxide formed in the absorption of oxygen from 
air was remarkable, as little as 0.002 to 0.003 pet 
cent being generated. 

In method 4 the pyrogallol crystals were added, 
while stirring, to the saturated potassium hy 
The solu 
tion remained in the bath until it reached room 
temperature. With 21 
solution vielded about the same amount of carbo 


droxide solution in an ice bath at 8° C. 
percent of oxygen, this 
monoxide as in method 3; with 80 percent ol 
oxygen, the yield was appreciably less. 

Method 5 was found to be the most convenient 
The absorption pipette containing the pyrogallo! 
crystals was nearly immersed in an ice bath, and 
the saturated solution of potassium hydroxide, at 
room temperature, was added. This solution re 
mained in the bath until it reached room tempera 
ture. The reproducibility in the amount of ca 


bon monoxide formed by solutions prepare: 


according to this method was very satisfactory 
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d the yield of carbon monoxide was low, about 
02 to 0.004 percent. None of the solutions 
owed variations in the amount of carbon mon- 
de formed that were greater than those that 
curred in different analysis made with a single 
ution. The reproducibility will be discussed 
more detail later. 


(b) Composition 


The composition for the solutions studied in 
his investigation are given in table 1, and refer- 
ence to the solutions of these compositions will be 
ade only by number. The table is arranged to 
show the manner of preparing the solutions as 
vell as the composition. 

\lthough the pyrogallol solutions represented in 
figure 3 were prepared according to the old method, 
they may be used to illustrate the general effect of 
changing the concentrations of potassium hydrox- 
ide and pyrogallol in solutions, if the obviously 
erratic runs are eliminated. In solutions S-1 
through S-5, the pyrogallol remained constant 
while the concentration of potassium hydroxide 
Solution S-1 yielded 
().002 to 0.006 percent of carbon monoxide, and as 


became progressiv ely less. 


the alkalinity decreased, the carbon monoxide 
produced became greater until in S—5 about 0.01 
to 0.02 percent of carbon monoxide was formed 
during air analyses 

The concentration of potassium hydroxide 
remained constant in solutions 6, 7, X, and 10, 
but the weight of pyrogallol was varied (see table 
|). The changes in the concentration of pyro- 
gallol (in the solutions used) apparently did not 
affect the amount of carbon monoxide produced 
during absorptions of the oxygen in air. This 
fact was also supported by the analyses in figure 5, 
which represents solutions prepared according to 
method. The 
sented by blocks 1 and 2 contained 17 g of pyro- 


the improved soultions repre- 
gallol per 100 ml of saturated potassium hydroxide 
solution and were prepared in the same pipettes 
is those in blocks 5 and 10, the solutions of which 
contained only 14 g of pyrogallol per 100 ml of 
iturated potassium hydroxide solution. 

Increased viscosity and clotting occurred sooner 
having higher concentrations of 
The clotting impaired the efficiency 


solutions 
vrogallol. 
' the distributor pipettes by reducing the number 
openings through which bubbles could escape, 
trapping the fine bubbles in the viscous solu- 
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tion, and by reducing the rate of flow of the sample 
into the solution. Finally, the distributor open- 
ings became plugged and no gas at all could be 
bubbled into the solution 


in figure 5, except block 1, oxygen was absorbed 


In all solutions shown 


by the solution until the distributor became so 
badly plugged that the sample could not be 
passed into solution at a rate greater than 40 
ml/min. At least twice as much oxygen was 
absorbed by S-Z (blocks 5 and 10) as by S-X 
because of the delay in clotting, which resulted 
from a decrease in pyrogallol content, 

It was previously stated that in solutions 
prepared by the old method no significant differ- 
ence was found in the amount of carbon monoxide 
evolved during oxygen absorptions when different 
However, when 

the improved 


platinum-disk pipettes were used. 

the solutions were prepared by 
method, the data shown in figure 5, blocks 1, 2, 5, 
and 10, where pipettes SG—1 and SG-—2 were used, 
indicate that in the initial absorptions of 80 per- 
cent of oxygen, 0.012 to 0.017 percent of carbon 
monoxide was formed, and in all the other pipettes 
the corresponding absorptions vielded only 0.004 
All of the 
sintered-glass bubbler pipettes were supposed to 
be identical, but pipettes SG—1 and SG-2 were 
purchased at a different time from the others and 


to 0.008 percent of carbon monoxide. 


there may have been some difference in the 
sintered-glass bubblers. 

Two or three drops of an antifoamer, 2-ethyl 
hexanol, were added to the solutions represented 
in blocks 3, 5, 7, and 9. No undesirable effects 
were observed as a result of its use, and it reduced 
the excessive foaming present when the solutions 


were new and of relatively low viscosity. 
(c) Age 


The age of a solution may refer to its age in 
time or to the extent to which it is expended. It 
may be seen in figure 6 that the aging of a solution 
before oxygen absorptions were begun had no 
effect upon the amount of carbon 
monoxide produced. 


significant 
The graphs on the left 
represent S-X, that prepared 
according to the old method in a plain bubbler 


solutions, were 


pipette. No significant variation in carbon mon- 
oxide formed is noted from 1 to 8 days. The 
erratic result at 6 days is an illustration of the 
vagaries to be expected when solutions are pre- 
pared from the same formula by the old method, 
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The plots on the right-hand side of figure 6 
show more conclusively that age has no effect. 
Two liters of solution B7 was prepared in an 
Erlenmeyer flask by the old method. The flask, 
on which the glass cap was sealed with paraffin, 
was stored out of direct sunlight for 17 weeks. 
About 150 ml of the solution was then transferred 
to a plain bubbler pipette, and the flask was again 
sealed and stored. This procedure was repeated 
when the age of the solution was that indicated 
in the graphs. Air analyses were made with each 
portion of this solution, and it was found that no 
significant change had occurred in the amount of 
carbon monoxide produced by solutions that were 
stored from 17 weeks to 21 months. 

The relationship between the carbon monoxide 
produced and the total volume of oxygen absorbed 
at the point of analysis is shown in figure 5. At 
first there is a gradual increase in the amount of 
carbon monoxide formed as more oxygen is ab- 
sorbed, and, in most cases, there is a rather sharp 
increase in the amount of carbon monoxide pro- 
duced after about 800 ml of oxygen absorbed by 
S-X and after about 2,500 ml of oxygen absorbed 
by S-Z. 


2. Effect of Partial Pressure of Oxygen 


The graphs of figure 7 show the amount of car- 
bon monoxide produced by five partial pressures 
of oxygen when absorbed by pyrogallol solutions 
prepared according to the old method in three 
types of pipettes.‘ Since the old method of pre- 
paring soluiions could not be relied upon to give 
consistent results from one run to another, these 
graphs may not show exactly the effect of the 
type of pipette upon the amount of carbon mon- 
oxide produced, although the Orsat contact pipette 
does show the characteristically high production 
of carbon monoxide. The graphs may be used to 
indicate the relative effects of several partial pres- 
sures of oxygen ranging from 21 to 100 percent of 
oxygen at various stages in the life of the solution. 
The broken lines in figure 7 separate the points 
into groups according to the partial pressure of 
oxygen absorbed. They show more clearly the 
fan-shaped effect that results from plotting 
amounts of carbon monoxide produced during 
absorptions of various partial pressures of oxygen 


versus the degree to which the solution was ex- 


? The letters in the upper left-hand corner indicate the type of pipette used 
For the key to abbrevatiations, see the legend of figure 2 
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pended. The same effect was shown in figure 5 
although it was less well defined with respect to 
the total amount of oxygen absorbed in the solu- 
tions prepared by the improved method. 

In figure 8 is a comparison of the carbon mon- 
oxide produced in the absorption of various partia! 
pressures of oxygen by a solution prepared by th 
old method and one prepared by the improved 
method. In the latter, the amount of carbon 
monoxide produced was directly proportional to 
the partial pressure of oxygen absorbed. In the 
solution prepared by the old method the carbon 
monoxide produced increased with the partial 
pressure of oxygen, but there was not the simple 
linear relationship that was observed with the im- 
proved solution. 

Since the volumes of the samples used were con- 
stant, the amounts of carbon monoxide formed 
were directly proportional to the total amount of 
oxygen absorbed and independent of the partial 
pressures of oxygen in the samples. This permits 
a close estimate of the carbon monoxide yield and, 
therefore, the proper correction for any determined 
analytical condition. 


3. Effect of Temperature of Reaction 


Figure 9 illustrates the effect of temperature 
variations in four different solutions that were pre- 
pared according to the improved method. The 
numbers next to the symbols in these graphs indi- 
cate the order in which the analyses were per- 
formed. A variation of 7 or 8 deg. C in the 
temperature of reaction was necessary to change 
significantly the amount of carbon monoxide 
formed in the absorption of oxygen by pyrogallol 
solutions. There was a definite relationship be- 
tween the temperature of reaction and the amount 
of carbon monoxide formed, although under most 
laboratory conditions the effect of temperature 
changes would be negligible. 


4. Effect of Rate of Flow 


The formation of carbon monoxide was signifi- 
cantly increased only when the rate of flow of the 
original sample was less than 20 ml/min, at which 
rate about 5 minutes was required to make one 
pass into the solution. So slow a rate would be 
considered impractical for most volumetric gas 
analyses. Figure 10 shows the amount of carbon 
monoxide formed at various rates of flow. 
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oxzygen into pipette SG-1 containing solution S—X. 


lhe numbers next to the symbols indicate the order in which the analy 


e performed 


It will be shown later in section 6 that the con 
tact type of pipette caused a much greater produc- 
tion of carbon monoxide than did the bubble: 
tvpe. This indicated that carbon monoxide was 
formed when high partial pressures of oxygen cam 
in contact with expended pyrogallol solution on 
the glass walls. If this assumption were true, 
would follow that a reduction in the rate of flow 
of the sample would cause a greater production o! 
carbon monoxide in the bubbler type of pipett 
because the oxygen would be in contact longe: 


~te28 5 5 6 5 £5 3s BD 











with the expended solution adhering to the inlet 
capillary tube leading to the bubbler. This seems 


Ficure 9 Percentage of carbon monoxide produced at ; ; 
to be a logical explanation of the data present: 


various té m peratures of reaction during the absorption 


from mixtures containing SU percent of ol ygen in figure 10. 
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5. Effect of Number of Passes Into Pyrogallol 


The efficiency of oxygen absorption in the 


itinum-disk bubbler pipette has been previously 


termined [5]. In the work referred to it was 
ind that in the analysis of a mixture containing 
percent of oxygen, more than 70 percent of 
total oxvgen in the sample was removed during 
first pass and about 93 percent had been 
moved after the second pass lt may therefore 
assumed that any significant amount of carbon 
onoxide was formed in the first two passes, since, 
” any one reagent at one temperature, the 
irtial pressure of oxvgen is the most important 
ctor in the formation of carbon monoxide 
Che graphs on the left of figure 11 show the amount 


of oxygen absorbed during the number of passes 


PD without bubbler 


cant amount of carbon monoxide after the second 
pass when about 6 percent of oxygen remained in 
When a 


percent of oxygen was absorbed by this solution, 


the residue. mixture containing 80 
erratic results were obtained in the two analyses 
made with five passes. This may be explained 
by the fact that about 1,200 ml of oxygen had 
been previously absorbed and the breaking point 
in the absorption curve had been reached, causing 
less efficient absorption of oxygen as well as a 


larger vield of carbon monoxide 


6. Effect of the Design of the Absorption Pipette 


For the determination of the effect of several 
tvpes of absorption pipettes upon the amount of 
carbon monoxide produced, the improved method 











4 








Beef 








" rm 





i. 
24 26 26 32 34 MM 8 40 42 446 SB DO 
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reentage of 


carbon monoxide produced and oxygen abse 


co PRODUCED, % 


ng varying numbers of} passes of sample into 


p fjrogatulot, 


ndicated The graphs on the right show the 
amount of carbon monoxide formed during those 
absorptions In order to exaggerate the effect, 
bubbler was not used and the sample was passed 
nto the pipette on top of the solution. The 
eakly alkaline solution S—5 in an absorption from 
ras contaming SO percent of oxvgen absorbed 
bout 10 percent of oxygen on each of the first 
hree passes, and the amount of carbon monoxide 
formed was almost directly proportional to the 
vvgen absorbed. Furthermore, absorptions of 
gen in air made with this solution showed no 
vnificant amount of carbon monoxide produced 
ter the fourth pass when only about 2 percent 
| oxvgen remained in the residue. 
In the absorption of oxygen from air, the 
ongly alkaline solution, S-7, vielded no signifi- 


Carbon Monoxide From Alkaline Pyrogallol 


e@. 


was used in preparing 2 liters of solution S-Z in 
an Erlenmeyer flask. Two hundred ml of this 
solution was siphoned, in an atmosphere of nitro- 
gen, into each of the pine pipettes illustrated in 
figure 2, and absorptions from mixtures of 21 and 
60 percent of oxvgen were made in each of the 
pipettes. The carbon monoxide evolved is shown 
in figure 12. The Orsat contact pipette was the 
only one that showed any great deviation from 
the others in the amount of carbon monoxide 
The sintered-glass bubbler, Branham- 
bubbler 


pipettes produced the smallest amounts of carbon 


produced. 


Sperling bubbler, and platinum -disk 


monoxide. The Dennis-Friedrichs pipette may 
also be included in this group with the note that 
the carbon monoxide formed might have been 
slightly less, except that the rate of flow of sample 
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Figure 12. 
the absorption of oxygen by solution S-BZ in various types 
of absorption pipettes (see fig. 2). 


Length of stem of bubbler: SG, 18.5 em: BS, 19 em: PD, 20em: DF, 19 em: 
FV,i5em; FA,l4em; PB,13em; T,10.5em , 20% On (air); @. 60% Osin Ny 
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was reduced to about one-half of the usual rat: 
This had to be done to prevent the sample fror 
passing to the atmospheric arm of the pipett: 
Figure 12 shows, in general, that the amount o 
carbon monoxide formed is a function of th 
depth of the pyrogallol solution above the bubbk 
and/or the size of the gas bubbles. 

It appears that there may be a correlation bi 
tween the formation of carbon monoxide and 
the depth of pyrogallol solution through which 
the bubbles travel. In the first four pipettes 
the distance from the bubbler tip to the top of the 
solution was 18.5 cm or more, and the carbon 
monoxide formed in the absorption from a mixture 
containing 60 percent of oxygen was 0.006 to 
0.009 percent; in the next four pipettes where the 
bubbles traveled from 10.5 to 15 em in solution, the 
carbon monoxide produced during similar absorp- 
tions was 0.010 to 0.014 percent. When the depth 
of the solution above the bubbler was decreased 
the effect became more pronounced, therefore, it 
appears to be desirable to have at least 18 em of 
solution above the bubbler. 

For a minimum production of carbon monoxide 
it is probably also desirable to distribute the 
sample in fine bubbles. The sintered-glass bub- 
bler, Branham-Sperling bubbler, and platinum- 
disk bubbler pipettes, which produced the least 
carbon monoxide, distribute the gas in many 
very fine bubbles. The Fisher valve-bubbler 
and Francis auto-bubbler pipettes produce bubbles 
of an intermediate size, and the plain bubbler 
pipette has the largest bubbles. The Dennis- 
Friedrichs pipette is not comparable with the 
others because of the long spiral up which the 
bubbles travel, and the Tramm pipette combines 
the contact method of absorption with the bub- 
bling method. The large amount of carbon monox- 
ide formed in the Orsat contact pipette indicates 
the importance of distributing the sample in such 
a way as to prevent its contacting excessive 
amounts of expended solution. 

Shepherd has shown that distributors of fine 
bubbles are the most efficient for the absorption of 
oxygen [5]. However, since the bubble size may 
not affect the carbon monoxide formed by an 
amount significant to many analyses, some ma) 
find it desirable to avoid the trouble caused from 
plugged distributors by using the plain bubble: 
with a larger number of passes for the complet 
absorpt ion of oxygen. 
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VII. Conclusions 


(he accuracy, reproducibility, and rapidity of 
analytical determination of oxygen by absorp- 
1 in alkaline pyrogallol may be improved as 
icated by this experimental study. 
1) Solutions can be prepared that evolve no 
nificant amount of carbon monoxide during the 
ial analysis. This is accomplished by the direct 
dition of a saturated solution of potassium 
vdroxide to pyrogallol crystals under conditions 
wh that the heat evolved during the solution is 
lissipated rapidly. 
2) Samples having partial pressures of oxygen 


higher than 50 percent can be analyzed with the 


mproved solutions without dilution of the sample 

3) The solution should not be used for analysis 
after the absorption of 12 volumes of oxygen per 
volume of solution if high accuracy is desired. 

1) The temperature of the pyrogallol reagent 
should not exceed 30° C during absorption. 

5) The rate of flow of sample into the pyrogallol 
should be as rapid as possible and never less than 


20 ml/min. 


‘arbon Monoxide From Alkaline Pyrogallol 


(6) A pipette that disperses the gas in fine 
bubbles through at least 18 cm of solution is 
preferable. 

It should be emphasized that if the analyst is 
concerned with amounts of carbon monoxide that 
are less than about 0.01 percent of the sample, a 
determination of the carbon monoxide formed 


should be made after each absorption of oxygen. 


The author acknowledges with gratitude the 
helpful suggestions and advice contributed by 
Martin Shepherd during this investigation. 
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Tracing of Skew Rays 
By Donald P. Feder 


This paper presents a set of convenient recurrence formulas for the tracing of skew rays 


by means of a desk calculator 


Because of the increased speed and field coverage of modern 


lenses, it is frequently desirable to trace such rays in making a complete study of an optical 


svstem 


The formulas apply to systems of spherical refracting surfaces 


be centered 
Particularly convenient check formulas 


the final surface 


\ complete study of an optical system appears 
to demand some knowledge of the behavior of the 
general or skew A computational process 
for obtaining this information by tracing such rays, 


rays. 


which has proved convenient in practice, will now 
he described. 

It is contemplated that these formulas will be 
ised with a desk calculating machine. When such 
a calculator is used, twelve numbers per surface 
must be written down for each ray. 

An advantage of the method is the existence of 
a simple check that need be applied only at the 
final surface. 

Brief 
spherical 


of 


need 


consider a 


The 


We 


refracting surfaces. 


derivation system 
system 
not be centered 

We denote the incident and refracted rays by 
the unit vectors Q and Q’, and the radius vector 
to the center of curvature from the point of 
R Bold-faced type will be used to 
enote vectors 

We define a new vector P 


cidence by 


P=N Q@XxR 


The law of refraction becomes 


P—N’ QR, 


N and N’ are the refractive indices of the 
irrounding media 
We can now choose our independent variables 


m the set of vectors (P, Q, R). Any two of 


here 


lracing of Skew Rays 


are 


The system need not 


available, which need be applied only at 


these determine the direction and location of the 


ray and may be used as independent variables. 
An obvious choice is the pair Q@ and R, but Pand Q 


vield simpler equations, particularly a simpler 


transfer. 

We eliminate R from eq 1 and 2 in order to obtain 
Q’ asa function of Pand @. Multiplying eq 1 vee- 
torially by Q, one can obtain R, which is then 
substituted in eq 2. After solving the resulting 
equation for Q’, we obtain 

Q’=—aQ+bP«Q, 
where 
P PP NN'(Q’-R)(Q-R) 
NN'UR 


N’'(Q’-R)—N(Q-R 
NN'UR 


and (Q’-R) may be evaluated by squaring 
We find 


(Q-R 
eq l and 2 respectively. 
P- 


N7(Q-RP = NPR 


N”"(Q’-R)?=N’ R?-—P 

Given P and Q, these equations enable us to 
find Q’=@Q It remains only to find P, for the 
following surface in terms of P and @ at the 
preceding one 

In figure 1 O and OQ, are the respective centers, 
and A and A, are the points of incidence of the 
ray. In general the quadrilateral is not plane 
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The signs of g and g’ are identical and are pos 
tive for surfaces convex toward the entering ray 





Equations 10 through 16 inclusive, form the 
complete set for systems of spherical surfaces. 

One notices that P, is an invariant of the ray 
P,=0 implies that the ray lies in a meridian plan 

Plane surfaces. At a plane refracting surfac: 
normal to the optical axis, the refraction equa- 
tions are greatly simplified. 


From eq | and 2 by taking R=i one obtains 
FiGtre | 


R, Q;, and R, denote vectors. N ) N ) 


on ee ee ee 


Denoting the vector from O to O, by T. one N’Z’ NZ | 
has immediately win x’ a1-y"_-z". 


R,+AA,@,—R+T () 


and .\’ is always positive. 

Multiplying vectorially by V’Q’=.N,Q,, the term Since P does not define the location “ ws , 
in AA: vanishes. and we have at a plane, we introduce a vector F, lying in the 
plane and extending from the axis to the inter- 


P,=P+@Q,> section of plane and ray. 


where t=N,T P 
bd | 
P cy 


| TOR <P. enact w Reape Go” 





Summary in scalar form. Since the vast ma- 
jority of optical systems are centered, let us sup- 
pose that OO, lies on the optical axis. We may 
take the zs axis along this line so that t,=?,=0. 

Suppose that @=CY, ¥, Z) and P=(P,, P,, P) 
are given in scalar form. Defining the auxiliary 
constants r= N|R) and r’=N’ RR}, and writing 


DM, Shae R eae Be? 





OG 9 et A RITTER RRC 


(10) 


(11) 





Ficure 2. 


R, GQ, F, and T denote vectors 


From figure 2 we obtain AA,Q,=R+ T+ F and 
multiplying by N,Q, as before, we have 


P=N\FXQ,+txQ, 


ee SOR A Lerwes ve 


where f=N,T. Inscalar form, if we take t, 
P,=N,|F,Z—F.Y}| 
P,=N,FL.X,—tZ, 

aZ+6|P,Y—P,.X] P,=—-N,F,X,+tY 


aX+b[P,Z 


ay +b[(P,.X—P,Z| 
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Vara = 
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+ See Eten 





dently eq 20 may be used to obtain F in terms 
P or conversely. Frequently they are used as 
ial and final equations as well as at refracting 
races 


hecking equations At each surface 


P-Q—P-Q, 0| 
Q@=—-@=1 | 


(2] 


meqlor2. It is important to note that these 
cks need be applied only at the final surface. 


For we have from eq 3 by scalar multiplication 


P-Q’ —aP-Q 
d similarly from eq 9 
P,-Q,= P-Q,=aP-Q. (29) 


} 


Furthermore, since @ can be shown to be always 
onzero, it follows that if P-@ is nonzero at any 


nt, it remains nonzero. This implies that an 





error that causes P-Q to differ from zero at any 

surface will appear at the final surface. 
Similarly, if any error is made so that Q 

the error will appear at the final surface 


140, 


In addition one must have at each surface 
a- i?P* l 0) 


P—P—P—P=0} 

These checks should guarantee the accuracy of 
the final data, provided no errors are made in the 
constants. 

Using this method, the accuracy of the results 
is not limited by a particular set of tables All 
that is needed is two square roots, and this may 
be obtained, with any desired accuracy, by the 
usual iteration formulas. It is very easy to ob- 


tain results to 10 significant figures by emploving 
at most a small table 


WASHINGTON, Mareh 30, 1950. 
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Strength of Heat-Resistant Laminated Plastics up to 
300° C* 


By Benjamin M. Axilrod and Martha A. Sherman 


The flexural properties of samples of glass-fabric laminates were determined for several 
conditions of heating and testing. The laminates tested were bonded with various resins, 
including unsaturated-polyester, acrylic, silicone, phenolic, and melamine types. Flexural 
tests were made under the following conditions: (1) at 25° C and 50-percent relative humidity 
after 200 hours of heating at a temperature 7; (2) at a temperature T after 0.5 hour at the 
temperature 7; (3) at a temperature 7’ after 200 hours at the temperature 7; and (4) at 
25° C and 50-percent relative humidity without heating prior to testing. The temperature 
T was 150°, 200°, and 250° C for condition 1; 150°, 200°, 250°, 300°, and 375° C for condi- 
tion 2; and 150°, 200°, 250°, and 325° C for condition 3 

The one sample of silicone resin laminate tested was superior to the other laminates in 
retention of flexural properties at temperatures of 250° C or higher. The single phenolic 
laminate tested showed good retention of flexural properties when tested at elevated tem- 
peratures after 0.5 hour of heating After prolonged heating the behavior of the phenolic 
sample was much less satisfactory, the strength being almost negligible at 250° C. The 
sample of melamine laminate was superior to the phenolic in retention of flexural strength 
after prolonged heating at 200° and 250° C. The polyester laminates lost at least 80 per- 
cent of their flexural strength when tested at 250° C. For each laminate and each test con- 
dition loss in modulus of elasticity correlated with loss in flexural strength, but the former 
loss was less than the latter. The results of one test in which duplicate sets of specimens 
were heated in still air and in rapidly moving air indicate that the loss of flexural properties 


during prolonged heating may depend on the method of exposing the specimens 


I. Introduction mation on the mechanical properties of laminated 
-_ : , oe plastics at temperatures much higher than 100° C, 
lhe mechanical properties of glass-fabric lam- 


inates at elevated temperatures are of interest to The present investigation was undertaken to 


aircraft designers, as these materials have prop- determine the flexural properties of promising 


: ‘ : ‘ inates at ’ rt ‘e 200° C C = 
erties that can be used to advantage in certain laminates at temperatures up to 3 C. = 


applications, such asradomes. Considerable work mercial resin manufacturers and laminators were 


had been done prior to 1947 on the mechanical informed of the project and invited to submit 


properties of various laminates at temperatures samples of laminates that could be expected to be 
in the range of 70° to 100° C [1 to 7],' and a few 


is were reported at temperatures of 200° C [8]. 


heat-resistant. Glass-fabric laminates bonded 
with various resins including silicone, phenolic, 
\\ th the development of higher speed aircraft melamine, acrylic, and unsaturated-polyester 
ind guided missiles. there was a demand for infor- types were rec ived ind tested. This report sum 
a marizes the results of the flexural tests, made at 
nvestigation was conducted under the sponsorship and with the - . oe —— ft ~ ; | = | | an en | 
inlidiaain Gk Qn Veit Abeta Guntiiiies te Beneath room temperature atiter prolonge eating ane 
\ir Matériel Command also at elevated temperatures after both brief and 
n brackets indic th iterature references at the end of thi 
prolonged heating. 
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As testing at temperatures of 300° C presents numbered. The numbered specimens were 


problems not encountered in testing at room tem- ranged in a random order and separated in 


perature or at 70° to 100° C, a considerable part groups of seven each. The plan was to t 


of the project was devoted to the development five specimens from each group with two sp« 
of suitable equipment and procedures. These mens for use in case of experimental difficulti 
are described in detail in this report The statistical analysis of the results can 

greatly facilitated by having the same num! 


II. Materials of test specimens for each condition. 
The flexural specimens were 4 in. long a 


The materials, submitted in sheets approxi- . , 
were machined to a width of 1.00 in., either by d 


mately ', in. in thickness, are identified in this Ag ‘ , wes 
erinding with a surface grinder or milling with 


report as follows 7 ; 
carbide-tipped cutter, 


2. Flexural Test Equipment ” 


Laminate samp tesin designation 


The flexural jig and the heated test enclosure a 
Unsaturated polyester | , . : oy : 
: Pa shown in figures 1 and 2 in place in the testi 
1 Laminae PDL-7- 261 2 4 . pes ° 
bk Plaskon 911-11 machine. The flexural jig and test enclosure are 
hs ) . : 
4 —— wo 1OA. ¢ both suspended from the upper cross-head of th 
araplex 3. : o- a 
Arerylic additior testing machine. The flexural jig, partly dis- 
( BCM monomer and poly assembled, is shown in figure 3. The base of th 
methyl methaery 
late 
Silicone bloc} t | | i , i : 
I \ re DO -, ) -10 »§ ) 
p- er ks may be bolted in place to give a span of 2, 
Phenolic 
i Redux 
Melamine 
] 


flexural jig has slots located so that the support 


Melmae 405 


All the laminates except the sample J were made 
with Fiberglas 181A fabric; Fiberglas ECC—128 
was used in sample J. The materials are de- 
scribed in detail, including fabrication conditions, 


in table 1 


Ill. Test Procedures 
1. Flexural Specimens and Sampling 


It is known that the sheet-to-sheet variability, 
especially in experimental laminates, may be large 
A proper sampling plan for the tests at various 
conditions could give information both as to the 
average behavior of the sheets of a given sample, 
and also the variations that might be expected 
both at standard test conditions and elevated 
temperatures. Since both the number of sheets 
and the sizes differed from sample to sample, it 


. Figure 1. Front view of flexural test enclos in plac 
was not practical to study sheet-to-sheet vari- pasaien dn de ota “ah veel 


a the te sting machine 
ability. :' 
. J . he recording deflectometer and selsyn motor are at the left below 
The following sampling plan was adopted: ontiaiiiad 
Approximately 80 percent of each sheet of each preg ERAT 
; ; : LL: ? Drawings of the flexural test equipment are available from the Ory 
sample was cut into specimens, which were lastics Section, National Bureau of Standards, Washington 25, D 
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Desc iplion of materials 


Sample A 


American Cyanamid ( 


Paraplex P-4 


Unsaturated p 


Thalid X-52s8-10A ¢ 


Plaskon 911 
saturated polyes 


Polyester 


Benzoyl! peroxide 


styrene 


Hydraulic pre 


Non 


i heating und 


ooling 
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TABLE 1. Description of materials—Continued 
Sample C Sample B Sample H Sample J 
Laminator of sample E. I. du Pont de Nemours & | Dow Corning Corp. Rohm & Haas Co Formica Co 
Co 
Resin used 
Trade name and number BCM monomer, 60%; poly DC 804 Redux Melmae 405 
methylmethacrylate HG—4, 
wr, 
Chemical type Acrylic addition type; thermo- | Silicone Phenolic Melamine-forma 
setting hyde 
Catalyst used with resin 
Name or type Residual benzoyl peroxide in None None 
polymethylmethacrylate 
Amount, ‘ 
Liquid added to resin to prepare it for application to fabric 
Name Methy] ethyl ketone roluene None Water, 95%; | 
ilcohol, 5+ 
Amount, parts per 100 parts of resin 135 67 67 
Resin content (by weight) of finished panel 
Manufacturer's data, i) 9 4) to 45 42 to 47 
44 ) 


National Bureau of Standards data, ‘ 62 
Fabric used 


Identity Fiberglas ECC-128 


sizing pyrolyzed t 


Fiberglas 181-38, 
Finish 114 


Fiberglas 181-A-13 Fiberglas 181-A-12 


light brown 


Number of plies 6 11 12 13 

Ply arrangement Cross Cross Cross Parallel : 

Treatment prior to application of resin (drying, etc None None : 

Method of applying resin Continuous dip and squeeze | Dipping Dipping Dipping : 
roll 

Drying of impregnated fabric Solvent evaporated in forced Air dried for 30 min 





Molding conditions 
Equipment 


draft oven at 43° C 


it room tempera- 
ture, 15 to 18 min 
at 104° © in forced 
draft oven 


Heat Steam Steam Steam Steam 
Pressure Hydraulic press Hydraulic press Hydraulic press Hydraulic press { 
Precure None None 
Temperature, °C 110 145 
Pressure, Ib/in.? Atmospheric 15 : 
Time, min 10 1 : 
. : 
Cure : . 125 to 124, 20 min 
Temperature, °¢ 110 175 145 aa 
(142 to 147, 20 min 
Pressure, Ib/in.? NM) 00 1) 1,100 
Time of heating under pressure, min 15 75 15 About 40 min 
Time of cooling under pressure, min } 10 None 15. 
(90° C, 15hr* 
j! ww? C,4hbr 
Aftercure None 10° C,4hr do None. 
[1 C, 6 hr 
200° C, 4 hr j 
Sample 
Area, in 14 by 19 11 by 12 12 by 12 10 by 10 21 by 36 
rhickness, in 0.10 to 0.13 0.10 to 0.13 0.13 to 0.14 0.13 to 0.15 0.12 to 0.13. 
Number of sheets 3 2 } i4 1 


>» Mechanical convection oven 
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FiGuRE 2. 


i, Loading stirrup: B, loading edge in loading stirrup: C, contact edges in st 


ray: //, stem of dial-type thermometer 


t, or 5 in. The loading stirrup, attached to 
the lower cross-head of the testing machine 


rod used for the loading edge extend into slots 





FicurE 3 Flexural J partly disassembled. 


left support block and clamp and the loading stirrup have 


land rotated through 90 
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Interior view of test enclosure. 


bee 


fig. 1), is centered by having the ends of the drill 


ipport block; D, loading port; F, hole for push rod; F, specimen in jig; G, speci 


in the flexural jig suspension (see figs. 2 and 3 
The flexural jig is made of steel, nickel-plated to 
avoid corrosion. 

The deflection of the specimen is measured by 
a Southwark-Peters deflectometer, model PD-—1, 
coupled to a recorder on the testing machine to 
produce load-deflection graphs. The deflectom- 
eter, located below the test enclosure (fig. 1), is 
supported by a platform attached to the flexural 
jig base. A glass push rod touching the lower 
side of the specimen operates the deflectometer; 
this rod is centered by a guide in the loading 
stirrup (fig. 3 

The test enclosure has an asbestos board 
exterior, a sheet aluminum interior, and Fiberglas 
insulation about 1.75 in. thick. The interior is 
divided into two compartments by an aluminum 
partition; the front compartment is the test 
chamber and the rear compartment contains a 
bank of five 250-w strip heaters. A fan operating 
continuously draws air past the heaters and blows 
it into the test chamber. After a number of tests 
had been made, it was found desirable to place 
two layers of '-in.-thick asbestos board in front 
of the partition (fig. 2) to improve the temperature 
distribution within the test chamber. 

The specimens are inserted with a spatula-like 
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holder through a loading port and then placed side 
by side on an aluminum specimen tray, visible at 
the right of figure 2. This tray holds five 1-in.-wide 
specimens. The specimens are manipulated into 
place in the flexural jig by push rods inserted into 
holes in the front and right side of the enclosure. 
In figure 2 one specimen is shown in place in the 
flexural jig and one in the specimen tray. After a 
test the specimen is removed from the flexural jig 
with a pair of tongs inserted through a port in the 
left side of the enclosure. The operation of the 
equipment is aided by a double-paned Pyrex glass 
window and a projection bulb located inside the 
enclosure. To avoid failure of the phenolic plastic 
socket or the bulb, the latter is operated at lower 
than rated voltage by means of a variable auto 
transformer, and the socket is attached to the 
outer wall of the enclosure. 

Preliminary temperature measurements were 
made with thermocouples located inside and on the 
surface of specimens in both the specimen tray and 
in the flexural jig. It was found that a specimen 
came to temperature equilibrium within 10 to 15 
min after insertion in the enclosure. Before the 
asbestos board was added to the partition, the 
temperature of a specimen in the tray was appre- 
ciably less than that of a specimen in the flexural 
jig, the differences amounting to 3 deg and 7 deg C 
at enclosure temperatures of 150° and 250° C, 
respectively. The difference was reduced to 3 deg 
C at 250° C by insertion of the asbestos board. 

A fine mesh screen (fig. 2) was wrapped around 
the upper part of the flexural jig to keep particles 
of Fiberglas or resin out of the guides for the 
loading stirrup. 

In testing at elevated temperatures the gaseous 
products given off by some of the specimens were 
a source of trouble, as the testing machine was 
located in a controlled atmosphere room in which 
the air was recirculated. The problem became 
more serious as the test temperature was increased. 
An exhaust pipe equipped with a blower was finally 
attached to the box; it was necessary to control 
the amount of air drawn from the enclosure by 
means of a valve in the exhaust pipe, as otherwise 
the temperature of the enclosure could not be 
maintained constant. 
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3. Flexural Test Procedure 


(a) Temperature Conditions 


Sets of specimens of the samples were tested at 
25° C after heating at temperatures of 150°, 200 
and 250° C for 200 hr. Other sets of specimen: 
were tested at an elevated temperature 7 afte 
0.5 br or 200 hr at the temperature 7. The tem 
peratures ranged up to 375° C in the 0.5-hr test: 
and 325° C in the 200-hr tests. Each sample wa: 
not tested at each set of conditions; in some in 
stances the samples were badly weakened at mod 
erate test conditions, and in other instances suffi 
cient material was not available. 

Two control sets of specimens of each sampk 


> 


were tested at 25° C and 50-percent relativ: 
humidity. One set was tested at the beginning 
of the testing program, the other set about 7 mo 
later, near the end of the tests. 

In heating the specimens for 200 hr, a cireu 
lating-air oven was used with each specimen in 
an open wide-mouthed glass bottle of 2-in. diam 
eter and 4-in. height. In one experiment speci 
mens were also placed on asbestos board laid on 
the oveo shelf in order to compare the effect of 
heating in still air with heating in rapidly moving 
air. For convenience, the 200-hr heating period 
was interrupted at about 96 hr, and the speci- 
mens were allowed to cool in the glass containers 
for about 18 hr; the heating was then continued 
for about 103 hr. The specimens were trans- 
ferred to the heated test enclosure 1 hr prior to 
testing; in making the transfer the specimens 
were transported to the immediate vicinity of the 
test enclosure in their containers to minimize 
cooling; transfer from the container into the test 
enclosure took a few seconds. During all the 
200-hr heating periods the oven temperature was 
recorded with a Brown recording potentiomete! 
and a thermocouple. 

In the tests at elevated temperatures after the 
specimens were heated 0.5 hr, it was considered 
that the thermal shock when a cold specimen 
was placed on the hot metal specimen tray might 
weaken the specimens. Accordingly, in som: 
tests the specimen tray was covered with a shee! 
of \\e-in.-thick asbestos paper. 
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25° C and 
specimens 


In tests at 
nidity, the 


50-percent relative 


were normally con- 


oned for 7 
osphere was controlled at the above values. 


days prior to test in a room whose 


(b) Flexural Test Conditions 
The 


tual span was determined at 25° C with screw 


lhe flexural jig span was set for 2 in. 
crometer calipers. The contact edges of the 

pporting and loading pieces were all of \-in.- 
diameter drill 
\intained at 0.04 in./min in all tests. 


rod. The rate of loading was 


(c) Calculation of Flexural Data 


In calculations of tests at elevated temperatures 
the specimen dimensions obtained at 25° C and 
50-percent relative humidity were used; however, 
the span was corrected for change with tempera- 
io 


The stress 


ture, assuming a thermal expansivity of 12> 
per deg C for the metal in the jig. 
range chosen for calculating the flexural secant 
modulus of elasticity values was about one-half 
to two-thirds of the flexural strength at the test 


condition concerned. 


(d) Resin Content and Loss of Weight on Heating 


The resin content was determined for three 
sets of flexural specimens, except for the laminate 
made with the silicone resin which could not be 
completely burned off. The initial weights were 
determined after the specimens were conditioned 
for at least 7 days at 25° C and 50-percent relative 


Strength of Laminates up to 300° C. 
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humidity. The resins were removed from the 
specimens by heating in a muffle furnace. 

The weight loss after 200 hr of heating was 
determined by weighing some of the sets of speci- 
mens prior to heating and reweighing after the 
specimens were reconditioned 7 days at 25° C 


and 50-percent relative humidity. 
(e) Tolerances on Temperatures and Times of Heating 


The temperature ia the circulating-air oven was 
controlled at the nominal value to within +5° C 
during the 200-hr heating period. In the heating 
at 325° C in a muffle furnace, the temperature 
fluctuations were +10 deg C relative to 325° C. 

In the test chamber the average temperature 
for each set of five specimens noted at the speci- 
men at the start of the flexural test was within 

+3 deg C of the test temperature reported in 
tables 2 and 3 with only one exception. The 
latter occurred in testing sample B at approxi- 
mately 375° C after heating for 0.5 hr, the average 
test temperature being 368° C. In more than 90 
percent of the tests the average test temperature 


was within 2 deg C of the desired value. The 
individual test temperatures were within —5 deg 


and +3 deg C of the nominal values except for 


the tests at 375° C for which the range was 367° 
to 381° C, 

In the tests at a temperature T after 0.5 hr at 
the temperature 7, the actual time of heating at 
the start of the flexural tests was 30 to 35 min. 
In the 200-hr tests, the total time of heating was 


200+ 1 hr. 
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‘lerural strength data for 


Special conditions 


FLEXURAL STRENGTH OF ¢ 


Oct. 1048 
Apr. 1049 


Control mean 


A. First set 
B 


Second set 


FLEXURAL STRENGTH OF 


Not reconditioned ¢ 


Not heated in bottles 4 


FLEXURAL STRENGTH 


A bestos-covered tray * 


Asbestos-covered tray * 
do 


FLEXURAL 


See footnotes at end of table 


ONTROL 


25.0 to 


SPF 


OF 


STRENGTH OF 


specimens of glass-fabric laminates tested under various conditions of temperature « 


heating * 
Unsaturated polyester resin laminates 


A Laminate E Laminate F Laminate G 


iminate 


Per 
cent- 
age of 
control 
mean 


Per- 

cent 

we of 
control 
mean 


Per- 

cent 

age of 
control 
mean 


Mean Range Mean Range Mean Range 


AT 25° C (AFTER AT LEAST 7 DAYS AT 25° C AND 530% REI 


HUMIDITY 


SPECIMENS 


10 ih 1 ih 10 ih/in2 
in? 
49.8 47.1 to 52.4 


48.4 to 52.9 


1 Lb/in2 108 Th/in.2 1 Ib/in2 


‘7 
49 


1 to (4.4 
Sto 53.4 


4M. ito 20.9 M6 to 61.8 


$1.2 8 


58. 2 to 60.6 


AFTER 200 HR OF HEATING © (RECONDITIONED 7 DAYS AT 25 


MIDITY, UNLESS OTHERWISE NOTED 


CIMENS AT 25° C 
RELATIVE Ht 


40. 1 to 49. 
1 17.6 to 46. ¢ 
24. 4 to 30 


18 


51.0 to 8.1 48.1 to 51.7 9 43.6 
5 to 4.4 49.7 4.7 14 «41 
7 to 44.9 3) 
21.4 


40.6 
m@.2 Mi 
).8 a5 
4.4 to 40.5 

$to 10.6 


6 


Sto? 


v6 


SPECIMENS AT A TEMPERATURE AFTER 0.5 HR AT THE TEMPER 


11.0 to 12.8 ) 4 
ito 


18.0 to 21.2 12.0 


9.9toll4 r 6, &. 9 


7.lto &0 
7.4to &4 
5.2to 4.6 


5. 2to 45.9 


RE T AFTER 20 HR AT THE TEMPERATURE T* 


SPECIMENS AT TEMPERATI 


27.2 to 30.6 23.6 to 25.0 


17.1 to 19.4 


&. 6to 10.1 


16. 6 to 18 


3 
2Z21to 24 
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Flexural strength data for specimens of glass-fabric laminates tested under various conditions of temperature and 
heating *—Continued 





Acrylic laminate C Silicone laminate B Phenolic laminate H Melamine laminate J 
Special conditions Per- Per- Per Per- 
centage centage centage centage 
Mean Range of con- Mean Range of con- Mean Range of con- Mean Range of con- 
trol trol trol trol 
mean * mean mean mean 


FLEXURAL STRENGTH OF CONTROL SPECIMENS AT 25° C (AFTER AT EAST 7 DAYS AT 25°C AND 50% RELATIVE 
HUMIDITY 





10 10 10 103 
lb/in 108 lh/in2 lh/in lbj/in. 18 lh/in? thjin2 105 Ib/in 
A. First set (Oct. 1948 244.3 23.7 to 25.1 28.9 7 44.9 5.0 to 50.6 51.6 W.2to 8.4 
B. Second set (Apr. 1949 26.8 25.1 to 28.4 7.6 48.2 42.0 to 33.6 51.5 5.0 to 52.8 
Control mean 29. 2 46. 5 1.6 


FLEXURAL STRENGTH OF SPECIMENS AT 25°C AFTER 200 HR OF HEATING > (RECONDITIONED 7 DAYS AT 25° C 
AND 30% RELATIVE HUMIDITY, UNLESS OTHERWISE NOTED 





28. 2 26.8 to 29.4 116(A)) 27.7 26. 2 to 30.3 95 40. 1 6 Se 33.8 31.8 to 35.0 6 

Not reconditioned 1.7 27.6 to 34.9 130CA 2.0 2%. 6 to 29.1 OF 40.3 ; 7 11.6 1.3 to 32.1 6l 
15.0 8 32 26.4 5.7 to 28.4 51 

Not heated in bottl y 4 DD 2.0 25. 1 to 26.8 ”) 
26. 4 25. 2 to 27.7 ”) 0.2 17.6 16.8 to 18.4 ‘4 


FLEXURAL STRENGTH OF SPECIMENS AT A TEMPERATURE T AFTER 0.5 HR AT THE TEMPERATURE 7 





4/ 29to42 14(A)| 12.6 | 11.3 to 13.7 43 41.9 9.7 to 44.4 90 | 37.2 | 36.4to 381 72 
5.3 28. 2 to 40. ¢ 76 | 26.5 24.8 to 27.4 1 
10.9 10.1 to 11.9 7 | 2.9 21.5 to 28.7 | 22. 1 21.4 to 22.7 3 
Asbestos-covered tray ¢ 0.47 0.42 to 0. 51 1. 8(B 24.8 22.6 to 26.4 a 22. 4 21.8 to 23.5 45 
" 10. 1 9.6 to 10.9 85 19.8 17.4 to 21.7 43 20.5 19. 2 to 22.6 40 
: Asbestos-covered tray ° 10.5 | 10.3to 10.8 6 | 20.5! 17.0 to 25.3 4619.0 18 4t019.4 37 
do 9.8 ‘8.8 to 10.6 ‘34 14.7 12.3 to 16.4 2 2.0 l4to 2 ‘4 
FLEXURAL STRENGTH OF SPECIMENS AT A TEMPERATURE T AFTER 20 HR AT THE TEMPERATURE T° 
4.2 2.9 ¢t i 17(A 14 13.4 to 15.7 4 9. 1 36. 2 to 42.6 4 3 29. 4 to 32.3 sy 
10.9  10.4to 11.9 23 23.3 23.0 to 24.2 45 
13.4 12.7 to 14.0 4 0.55 O.5) to 0.59 1.2 15. ¢ 14.7 to 17.4 x0 
e 98 9.2t0 10.4 4 1.5 l.5to 1.6 i3 
* Specimens approximately 44 in. thick and 1 in. wide tested at 0.04 in./min Unless otherwise noted, all tests at elevated temperature were made with 
head motion with a span of 2.0in. Each value is the average for five the normal tray 
unless otherwise noted. Specimens taken at random from the f Specimens not randomly selected but taken from at least three sheets. 
nless otherwise noted « A muffle furnace controlled at 325° +10° C was used for conditioning these 
otherwise noted, all specimens that were heated for 200 hr were specimens 





open wide-mouthed bottles placed in a circulating-air oven b The flexural strengths for the first and second controls differ significantly 








. removal from the oven, specimens kept over calcium chloride for for this sample Ihe test results for various conditions were referred to the 
r to testing data for the control tested at nearly the same time; this is indicated by (A 
mens laid on asbestos board on oven shelf or (B) following the percentage 
normal tray is an aluminum tray he asbestos-covered tray is i Four specimens 

iluminum tray covered with asbestos paper about Me in. thic i Two specimens 
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IV. Results and Discussion 


The flexural strength data for the various 
samples tested under different conditions of tem- 
perature and heating are shown in table 2. Some 
of the data in table 2 are shown graphically in 
figures 4 to 6; the data that are plotted are for 
normal conditioning and testing, except when the 
only tests made were with the tray covered with 
The corresponding data for flexural 
secant modulus of elasticity are shown in table 


asbestos. 
3 and figures 7 to 9. The secant moduli selected 
for graphing are for a stress range from zero to one- 
half to two-thirds of the flexural strength; these 
secant moduli are expressed as percentages of the 
values for the controls for the corresponding stress 
range. The loss-of-weight data for the samples 
after heating for 200 hr at various temperatures 
are given in table 4. The precision of the data 
and other statistical points are discussed in the 
next section of this report. 

In the detailed discussion of the results that 
follows, it must be noted that, except for the 


polyester laminates, only one sample of a given 


tvpe laminate was tested. Hence only tentative 
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FicurReE 4 Flezural strength of the laminates at 25° C after 
200 hours of heating at various temperatures. 


A, E, F, G, Unsaturated polyester; B, silicone; C, acrylic addition; //, 
phenolic; J, melamine. 
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Figure 5. 
temperature T after 0.5 hour of heating at T. 


Flezural strength of glass-fabric laminates at a 


i, E, F, G, Unsaturated polyester; B, silicone; C, acrylic addition; H/ 
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Ficure 6. 
temperature T after 200 hours of heating at T. 


Flexural strength of glass-fabric laminates at 


A, E, F, G, Unsaturated polyester; B, silicone; C, acrylic addition 
phenolic; /, melamine. 
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ferences may be drawn regarding the behavior 
the various types of materials. 


1. Flexural Strength 


(he four polyester samples (figs. 4 to 6) show 
s than 20-percent loss in strength when tested at 

C after 200 hr of heating at 150° C. When 
ted at 25° C after 200 hr of heating at tem- 
strength of Laminates up to 300° C. 
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Figure 9. Flerural secant modulus of 
laminates at a temperature 
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4, E, 
phenolic; 


silicone, C, acrylic addition; /7/, 


F, G, 
J, melamine 


Unsaturated polyester; B, 


TABLE 4 


laminates after heating 200 hours at various temperatures* 


Loss of weight and loss of resin by glass-fabric 


N BS sampk 


designation Resin) Weight| Resin) Weight Resin) Weight, Resin 
: 


loss>| Joss |loss>| loss |loss*| loss | los 


saturated- 
polyester 
resin 
i 
I 
I 
G 
Acrylic addi 
tion resin 
« 
Silicone resin 
R 
Phenolic resin 
Hi 
Melamine res 
in 


* The initial weights were determined after at least 7 days’ conditioning 


at 25° C and 3% 
specimens had been reconditioned 7 days at 25° C 
200 hr of heating 
> Calculated on basis of NBS values for resin content except for laminate B 


relative humidity; the final weights were obtained after the 


and 50% relative humidity 


following the Each value is the average for 5 specimens 


the assumption is made that total weight loss is resin 


150° C, the 
increases rapidly, reaching about 90 percent at 
250° C. 
in strength ranging from 50 to 90 percent when 
150° C 200 hr of 


peratures above loss in strength 
These four polyester samples show losses 
tested at after either 0.5 or 


79 





heating at 150° C. For similar treatment at 
250° C the loss in strength is at least 80 percent, 
even for samples tested after only 0.5 hr of 
heating. 

The acrylic laminate C, when tested at 150° C 
after either 0.5 or 200 hr of heating, retained 
only about 15 percent of its initial strength. 

The flexural strength at 25° C of silicone 
laminate B was reduced only 10 percent after 
heating 200 hr at 250° C 
observed for any sample for this test condition, 
In tests at elevated temperatures, after heating 
the samples for either 0.5 or 200 hr, the flexural 
strength dropped considerably between 25° and 
150° C. The flexural strength at 150° C was 
about 13,000 Ib/in2 compared to a value of 
29,000 lb/in.? for the controls. Above 150° C the 
strength diminished very slowly with increase in 
temperature up to 250° C for the 200-hr and up to 
375° C for the 0.5-hr heating period, respectively. 
It is of interest to compare the results obtained in 
this laboratory with data obtained by the New 
York Naval Shipyard [9] on two samples of 0.5-in. 
thick laminate from two 


sources. 


; this was the least change 


glass-fabric silicone 
Some of these results, as well as data for 
sample B, are shown in table 5. 

Taste 5. Comparison of results for sample B with other 
flexural data on silicone laminates 


Flexural strength at T after heating at T 


Temper- 

ature of 
condition 
ing and test 


Sample J [9] Sample J/ [9] Sample B 
Heated Heated Heated Heated Heated Heated 
1 hr 192 tir 1 hr 192 hr 0.5 hr 200 hr 


10° lh/in ibe 2 10° Dh/in2 10° TA/in.2 10° Th/in2 
7.2 i 9.8 12.0 12.6 14.3 


12.8 5.0 2 10.9 13.4 


LO‘ ih/in 


Control val 


ues 


‘ 


Sample J was reported to be softened at 250° ¢ 
when heated only a short time at that tempera- 


ture; in view of the improvement on further heat- 


ing, this behavior was believed to be caused by 
further cure [9]. To test this idea they report [9] 
that specimens of sample / were heated for 72 hr 
at 250° C, reconditioned for 72 hr at 25° C and 
50-percent relative humidity, and then tested at 
250° C after 1 to 4 hr at 250° C; flexural strengths 
of 11,000 Ib/in.? were obtained. These results 
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give some indication of the variability in behavio 
of different silicone laminates and also show th; 
possibility with such laminates of having an under 
cured material. 

The sample of phenolic laminate /7 was out 
standing in retention of flexural strength at el 
vated temperatures when heated only 0.5 hr at 
the test temperature. For example, the flexura| 
strength at 375° C for this test condition was 
15,000 lb/in.? compared to 46,000 Ib/in for th: 
control (fig. 5). The same sample, however, was 
seriously degraded by prolonged heating at 250° C 
after such treatment it exhibited a flexural strength 
much less than 1,000 lb/in.? when tested either at 
250° or 25° C. 

The melamine laminate J shows a nearly steady 
decline in strength from 25° to 250° C or more for 
the three heating and testing conditions. Th: 
flexural strength is reduced to about two-thirds, 
one-half, and one-third for heating at 150°, 200°, 
and 250° C, respectively. Ata given temperature 
T, the strength at 7 after 0.5 hr of heating is 
slightly greater than that for tests at 25° C after 
200 hr of heating, and the latter strength is 
slightly greater than that at a temperature 7 
after 200 hr of heating. These flexural strength 
values were compared with data obtained by the 
New York Naval Shipyard [9]. The latter 
laboratory made tests similar to those described 
in this report on a 0.5-in.-thick melamine laminate. 
The general trend of the strength-temperature 
data given in [9] is similar to that shown in figures 
4 to 6 for sample J. 


2. Flexural Modulus of Elasticity 


For a given sample the trend shown by the 
flexural secant modulus of elasticity versus tem- 
perature graph is quite similar to the plot of flex- 
ural strength, versus temperature for the same 
However, for 
flexural 


condition of heating and testing. 
each laminate the percentage loss in 
modulus of elasticity for a given test condition is 
always less than the loss in flexural strength for 


that test condition.* 
3. Correlation Between Changes in Flexural 
Strength and Modulus of Elasticity 
In figure 10 the flexural strength in percentage 
of control is plotted against the flexural modulus of 
elasticity in percentage of control for all the lami- 


3 There were exceptions to this statement when the loss in the flexural prop 
erty was 5% or less. Such losses are not considered significant. 
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Data for all conditions in tables 
and 3 were used, except that the values for 


tes except C. 


Special conditions’ were not plotted unless no 
her data at a given temperature were available. 
he modulus of elasticity data for the lowest 
ress range was used in the case of tests at 150° C. 
t is evident that a correlation exists between these 
wo quantities; furthermore, it appears that the 
orrelation equation is probably different for the 
olyester laminates as compared to the melamine 
aminate J. A detailed analysis of the correlation 
was not made, as only in the case of the polyester 
iminates was there more than one sample of a 


given type of laminate. 
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SITY , PERCENT OF CONTROL 
FIGURE 10. Correlation between 


flexural strength and 


flexu ral secant modulus of ¢ lasticity. 


Data taken from tables 2 and 3 DF, w; G, @ B, A; 


H, @. 


Sample A, . Ez, 


4. Relative Heat Resistance of the Laminates 


To give a concise comparison of the laminates as 
regards their retention of flexural properties at 
elevated temperatures, a table was prepared in the 
following manner: For each sample the tempera- 
tures in figures 4 to 6 were noted for which the 
flexural strength was 30 percent or more of the 
flexural strength for the control for all three con- 
ditions of heating and testing. The same pro- 
edure was followed for the corresponding secant 
nodulus of elasticity values, except that the per- 
entage was 50 instead of 30. These temperature 
mits are shown in table 6. Since it was often 
ecessary to interpolate the stress range as well as 
he modulus of elasticity, the temperatures given 
n table 6 for the modulus data are less accurate 


Strength of Laminates up to 300° C. 


Relative heat resistance of glass-fabric laminates, 
based on flexural properties 


TABLE 6. 


remperature T at which flexural 
secant modulus of el 
at least 50°; of control value for 
the three types of test * 


Temperature 7° at 
which flexural 
strength is at least 
30%, of control 
value for the three 
types of tests * 


asticity is 


NBS sample 
designation 


Pemperature Stress range 


Unsaturated 
polyester resin ; . thin 
4 1») 1”) 0 to 6,000 
E 150 200 0 to 7.500 
F 140 (estimated 200 
Gi 130 (estimated Less than 150 
lerylic addition 


0 to 4.000 


0 to 4,000 


reain 
Cc Much less than 150 
Silicone resin 
B 5 325 (for 65%, esti- | 0 to 7,5 
mated 


Much less than 150__| 0 to 5,000 


Phenolic resin 
H 220 (estimated 0 to 5,000 
Melamine resin 


J 25 220 (estimated 0 to 10,000 


* Tests at 25° C after 200 hr of heating at temperature 7, and tests at tem 
perature 7 after 0.5 or 200 hr of conditioning at temperature T 


than the temperatures for the flexural strength 
It must be emphasized in addition, that for the 
modulus of elasticity data, temperatures different 
from those given in table 6 would apply for other 
stress ranges; when the stress-strain diagram is 
nearly linear for the ranges considered, the effect 
of differing stress range is slight. 

It is of interest to compare the values in table 6 
with data obtained by Sieffert and Schoenborn [10]. 
The latter authors determined the weight loss as a 
function of time for a number of laminates at 
various temperatures; from these data curves were 
plotted of weight loss versus temperature. It was 
found that for each material the weight loss as a 
function of temperature began to increase rapidly 
beyond a certain temperature. By using an 
arbitrary criterion, a “critical thermal instability 
temperature’ was determined from the plotted 
The 


ported for several laminates are shown in table 7. 


data. values Sieffert and Schoenborn re- 

It is evident that there is close agreement be- 
tween the critical thermal instability temperature 
and the temperatures given in table 5 for retention 
of 30 percent of original flexural strength. 


5. Effect of Speed of Heated Air on Changes in 
Flexural Properties 


The data in tables 2 and 3 for tests at 25° C 
after heating the specimens at 200° C either on 


the oven shelf or in open bottles were examined 
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TaBLe 7. Critical thermal instability temperature and resin 
loss data reported by Sieffert and Schoenborn [10] for glass- 
fabric laminates 


Critical | 
thermal Temper- 
Resin instabil- ature for | 
Glass-fabric laminate con- ity tem- 5° resin | 
tent perature loss (72 | 
M=4,72-hr hr) 
data 
Percent Cc ¢ 
Polyester 31 10 140 
Silicone A 52 350 375 
Silicone B 40 355 250 
Low-pressure phenolic 2s 210 2) 
Melamine -formaldehyde 
A 40 235 115 
Melamine -formaldeh yde 
B 45 240 100 


statistically (see section V). It was found that 
both the flexural strengths and flexural moduli of 
elasticity for the polyester laminates and the 
phenolic laminate were significantly lower for the 
case of heating on the oven shelf, i. e., in rapidly 
moving air, as compared to heating in the open 
bottles. However, for the melamine laminate J 
the differences were not significant. 

These results indicate that information on the 
effect of prolonged heating on the glass-fabric 
laminates under different atmospheric conditions 
is needed. The degradation of the resin may 
depend on the rate of removal of degradation 
products as well as the rate at which heated air 
containing oxygen, an important agent for some 
degradation reactions, is swept over the surface of 
the material. It follows that the use of data of 
the type presented in this report on the properties 
of laminates after prolonged heating should be 
used for design purposes with caution, since in 
practice the conditions of heating may differ 
widely from those that apply to the data given 
here. Likewise, interlaboratory comparisons of a 
given laminate may give discordant results if the 
heating conditions are not duplicated. 


6. Loss of Weight on Heating 


After 200 hr of heating at only 150° C, all the 
samples except the silicone exhibited losses in resin 





of 2 to 8 percent of the initial amount, with the 
melamine laminate showing the greatest loss. As 
the temperature of heating was raised, the loss in 
amount of resin increased very rapidly for all the 
polyester laminates and the one phenolic laminate 
For these samples the resin loss is roughly four- 
fold at 200° C and twelve- to thirty-fold at 250° C, 
compared to the losses at 150° C. In contrast to 
this behavior the melamine laminate showed much 
less rapid increases in resin loss (see table 4). 


V. Statistical Analysis 


l. Precision of Results 


The standard deviation was determined for the 
data in tables 2 and 3 for tests at various con- 
ditions. For a given sample it was observed that 
the standard deviation tended to be proportional! 
to the mean value. For this reason, further 
analysis was made with the coefficient of variation 
instead of the standard deviation. The coefficient 
of variation is defined as 100 times the ratio of the 
standard deviation to the mean value. 

The coefficient of variation was not reported for 
each flexural strength value in table 2 or each 
modulus of elasticity value in table 3 because such 
statistics, based on only five observations, are 
subject to wide variability. It was decided to 
analyze the individual coefficient of variation with 
respect to its variation from sample to sample and 
between the different temperatures. This was 
done to determine in what manner these data 
could best be combined for more reliable estimates 
of precision. The analyses are discussed in the 
following paragraphs. 


(a) Flexural Strength 


The coefficient of variation (C,) was determined 
for all the values in table 2. A statistical analysis 
of some of these data indicates some variation of 
the value of C, between samples. No consistent 
trend for the effect of conditioning or test tem- 
perature on the value of C, for the samples was 


noted. Therefore, the estimates of C, within 
each sample were pooled with the following 
results: 
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Laminate Cc. 


Unsaturated polyester: Percent 
d 6. 6 
E 5. 1 
, 7. 3 
G 8.4 
Acrylic addition 
Cc 8.8 
Silicone 
B 5. 2 
Phenolic 
H 10 
Melamine 
J 2. 9 


In a few instances, usually for the highest testing 
or conditioning temperature and when the flexural 
strength was a small fraction of the initial value, 
a very large value of C, was obtained; such values 
were omitted in obtaining the above estimates 
of C, 

(b) Flexural Modulus of Elasticity 


The coefficient of variation was determined for 
table 3. A_ statistical 
analysis was carried out for the tests made either 
at 25° C or at a temperature 7 after 200 hr of 
heating at 7. It was found that (1) the effect of 
conditioning temperature on C, varies significantly 


the modulus values in 


from sample to sample; (2) C, cannot be assumed 
to be the same for all samples; and (3) the effect 
of testing temperature, i. e., 25° C compared to 
T, on C, varies significantly from sample to 
sample. The above effects do not seem to be 
important except the sample-to-sample variability. 
It was therefore decided to combine the estimates 
of C, for the 


values of C, for the modulus of elasticity are: 


various samples. The estimated 


Laminate Cc. 


Unsaturated polyester Percent 
4 i. 0 
E 5.1 
F 8.8 
G 14 
Acrylie addition 
C 7 
Silicone 
B 1.5 
Phenolic 
H 10 
Melamine 
J 3. 0 


Strength of Laminates up to 300° C. 


In a few instances, usually for the highest testing 
or conditioning temperature and when the mod- 
ulus of elasticity was a small fraction of the 
initial value, a very large value of C, was ob- 
tained; such values were omitted in obtaining the 


above estimates of C,. 


2. Thermal Shock in Tests After 0.5 Hour of 
Heating 


A statistical analysis was made of the flexural 
data for the tests after 0.5 hr of heating in which 
duplicate sets of specimens were tested, one set 
with and one set without asbestos paper on the 
metal specimen tray. All the samples do not 
behave alike with respect to differences in flexural 
strength or modulus of elasticity for the two tray 
conditions. Although in some cases there appears 
to be a significant effect, the data as a whole do 
not support the assumption that the thermal 


shock is serious. 


3. Effect of Heating in Still Air Compared to 
Heating in Rapidly Moving Air 


A statistical analysis was made of the flexural 
data for the tests at 25° C after 200 hr of heating 
at 200° C with the specimens exposed (1) in open 
bottles, and (2) on the oven shelf. It was found 
that flexural strength and flexural modulus of 
elasticity were similarly affected for specimens 
contained in open bottles compared to those 
placed directly on the oven shelf during the heating 
The data for each sample were tested 
separately. Samples A, G, and fH all 
significantly greater losses in flexural strength and 
modulus of elasticity for shelf heating compared 
to bottle heating. For sample EF the effect was 
significant for flexural strength but not for modulus 
For sample J no significant effect 


period. 
show 


of elasticity. 
was noted, 


VI. Conclusions 


The silicone laminate is superior to the other 
laminates tested in retention of flexural properties 
at temperatures of 250° to 300° C. When tested 
at about 325° C after either 0.5 or 200 hr of 
heating at that temperature, this laminate retained 
at least 30 percent of its initial flexural strength 
and over 50 percent of its initial flexural modulus of 
elasticity. 
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The phenolic laminate showed good retention of 
its flexural properties when tested at elevated 
temperatures after heating for 0.5 hr at the test 
temperature. When tested in this way at 375° 
C, the flexural strength was 15,000 Ib/in.*, and 
the modulus of elasticity was 1.010* Ib/in., 
values roughly one-third of the corresponding 
values for the controls. The behavior of the 
phenolic laminate after prolonged heating was 
much less satisfactory; at 25° or 250° C after 
200 hr of heating at 250° C the flexural strength 
was much less than 1,000 Ib/in.’. 

The melamine laminate, although much inferior 
to the phenolic laminate when tested at a tem- 
perature of 375° C after 0.5 hr of heating, was 
superior to the phenolic in retention of flexural 
strength after 200 hr of heating at 250° C. 

The four polyester laminates lost at least 80 
percent of their flexural strengths when tested at 
250° C after either 0.5 or 200 hr of heating. 

The areylic addition laminate lost at least 80 
percent of its flexural strength when tested at 
150° C after either 0.5 or 200 hr of heating. 

For each laminate and each test condition, loss 
in flexural modulus of elasticity correlates with 
loss in flexural strength; the former loss is less 
than the latter. 

The loss in flexural properties during prolonged 
heating may depend on the method of exposing 
the specimens in the oven. For specimens of the 
polyester and phenolic laminates, heating in a 
circulating-air oven in open bottles caused signifi- 
cantly less loss in flexural properties than did 
heating on the shelves of the oven. 

These 
laminates, are based on data for only one sample 
of each type of glass-fabric laminate and hence 
tentatively as being 


conclusions, except for the polyester 


can be considered only 


representative of the various types of material. 
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Attenuation of Oscillatory Pressures in Instrument Lines’ 
By Arthur S. Iberall 


A theoretical investigation has been made of the attenuation and lag of an oscillatory 
pressure variation applied to one end of a tube, when the other end is connected to a pressure- 


sensitive element. 


An elementary theory based on incompressible viscous-fluid flow is first developed. 
The elementary solution is then modified to take into account compressibility; finite pressure 
amplitudes; appreciable fluid acceleration; and finite length of tubing (end effects). Account 


is taken of heat transfer into the tube. 


The complete theory is derived in an appendix. The results are summarized in eight 


graphs in a form convenient for use in computing the lag and attenuation of a sinusoidal 


oscillation in a transmission tube. 


l. Introduction 


In many industrial processes, it is necessary to 
know or to utilize the pressure at one or more 
points in a fluid conduit. It is not always possible 
to connect an instrument directly into the conduit 
Instead, recourse must be had to 
In the case that a 


at those points. 
remote indication or control. 
fluid is used for transmitting the pressure, it is 
often of interest to the designer or user of such 
systems to know their response to variations in 
pressure. At the present time, the only solution 
easily available to the engineer is generally based 
on an elementary theory that considers the system 
as equivalent to an R-C electrical network. (See, 
for example, NACA Technical Note 593, Pressure 
drop in tubing in aircraft instrument installations, 
by W. A. Wildhack.) The main defect of the 
theory is that it does not provide criteria for the 
limits of its applicability. 

In the present paper, a relatively complete treat- 
ment is given for the transmission of oscillatory 
pressures in tubing. Primary consideration is 
given to simplifying the design of high-quality 
transmission systems for relatively low frequencies. 

The elementary solution is derived and then 
extended to apply for oscillatory pressures that are 
an appreciable fraction of the absolute mean 
pressure, for appreciable frequencies of oscillation, 
and for tubing short enough to require end cor- 


This work was supported by the Office of Naval Research under a project 
‘Basic Instrumentation for Scientific Research.” 


Attenuation of Pressure in Tubes 


rections. The effect of heat transfer in modifying 
the oscillatory response of the tube is also dis- 
cussed. 

The chief utility of knowing these corrections is 
that it permits the designer to choose the size of 
tubing for specific applications with greater con- 
fidence than can otherwise be done. 

In the next section, the elementary theory of 
transmission lags is developed, and the corrections 
are discussed. The complete theory is presented 
in graphical form for the convenience of the de- 
A number of examples of the use of the 


signer. 
This section is then 


design charts are also given. 
followed by a mathematical appendix in which 
the more exact results are derived. All math- 
ematical symbols used in this paper are defined in 
section II and also when they are first used. 


II. List of Mathematical Symbols 


A=tube area. 
velocity of sound. 
D=inside diameter of tube. 
E-=elastic modulus of tube. 
F=correction functions. 
K=thermal conductivity of fluid. 
L=tube length. 
M=mass flow. 
N=dimensionless parameter of fluid regime. 
()=volumetric flow. 
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R=a volume ratio. 
R,= Reynolds number. 
T= absolute temperature. 
V =instrument volume. 
b=compressibility factor for liquid. 
e=any arbitrary constant. 
c,»—specific heat of fluid. 
f=any arbitrary function. 
q=a Bessel function argument. 
h=a Bessel function argument. 
k=compressibility of a liquid. 
/=entrance length. 
m=exponent of “polytropic’ 


expansion in in- 
strument volume. 

n=exponent of ‘‘polytropic” expansion in tube. 

p= pressure. 

s=tube wall thickness. 

t=time. 

u—=axial velocity. 

r=axial distance along tube. 

y=a dimensionless axial distance variable. 

dimensionless parameter of fluid regime. 


‘ 


y=ratio of specific heats. 
5= phase angle. 

n= density ratio. 

\= time constant. 
u=fluid viscosity. 


v= kinematic viscosity. 


= fractional pressure excess. 
p= fluid density. 

o= Prandtl number. 

¢= velocity potential. 

x =attenuation factor. 
¥—attenuation parameter. 
w=angular frequency. 


Ill. Elementary Theory 


Figure 1 is a schematic drawing of the system 
that will be discussed throughout the paper. A 
tube transmits fluid pressure from a conduit to 
a pressure-sensitive instrument. The conduit 
applies an oscillatory (sinusoidal) pressure to the 
entrance of the transmission tube. The tube, 
which transmits the pressure, is characterized by 
a constant cross-sectional area and its length. 
The pressure-sensitive instrument, which receives 
the pressure, is characterized by its enclosed 
volume. It is assumed that if the walls enclosing 
the instrument volume are flexible (either elastic 
or piston-like), the enclosed volume can be re- 
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placed by a larger equivalent rigid volume th; 
will store the same mass of fluid per unit pressu; 
change. It is further assumed that the pressur: 

sensitive instrument will be so chosen that i 

indication is independent of the frequency « 

expected pressure oscillations. 

In deriving the elementary theory, it is assum: 
that Poiseuille’s law of viscous resistance holds ; 
each point in the tube; that the fluid is incom 
pressible in the tube; that the sinusoidal pressu: 
oscillations at the beginning of the tube are of 
small amplitude compared to the mean absolute 
pressure; and that, if the fluid is a gas, it expands 
and contracts isothermally in the instrument 
volume. 


aa 




















Pp —¢ - a 
_— 
! 2 3 


Figure 1. 
system (1-conduit, 2-transmission tube, 3-pressure instrv 


Schematic diagram of a fluid transmission 


ment). 
p=potAp cos wt 


The same assumptions applied to an incom- 
pressible fluid (e. g., a liquid) lead to the conclusion 
that there is no loss in amplitude or lag in a liquid- 
filled system as a liquid would not expand o 
contract in the instrument volume. 

We may write 


Op _ 128 po 
"= Bp ( 


for Poiseuille’s law, and 


oM aii Op 19 
Or ~ Ot ; 
for the equation of continuity. Here 
p=instantaneous pressure at any point in 
the tube 
z=distance along the tube measured from 
its entrance 
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up mean fluid viscosity 
D=tube diameter 
(}=volumetric flow at any point in the tube 
p=instantaneous density at any point in the 
tube 
\f/=mass flow at any point in the tube 
A=cross-sectional area of the tube 
t=time. 
We infer from the equation of continuity and 
e assumption that the fluid is incompressible in 
tube (i. e., 0p/Ot=0) that the mass flow, and 
erefore the volumetric flow, does not vary along 
he tube, but at most varies only with time (the 
iid motion is piston-like 
sv differentiating eq 1, we then obtain 


O-p 
Or" 


0 3) 


long the tube 
Qur boundary conditions are that at 7=0 


Pp p p Ape ade } 


sinusoidal pressure variation about the mean 


pressure, and that at r= L 


V op 
i) 
« p or 
Op 128 po 
Ou s DS @ 


The first line of eq 5 expresses the rate at which 
a compressible fluid entering a rigid volume builds 
ip pressure, whereas the second line of eq 5 states 
that the flow into the volume is limited by the 
pressure gradient at the end of the tube. Here 
Po=Mean pressure at the entrance 
Ap=amplitude of the pressure oscillation at 
the conduit 
V=instrument volume 
w angular frequency of the pressure oscil- 
lation 
L=length of the tube. 


It is convenient to introduce a new variable &, 


the fractional pressure excess, defined as 


Pp Po. 6) 
Po 


¢ 
s 
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so that eq 3, 4, and 5 become, respectively, 


Ore i 
a 0 ri 
Or 
at r=0 
é Ene sS) 


and at z=L, 


Ou r po D' ot 
9 
Ay OF 
L ot 
where 
128 wy a 
Xo r py D | 
10 
an (LY V 
32) AL 
Here 


¢ —fractional pressure excess 
amplitude of the fractional pressure ex- 
Ap Po) 


ft 


cess at the origin 
Ay=a time constant of the system. 

It is of further convenience to separate the 
pressure excess into a part that varies with + and 
one that varies with ¢. 

Let 


E— fer! 11 
where & is the maximum amplitude of the pressure 


excess at any point of the tube. 
Our equations then become 


i ke t 
> =( 12 
dz’ 
ats 0), 
E=§ 13 
and at z=L 
Fs iia 
“* ee ye (14 
dz L 7 
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The solution of eq 12, which satisfies eq 13 and IV. Discussion of Corrections 


14, is 
The assumptions made in the elementary theory 





. Anwo( - )i are restrictive, and in the appendix we shall mod {y 
E=fo 1+ row) (15) them, one at a time, until finally we arrive at r 
, complete solution that accurately takes into ac- ~ 
The ratio of the amplitude of the pressure excess count all first-order phenomena, and partial!) . 
at the end of the tube to that at the beginning of takes into account second-order phenome: . 
the tube is then given by Complete results are presented in convenient q 
graphical form in figures 2 to 9. | 
fi l 
E) 1 + guy 
(16) 
I 
1+ XoJ 








Or . oo 
where Pi 
LATTE XE NT 
Pm NNT tl 
,=maximum amplitude of the pressure \ . 

excess at the instrument volume N \ 
Xo =an attenuation factor. . N LIN 


| pease 


Here 








ote: 
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The real part of eq 16 is the attenuation in 8 WJ 
amplitude of the pressure excess, whereas the Ficure 2. Amplitude ratio of the fundamental \Eox/toly iv ad 

imaginary part is the phase lag, or a volume terminated tube as a function of a paramet: 


»roportional to frequency (x7) for various ratios of inst 
prop Jreq Y (X10) | 


ment volume to tube volume (x1m/xro) with large damping 



































7 I ) (zS1). 
& {1 T xi]! . (18) ' 
tan do= Xo, 
where 4, is the lagging phase angle. Fe 
We will regard eq 18 as the elementary solution 
et : ome 100 Te” 
of our problem. It indicates that a transmission | / / 
system is characterized by a time constant », 4, i 
which can be computed from a knowledge of the ” sie 49 9 31 / pp 
dimensions of the tube, the internal volume of & ¢ To 
the end device, and the average conditions of the ry V4 Kh - 
gas in the tube; and an attenuation factor xo, for -S 40 4 aoa 
° oO | | 
each angular frequency, from which one can © re Wl Wi Y . 
compute the attenuation and phase lag ina tube. 2 vA we A sant! _2 
The tube dimensions and the instrument volume pz ws 
. ; . - = 
furnish the analog to the resistance and capacitance be Ol ) , 
of an electrical network. Xt J 





In principle, although difficult in practice, from , 4 
. Ficgure 3. Phase lag of the fundamental (50)9 in a volur 
a knowledge of the response to a sine wave, one 

_— terminated tube as a function of a parameter proportion 
can obtain the response to square waves, step to frequency (xr) for various ratios of instrument volur 


function, etc., by Fourier analysis. to tube volume (x10/x1) with large damping (z<1). 
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Ficure 8. Relative amplitude of the double harmonic dis- 


tortion (\E:1/Eole = in a volume terminated tube as a func- 





tion of a parameter proportional to frequency x7o) Jor 


various ratios of instrument volume to tube volume (x1/Xx ro 














t with large damping (z1) 
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Ficure 9. Phase lead (as measured on the fundamental time 
scale) of the double harmonic distortion (6;)5 in a volume 


terminated tube as a function of a parameter proportional 
to freque nmcey \XT0 for various ratios of instrument volume 


fa tube volume (x10/x79) with large damping (zS1 


The factors that must be taken into account 
are: 

1. Compressible flow in the tube. The effect of 
fluid compressibility is to introduce a time con- 
stant and corresponding attenuation factor (Ar, xr) 
depending on the tube volume in addition to the 
ones depending on the instrument volume. (The 
time constant and attenuation factor depending 
on the instrument volume will be referred to as 
h, and x, henceforth instead of Ay and x»). In 
terms of the electrical analog, the tube volume 
represents a distributed capacitance in addition 
to the equivalent capacitance of the instrument 
volume. 

2. Finite pressure excess. The effect of the ap- 
plication of a finite pressure excess to a compress- 
ible fluid in a transmission tube is to introduce 
harmonic distortion and to modify the mean pres- 
sure. However, the attenuation of the funda- 
mental is essentially independent of the magnitude 
of the pressure excess. The percentage of distor- 
tion is approximately proportional to the applied 
pressure eXcess, 

3. Fluid acceleration. The effect of fluid inertia 
is to modify the time constants of the system. 
Both the attenuation of the fundamental and the 
magnitude of harmonic distortion are affected. A 
dimensionless parameter z analagous to the “Q” 
of an electrical system characterizes the fluid 
regime and determines whether fluid inertia may 
or may not be neglected. 

When fluid inertia is negligible, a transmission 
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tube acts like a highly damped system; when flui: 
inertia is large a transmission tube acts like ar 
undamped system, and elementary acoustic theory 
is applicable. 

4. Finite length of tubing. The effect of flui 
acceleration at the ends of the tube results i 
further distortion of wave form, which must by 
taken into account in short tubes. 

5. Heat conduction. If there were no heat trans 
fer from outside the tube to inside, the oscillatory 
processes would take place adiabatically; if ther 
were perfect heat transfer into and through th 
tube, the processes would take place isothermally 
The effect of finite heat conduction is to make th: 
real process occur in between these extremes 
although in a rather complicated fashion. At low 
frequencies the process may be regarded as 
isothermal. 

Although an exact result is given in the appen- 
dix, it is advantageous to utilize the thermody- 
namic equation of condition, discussed in_ the 
following section, for elucidating the problem ol 
attenuation in tubing. 


V. Thermodynamic Equation of Condition 


In the case of an oscillatory variation of fluid 
flow, the equation relating the thermodynamic 
parameters of the fluid lie between the adiabatic 
and the isothermal equations of condition. For 
high frequencies, as in sound waves, it is well 
known that the adiabatic equation holds. How- 
ever, for viscously damped motion, the adiabatic 
relation is not, in general, attained. 

For a gas, we assume and justify in the appendix 
the processes can be described as “polytropic”’, 
that is, characterized by a constant exponent n, 
in the expression 


=p" 
/ p ae) 
1+§=7" ) 
with 
1Sjn| S7, 
where 


n=exponent of the “polytropic’”’ expansion in 
the tube 

y=ratio of specific heats 

n=density ratio (p/ po) 

po average density in the tube. 

¢ is used to indicate any constant. 

The viscosity of gases is independent of the 

pressure, and, as an approximation, proportional to 
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e absolute temperature. (The more rigorous 
proximation is that the viscosity is proportional 
ae 

[TP /01 


approximated by the temperature to a power 


-¢/7T] but over a small range this can 
osetoone. For example, for air at room tempera- 
re, a power of 0.8 fits experimental data quite 
ell. The difference from unity is unimportant 
1 our purpose 


Therefore, 


v4 cT, 7 
n—l 
iv 
l 7 &) ms 
u > (20) 
Me ’ 
ui , 
Mo J 





follows from the gas laws and eq 19. Here y is the 
nstantaneous fluid viscosity. and 7'is the absolute 
temperature 

Equations 19 and 20 thus express the variation 
of viscosity, density, and pressure in a polytropic 
process in a gas. At low frequencies, the poly- 
tropic exponent may be taken as equal to unity. 

For livuids, we assume that the equation of 


condition in a polytropic process is given by 


p Pot cp", (21) 


where 


For liquids, however, ¥ lies so close to unity 
that we may satisfactorily assume n=] 
Equation 21 can then be written in the form 
99) 


n=1-+4 be, (22 


where b=a compressibility factor | ko po) 
ky=liquid compressibility at average condi- 
tions in the tube 
The variation in viscosity of a liquid over a small 
range of temperature cun be neglected, so that in a 
poly tropic process 


= Mo. (23) 


Actually the implication in eq 22 and 23 is that 
n a liquid-filled transmission line, the effect of 
conditions appreciably different from isothermal 
s negligible. 
It is also necessary to take into account heat 
exchange at the pressure element. 
For an isothermal process with a gas in the 
nstrument volume, we previously assumed that 
V Op 


@ Po Or (5) 
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represents the influx of fluid. If, instead, a 
polytropic process in the instrument is assumed, 
characterized by an exponent, m (the heat ex- 
change may differ in the tube and instrument 
volume so that m is not necessarily equal to n 
then eq. 5 should be modified to 


Vo Op 
mp ot 
24) 
Nn V Op 
} 
& m p Ot } 
in the case of vases: or to 
V7 
Qu SF, (25 
Po Ot 


for liquids. 

If the fluid is regarded as a spring, the exponent 
of the polytropic process for a gas, or the com- 
pressibility of a liquid may be viewed as quanti- 
ties that make the fluid spring stiffer in the case 
of gases, or almost infinitely stiff in the case of 
liquids. It is shown in the appendix that these 
polytropic exponents modify the time constants of 
the tube and volume. 


VI. General Procedure, with Examples, 
for Computing Transmitted Pressure 


The computation of the attenuation and phase 
lag at one end of a transmission tube of a sinusoidal 
pressure variation imposed at the other end can be 
These 
figures are based upon the theory largely developed 
in the appendix 
primarily for the attenuation at the fundamental 


carried out with the aid of figures 2 to 9, 
The computations are made 
frequency. An estimate of the distortion arising 
from finite input amplitudes with high damping 
is made in the appendix. The computation for 
the first harmonic in the distorted output can be 
made with the aid of figures 8 and 9. An outline 
of procedure for making computations follows 
1. Compute 


(26) 


a dimensionless parameter of the fluid regime that 
characterizes the amount of damping present. 
When this parameter is less than | (large damp- 
ing), use figures 2 and 3; when greater than 100 


(small damping), use figures 6 and 7. For inter- 
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mediate values of this parameter, use figures 4 and 
5 as an aid to interpolation. 
2. Compute the attenuation factors 


32u0 (L\? 
Xr ( ) ? 
0 Po 


Xe 1/V 
= m (AL) 


D 


for a gas, or 


Xr, =32yokow ( ‘ ; +hh EE - ), 


Xi V l 


Xm (14 — LE *) 


for a liquid. These quantities, xyz and xo, are 
factors based on the tube volume and instrument 
volume, respectively. The zero subscript means 
that they are values for the case of large damping. 


3. Compute the input pressure excess 
_ Ap 
Po 


4a. For values of z less than 1, 


(31) 


enter figure 2 
with xr and xz/x7r to find the amplitude ratio 


Eor/t|9 and enter figure 3 to find the lagging phase 
angle (do). 
4b. The output pressure excess is then com- 


puted from 
(32) 


5a. For values of z greater than 100, compute 


wh _[ 2xr0 ph, (33) 
c Sy iat 


_# r( an ): (34) 
Xr XTo 
For liquids, assume Y= 1. 
The quantities x, and x, are the attenuation 
With low 
damping, it is convenient to use the dimensionless 
parameter wL/C, which is proportional to fre- 
quency, as the independent variable. 

5b. Enter figure 6 with wL/C and x,/x,r to find 
the amplitude ratio |fo,/&! and enter figure 7 to 


factors for the case of low damping. 


find the lagging phase angle 49. It is necessary 
to estimate the phase angle by interpolation. For 
very small values of yr) compared to 1, the lagging 
phase angle is zero up to the first resonance. In 


figure 7, curves have been presented to indicate 
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the phase angle for [yx7o/16]'?=0 and [yx 70/16)!” 
One may linearly interpolate between these curvy: 
on the basis of [yx7o/16]'? for values lying betwe: 
0 and 1. 

6a. For values of z lying between 1 and 10) 
one may interpolate between the values of amp! 
tude ratio and lagging phase angle obtained | 
step 4a and those obtained in step 5b by the us 
of figures 4 and 5. Enter figures 4 and 5 with , 
and x79/xro to find the amplitude ratio Sos & | and 
lagging phase angle 49. These are the values fo 
In order to interpolate, plot a logarithmi: 
graph with 2 as abscissa and the amplitude rati: 
Plot the values from 
1, from step 6a at 2=6.25, and fron 

100, draw a curve through thes: 


2=6.25. 


or phase angle as ordinate. 
step 4a at 2 
step 5b at 
three points, and interpolate on this curve for th 
intermediate value of 2. 


1. Computation of Double Frequency Distortion 


1. This computation represents only an esti- 
mate of the double frequency distortion and _ is 
strictly valid only for values of 2 less than | 

Compute Xr and xX,/Xrm. Enter figures 8 and 
9 to obtain the relative amplitude ratio | £),/E\9/ 
and leading phase angle (4,)) for the double fre- 
quency wave. The leading phase angle is meas- 
ured on the time scale of the fundamental, where 
both the fundamental and double frequency waves 
are cosine terms. 

2. Compute the pressure excess of the double 


frequency £,.9 from 


: l 

t 

SIL 0 (; 
s 


The various quantities in the above section are 


fin ) ¢2 
Fag $o0- 
s 


‘0 0 0 


defined below: 
uo mean fluid viscosity. 
vy==mean kinematic viscosity. 
mean fluid pressure. 
ambient pressure external to the tube. 
amplitude of the applied sinusoidal pressure 
applied fractional pressure excess. 
pressure excess of the fundamental at the 
instrument volume. 
pressure excess of the double frequency at 
the instrument volume. 
lagging phase angle of the fundamental at 
the instrument volume. 
leading phase angle of the double frequency 
at the instrument volume. 
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mean fluid compressibility. 
\’=equivalent rigid internal volume of the 
instrument. 
A=internal cross section of the tube. 
= internal diameter of the tube. 
s=wall thickness of the tube (assumed small 
compared to the diameter). 
E=elastic modulus of the tube material. 
length of the tube. 
(’=velocity of sound in the fluid. 
y=ratio of specific heats of the fluid (assumed 
to be one for liquids). 
m=coefficient of the polytropic process in the 
instrument 
information, it may be assumed to be 
ope.) 
dimensionless parameter characterizing the 
fluid regime. 
w=angular frequency applied. 
yrp=attenuation factor based on the tube 
volume. 


attenuation factor based on the instrument 


a 


I 
volume. 

Subscript 7’ refers to parameters based on tube 

refers to 


volume; subscript J parameters 


based on end volume; O or 1 following a 
T or J denotes the fundamental or first 
harmonic; an end subscript of 0 denotes a 
value for the case of large damping. 

The attenuation of the fundamental may be 
validly computed from the formulas developed in 
this paper when 

y% 
CD 
(36) 
VyW 


(? 


The second harmonic distortion, which was only 
stimated approximately, may be validly com- 
puted from the formulas developed when 


wl)? 
4, 


l, (37) 


ind when the applied pressure amplitude is suffi- 
iently small at the applied frequency to permit 
aminar flow. 3. Compute the mean pressure in 
he instrument volume, which is larger than the 


lean pressure at the tube entrance by 


E,Ap[1 —| Eox/Eol3]/4 
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volume. (In lieu of other 


2. Examples of Computations 


The calculation of attenuation by the general 
procedure outlined above will be illustrated by a 
number of examples. 

(a) What is the longest length of *{.-in.-inside- 
diameter tubing that can be used to transmit air 
pressure to a Bourdon pressure gage (equivalent 
internal volume assumed negligible) up to a fre- 
s with a loss in amplitude not 
What will be the double 
For air assume py=2* 10™4 


quency of % ¢ 
greater than 25 percent? 
frequency distortion? 
poise, w—1/6 stokes, m=1, y=1.4, po=10° 
dynes/em? (atmospheric pressure), angular fre- 
quency @=f. 

Using eq 26, z=1.1 (computed in consistent 
units). This value is sufficiently close to unity to 
permit the use of figures 2 and 3. Enter figure 2 
with koe £olo=0.75 and xr0/xrm=0, since the in- 
strument volume is negligible, to find yx 7)=2.1. 
Compute Z in eq 27 to be 160 feet. 

2.1, to find that the 
maximum phase lag will be 53 degrees. 

Entering figure 8 to find that the relative ampli- 
tude of the double frequency as £5 9/Eo = 0.35. 


For initial pressure excesses of 0.1, 0.3, and 1, 


Entering figure 3 with x7» 


respectively, the double frequency amplitude, rela- 
tive to the input amplitude, will be 3', 10%, and 
35 percent, while the mean pressure will increase 
0.0010, 0.010, and 0.11 of an atmosphere, respec- 
tively. 

(b) What lengths of 0.1-in.-inside-diameter tub- 
ing (nominally *{,-in.-outside-diameter tubing) can 
be used for quality transmission of air pressure for 
frequencies up to 1, 10, 100, 1,000 c/s into pressure 
instruments with equivalent rigid volumes of 0.1 
and 1 in.*? 

We will define quality transmission as that in 
+ 5-percent change 


+ 30° 


which there is no more than 
in fundamental amplitude or more than 
phase shift (whichever is more stringent). 
Assume that po=2107* poise, vy=1/6 stokes, 
1.4, D=0.1 in., A=0.0079 in.*, po=10° 


0.0012 g/cm’. 


m=1, 4 
dynes em’, Po 

We will calculate for each frequency separately. 
b (1). f=le/s: 


Using eq 26, ) 


0.61; therefore, use figures 2 
and 3. 
Assume AL 
Enter figure 2 for |&z/t\o=0.95 to find xn 
0.80. 


©, therefore, by eq 28, xs0/xro= 9. 
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Enter figure 3 for (é9))>=30° to find x~w=1.1; 
use 0.80 since it is more stringent. 

Calculate L from eq 27 to be 450 in. 

Calculate AL to be 3.5 in. 

Compute x7o/x7ro from eq 28 to be 0.029 forV=0.1 
in.?; 0.29 for V=1 in 

On figure 2, xr is modified negligibly for V 
0.1 in. 

Therefore, L 

Reenter figure 2 
0.29 to find xro=0.5 

Calculate L to be 350 in. 
b (2). f=1,000 e/s: 

Using eq 26, 2=610; therefore, use figures 6 


37 ft for V=0.1 in’ 
0.95, 


450 in. 
for |£or/Eolo X10/Xro= 


29 ft for V=1 in. 


and 7. 

Assume x;/xr=49 (the line volume will probably 
be small). 

Enter figure 6 for | ,/&|—1.05 to find wL/C 
0.031. 

Calculate L from eq 33, 26, and 27 to be 0.066 in. 

Calculate AL to be 0.00052 in’ 

Using eq 34, it is seen that x s9/x 79 is greater than 
assumed, so that wl/C, and therefore L, is less 
One may note that 


so small that the 


than the previous estimate. 
the estimated length will be 
theory essentially predicts that no transmission 
tubing at all may be used. In fact, the acoustic 
impedance of the entrance orifice into the pressure 
instrument or the mechanical impedance of the 
pressure instrument itself will probably govern the 
response at this high frequency. 

b(3). f=10 e/s: 

Using eq 26, 6.1; therefore, use figures 4 
and 5. 

Assume AL= @, : 

Enter figure 4 with |£oz/Eo! 
find xro=0.12. 

Compute Z from eq 27 to be 58 in. 

Compute AL to be 0.45 in-’. 

Compute x/xro from eq 28 to be 0.22 for 
V=0.1 in.’, =2.2 for V=1 in? 

In figure 4, XTo is modified to about 0.07 for 

0.1 in.® 
Therefore, LZ is reduced to about 4 ft for V=0.1 


} 


therefore, xso/x7=0. 


1.05 and y=1.4 to 


in. 
Enter figure 4 for tou &|=1.05, and xz0/xro=2 
to find xro=0.018. 
Compute L to be 22 in. for V=1 in. 
Compute AL to be 0.17 in* 
Compute xz0/xro=6. 


Enter figure 4 to find xz 


8 


0.007. 
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Compute ZL to be 14 in. 
Compute AL to be 0.11 in.* 
Compute xr0/xro=9. 
Enter figure 4 to find x79 
Compute LZ to be 11 in. 
Compute AL to be 0.09 in.* 
Compute xyo/xr=11. 


0.004. 


In figure 4, x79 is modified negligibly. 
Therefore L=about 1 ft for V=1 in.* 
To check the phase angle, enter figure 5 \ 
x10o/Xro= 11, and x7—0.005, to find 4°. 
b(4). f=100 e/s. 
z=61 (interpolation is necessary). 
First estimate from figure 6 and 7. 
Assume x;/xr=9. 
Enter figure 6 to find wl/C=0.068 
As in b (2), compute Z to be 1.5 in. 
Compute AL to be 0.011 in. 
Compute x;/xr=9.1 for V=0.1 in.’; =91 for V=1 
in.* 
By figure 6, wLZ/C is negligibly modified for 
V=0.1 in; 
Therefore L 
estimate. 
For V=1 in., 
small tube is predicted, so that the impedance of 


1.5 in. for V=0.1 in is our first 


we find again that an extremely 


the entrance orifice will probably govern. 
For V=0.1 in’ and LZ 
to be .O1. 

From 


negligible. 


1.51n., estimate [yx 79/16] 


figure 7 we find that the phase lag is 


0.0009, from (8.2) for w=200r. 


9 to find |z/E| = 1.00 


Compute x79 
Enter figure 4 for xz0/x70 
Interpolating between |£oz/f|—1 at 
tor/&| = 1.05 at z=100 for z=—61, 
is negligibly affected. 

Therefore L=1.5 in. for V 


=6.25 and 


we find |£o; 


0.1 in 


VII. Appendix. Development of the Theory 
1. Introduction 


The difficulties of deriving, elucidating, and 
comprehending the mathematical results of trans- 
mission in tubing from a rigorous point of view 
have led the author to treat the problem in a series 
of somewhat artificial steps. Thus in the previ 
ous sections, the elementary solution was pre- 
sented, to give the reader a general view of th 
problem, even though many of the details of th 
solution were slurred over. Here steps are taken 
one at a time, to remove the restrictive assump- 
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ons made in deriving the elementary solution. 
evertheless, a complete solution to the problem 
not obtained. All first-order effects are treated 
the point 


equencies well into the sonic region. 


where the solution is correct to 
However, 
ily an elementary treatment is given for the 
eond-order distortion effects. It is felt that 
hen these second-order effects become appreci- 
ble, the solution presented is of no quantitative 
ility to the instrument system designer or user, 
suit is only indicative as to order of magnitude. 


2. Theory Corrected for Compressibility 
(Infinitesimal Oscillatory Pressures) 


In this approximation, the assumptions are 
Poiseuille’s law of viscous resistance; small frac- 
tional pressure excess; and that density, pressure, 
and viscosity are related by the equation of 
condition. 


For gases one can then write 


Op 128 pu 
Or x D @, 


Op _ 128 im 
p P1284, 


for Poiseuille’s law, and 


Ox colt 
for the equation of « ontinuity. 
One can eliminate the mass flow ./, to obtain 


32 Op 


D? dt (38) 


ae 


Or \ u Or 


By virtue of the assumption of small pressure 
excess, and the equations of condition (eq 19 and 
20), we can disregard the differentiation of p/m in 
eq. 38, and replace it by its mean value. Equa- 
tion 38 then becomes 


O*p 32u Op 
Or? = npol?® Ot’ 


O7— 32m Of 
npolF ot 


Utilizing the previous definition of X» (eq 10), 


eq 39 becomes 
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oe ( Ao \ 1 dE 
oz? \V nn \p ot 


OE Aro OE 
ox? L? ot 
where 
—) ro, (41) 
The significance of the new time constant Ago 
can be understood by inspection of the definition 
of A» (eq 10). One may note that Ago is a time 
constant based on the tube volume, AZ, instead 
of the instrument volume, V; and that it gives 
weight to the exponent of the polytropic process 
in the tube. It is thus related to the equivalent 
distributed electrical capacitance of the tube. 
The weighting by the exponent, n, arises from the 
fact that it represents the additional “stiffness” 
of the air column in the tube as a polytropic 


spring. 


If, as in the elementary solution, we separate 
our pressure variable into a space and time part 


é be ead” (1 1) 
eq 40 becomes 


dé ArqQW +: 


de L? J 


ad’ _XT0 2 
» DE 
dr IL? 
where 
XTo Arrow (43) 


The quantity x79 is an attenuation factor based 
on the tube volume. 

Equation 42 may be compared with the corre- 
sponding equation of the elementary solution, eq 12. 
It may be noted that it is necessary that x79 be 
small in order for the elementary solution to be 
valid. 

Referring now to eq 42, the boundary conditions 
are 
at r=0 


and atrz=L 


di 


dr 


(see eq 5, 9, and 24). 








We may redefine a time constant and attenua- 
tion factor for the instrument volume, which takes 
into account the polytropic process as 


bid No 
m (45) 
X70 hol 


At r=L, eq 44 therefore becomes 


dé Xi «> ; 
de L?* (46) 


The solution to eq 42, which satisfies boundary 
conditions (eq 13 and 46) is 


i 


r 
« P ¥TO , v 
g é ¥T0(Wro Wr0)¢ L-+-@"1 Vro T W10)€ we 
£0 eT Yo tro) +e °™(Wro—Yr0) 
(47 ) 
where 
XTo . - 
Vro I - ] (1+ 7) 
- . (48) 


( 
¥r0 ]X0 ) 


The new y's, which shall be referred to as 
attenuation parameters, are 

¥ro an attenuation parameter depending on 
tube volume; 

¥ro an attenuation parameter depending on the 
instrument volume. 

The ratio of the fractional pressure excess at 
the end of the tube & to that at the beginning 
of the tube & is then 


L ¥ro 


= : 9 
) Vro cosh Vro + Wr sinh Vro (49) 


fer fer 


It is instructive to examine the limiting values 
of this equation. For small pro, the attenuation 
approaches 


1 Seat ™ (50) 
& Ll+¥nm 1+ 7Xn 
the same result as in the elementary theory 
(see eq 16). 

For small values of Yo, the attenuation ap- 
proaches 


l 
: , (51) 
», cosh Vro 


~ 


c 
s 
¢ 
s 





which for small Yr). becomes 


frre fre 
& 
_ 
_ 


+- 
+ 


The form of eq 50 and 52 is similar. In fa 
for small values of both Wry and pp it is possib 
to define a composite attenuation factor x by t! 
relation 


XTo 
x Xo t+ Sr 
t ’ 
or | ] ) 
Aro 
r dot 6p 


such that the real magnitude of the attenuation 
is approximately 


[1+ x7] 


(18 


fer \otri 
| 


which preserves the form of the elementary 
solution. 

Equation 18 can be interpreted as meaning that 
the “‘proper”’ time constant of the system can be 
obtained by adding to the n weighted volume of 
the instrument, 1/[6]'* of the m weighted volume 
of the tube, and substituting this in the elementary 
formula for the time constant of the system. 

In principle, for larger values of Pro or Pm, a 
coupling coefficient (of approximately unity) could 
be introduced as an addition to the coefficient 
1/(6}'”, which would vary somewhat with the rela- 
tive magnitude of Yr and Yo, to permit strict 
preservation of the elementary form. It is, how- 
ever, simpler to compute attenuation from eq 49 

For liquids, we start from eq 38. 
ra) (? Op 32 Op. 2g 
or op D? dt 


As before, with the aid of eq 22 and 23, we 
obtain the result 


DE —32uob dE | 
Or? = pol ® Ot 





O7E Ary VE - 

oz" «—«L? Ot { » 
or 

dt - XrTo jé 

dx? L? 7% 
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vhere 


AL 
Aro ? bry 
i sai 
Do) 
XTo Arrow. 
The boundary condition at r=0 is 
g So (13) 
idatz ta 
dé (A bw) } é 
ile “ered T 
(56) 
di 3 
dx Xa] L 
where 
Axo hod ) 
(57) 
X70 yoo) 


see eq 5, 9, and 25). 


The form of eq 54 and 56 is identical with eq 42 
and 46, with the difference that the coefficient in 
the \’s is the very small compressibility factor 
rather than the reciprocal of the exponent of the 
polytropic process. Physically, this simply means 
that the liquid is a spring of almost infinite stiff- 
ness compared to the gas 

Because of the formal identity of the equations, 
the previous solution holds in toto, with the modi- 
fied value of Av. The following interpretation is 
now possible for the elementary result that there 
is no attenuation with liquids. The dX» time con- 
stant of elementary theory did not take into 
account the effect of liquid compressibility, which 
is small. If, however, \) is weighted by 5} (i. e., 
\=6bX») then the same attenuation curve holds 
for both liquids and gases, but with liquids we 
operate on the very beginning portion of the atten- 
uation curve for gases. 

There is one complication that should be con- 
sidered in liquid tube attenuation. Because of the 
small compressibility of liquids, it is often possible 
that the flexibility of the tube gives rise to a com- 
pressibility comparable to that of the liquid. The 
simplest way of taking into account the flexibility 
of the tube is to define and replace the compressi- 
bility factor of the liquid by an effective value 


band Ff. 
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, (58) 

Ep Po—PaD 

bakers 


where 
Po=mean liquid pressure; 
Po=ambient external pressure (usually atmos- 
pheric) ; 
E=elastic modulus of the tube material: 
s=wall thickness of the tube. 

In the derivation of eq 58, the assumption has 
been made that the thickness of the tube wall is 
small compared to the tube diameter. 

3. Theory Corrected for Finite Oscillatory 
Pressures 

In this section, we will determine the effect of 
finite fractional pressure excess on the attenua- 
tube. We assume only that the 
distribution holds. We will 
show that the effect of finite pressure excesses is 


tion in a 
Poiseuille velocity 


to excite higher harmonics, resulting in a distor- 
tion of wave form, and to raise the mean pressure 
along the tube. The higher harmonics are excited 
because of the nonlinearity of the equations. 

The method of solution selected will be that of 
expansion in harmonic series in which the excita- 
tion of sum frequencies only are considered and 
the difference frequencies are neglected, so that 
the solutions obtained are only valid for the lead- 
ing term of each harmonic. The second order 
term in the variation of the mean pressure will be 
estimated separately. We will assume open func- 
tions of the distance coordinate for the coefficients 
of each harmonic term of the series and show that 
the expansion is valid for moderate values of the 
initial pressure excess. It is obvious that these 
distance dependent coefficients must be the 
solutions of second-order differential equations 
in order to provide two sets of adjustable con- 
stants to satisfy the boundary conditions at the 
two ends of the tube. However, by considering 
the solution for an infinite tube (for which only 
one set of boundary conditions is required ) we 
shall be able to discuss the question of convergence 
of the solutions. 

For the purposes in view, it will turn out to be 
convenient to derive the equations on a density 
basis. Density and pressure are, of course, re- 
lated through the equation of condition. 
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For gases, we start from 
re) (2 Op 32 Op 
Or \ uw Or 1F ot 
By algebraic manipulation, in which eq 19 and 
20 are used to eliminate viscosity and pressure, we 
obtain 


(38) 


ra) 
or 


On 
” Or 


Aro On 
L? ot’ 


07’ = 2X TO On 
oz? —sCdL?:séOt’ 


a nonlinear partial differential equation. 
Assume as solution for the density ratio 
9= 1+ me" + 92e7** + age 
where 7; is the fractional density amplitude ratio 
of each harmonic (functions of z). 
At the moment assume that the applied pressure 


SJwt 1 
l 


, (60) 


wave has all the mathematical properties necessary 
to make the Fourier expansion of eq 60 valid. We 
will discuss this point again. 

When eq 60 is substituted in eq 59 and the 
coefficients of like terms in the respective har- 
monics equated, the following system of differential 
equations result for the coefficients 1,. 


TE) 03 de 


D gp? (mot mm+ NoN2) 
° ¥m\ 
. L/™ 


(mo + Non ) 


1 @? . 
“> dz? (3% mM TMM + nons) 


1 d@ ; 
5 dy? (1M +maumtes s+ 


i(*)n 





1-1 + MoM) 4 


The coefficient 4)(=1) has been added for com- 
pleteness. 

From eq 61 it can be shown that the coefficients 
n: for an infinite tube are equal to 


m- ime (+ te—11')y 


Ret a 


m1 


(2+ (4-2) "yw 


: Hi-2M€ 
(1X0+1X2[¢—2]'”) (21+ 


Nit-3 The 


~ (X0+1X 2t— y+ 


2X 2[i—2]') 


—(3+ [4 
(1X0+1X 2[1—3]'”) (2K 1+2X 2fr—3]'”) (8 K24+-3 XK 2fi—3 


The solution of these equations consists of 
complementary part that introduces two new co 
stants for each coefficient ;, and a particu! 
solution that depends upon the solutions for e- 
efficients with lower values of 7. The latter part 
represents the excitation of higher sum frequen: 
modes. 

Detailed investigation of the convergence of t! 
solutions for the various coefficients leads to the 
following conclusions: 

The differential equations of eq 61 will admi 
physically admissible and convergent solutions for 
any bounded periodic pressure or density wave 
at the origin because (1) either the wave at the 
origin has a derivative that is of limited variation, 
in which case the solution of eq 59 in series (eq 60 
is always valid; or (2) if it does not, spatial 
attenuation of viscous waves occurs so rapidly 
for higher harmonics that the wave will have a 
derivative of limited variation at a short distance 
beyond the origin, so that as far as effects down- 
stream are concerned, the given input wave can 
be replaced at the origin by a similar looking 
function (i. e., a finite polynomial instead of an 
infinite Fourier derivative of 
limited variation. As illustrations, we can replace 
a square wave by its first few harmonics, or a 
Weierstrass function by a smooth integrable 
function. In simpler language this means that 
in a viscous transmission tube, detailed or sharp 
wiggles in the initial disturbance (high harmonics 
are not transmitted. 


series) whose is 


The more practical question as to the rapidity 
of of the solution for the 
attenuation of a given entrance disturbance can 
be answered approximately by recourse to the 
solution for an infinite tube, for which only one 
set of boundary conditions must be satisfied. 


convergence series 


3) 2), 


Teo. 





|' 2) 
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The last line in eq 66 contains our desired con- 
y= zr (63) clusion. We can infer from it the maximum num- 
f TO e i ‘ 

L ber of terms that must be carried along in order to 
ee : ; know the distortion to within any desired ac- 
ere 7, is the constant of integration for the com- : a . 

‘ curacy. If we assume, for example, that we are 
ementary solution of each n;, and y is a dimen- 


: ; interested in only those harmonics whose content 
ynless distance variable. . 
Substituting the value of the coefficients from 
62 into eq 60, we obtain the result that, if at the 
igin (y=0), the density is written in the form 
> Ins! 
m 


a LF) 
4 243 —0.01. 
) 


4 


at the end of the tube is greater than 1 percent of 
the applied first harmonic, we can neglect all 
harmonics greater than the one for which 


e ¥ roll” 0.01, | 





For most practical problems, it can be shown 
that adequate information can be obtained from 





Sjwt 
ave Ww 5 )e T . . . 
+ 2X 2[1) y a knowledge of the first and second harmonic, and 
(64) rarely, the third harmonic. 
, or 1 To compute the harmonic distortion for a 
‘ i , Tr D oP " ‘ns , re Ww > ° . 
at any other point y, the density wave will be volume terminated tube, we go back to eq 61. 
edt 4 7 The solution for the density wave becomes 
+(1)'")y n=1+ (",,e7+7,_e7")e™@+ 
2jwt 
: ) € 
xo} :) = » = yus 


(n,, el) Y—1 eH L- ely 


t+ [2] "wy . , >is - : 
Layer + (65) me * eM +( H,,c881" — 
3+ 2[2]” 


zor ) ee" 4 312)" ye (2+ y 4 





J 14 , eV 41s _¢ 


Our previous conclusion permits us to assume 
that expression 64 is manageable (i. e., of limited 
variation with a time derivative of limited varia- 
tion) so that it must converge. We may therefore 
infer the following relations: 


For large enough 7 





Here 444, 7; represent the two sets of integra- 
tion constants necessary to take care of an out- 
(66) going and reflected wave in the tube. They are 
fractional density excesses. 
We will consider the boundary conditions to be, 
for the moment, at r=0 (y=0) 





9=1+9,c + 9.0" +mem+ ... (69) 
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where 7, is the amplitude of each input harmonic 
in the density wave (it differs from 7;, which 
represents both the input harmonic and the excita- 
tion amplitudes of that harmonic), and at 


z=L (y=v¥n) 


: - ¥ 
iat E Vro COSh Wro+Vp sinh =| — 


: [2]""*Wao 
Ee *VYro cosh [2]'*Yro+ 2% sinh [2]! Prot 


[2]'"*Wn* (cosh 2pro 
Ure 9/92 are 9)1/2 
2((2]"*Wro cosh [2]'*Yro4 


This result is for a given input density wave. 
It will be shown later that the results shown in 
eq 71 are only valid when the equation of condi- 
tion is isothermal (i. e. the “ polytropic”’ coefficient 
is unity). It therefore follows that if the input 


TO 


cosh [2]"*Yro) + 2b w?Wro((2]'? sinh 2¥r)—sinh [2]'“Y7) + 
- [2] "Yrown? (cosh 2vr,.—cosh [2]7) 
2¥79 sinh [2]'*Wro) (Wro cosh Prt sinh Pr»)? 


On om On 
"or wh Ot 


(see eq 1, 19, 20, and 24). 
The application of these two boundary con 
tions leads to the result that at the end of the tu! 
5 





pressure wave is given by 


b ,jwt 2h 
g= fe"! + Ee” 


the wave at the end of the tube is 


[2}'"Yro 


e=[f Vro COSh Pro+ Yo sinh = sili [2]? Yr cosh [2]? Pro t+ 210 sinh [2]'*Yro 


¢ ¢2 
$! so 


For liquids, we can start from eq 38 


32 Op 


re) ra) 
(? r) I? dt 


Or wor 
By the use of eq 22 and 23, we obtain 
O*n* = 2A On 
or? OL? Ot’ 
where 


AL 
Ar. V bry (55) 


The equation is exactly the same as before with 
the single modification that 1/6 is substituted for 
n and m, so that our previous result (eq 73) holds. 

The change in mean density along the tube can 
be estimated from eq 59 and 70. The equation 
of motion (eq 59) requires that the second deriva- 
tive of the mean square density vanishes, or that 
the first derivative is constant. However, the 
end boundary condition (eq 70) requires that the 
first derivative of the mean square density van- 
ishes at the end of the tube, and therefore along 
the entire tube, so that the mean square density 
and therefore the mean square pressure must 
remain constant along the tube. The leading 
part of the second order change in mean pressure 
arises, therefore, from the steady state portion 
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Vio(cosh 2y r~»9—cosh [2] 70) +[2]' 2W roWs0([2]'? sinh 2¥ ny — 
sinh [2] *Wr_) + ¥i.(cosh 2¥ m—cosh [2] yr) 
2(Wro cosh Yro+ Yo sinh Pro)? 


efet 4 ees 


of the square of the amplitude of the fundamental. 
It can be simply shown that the increase in mean 
pressure at the instrument is given by 


foAp[1 + | Eox/Eo\*1/4 
4. Theory Corrected for Acceleration 

In this section, we will remove the main re- 
strictive assumption—the assumed Poiseuille ve- 
locity distribution. In order to do this, it is 
necessary to go back to the equations of hydro- 
dynamics. Since the complete theory is too ex- 
tensive to be treated in this paper, we will simply 
state the results. 

It is possible to take the Navier-Stokes equa- 
tions of hydrodynamics (the equations of mo- 
tion), combine them with the equation of conti 
nuity, and with the energy equation, which repre- 
sents a detailed energy balance among thermal! 
and kinetic energies, to arrive at the Kirchoff 
(See Ravleigh, Theory of 
These equations 


equations of sound. 
sound, volume 2, article 348.) 
are valid to first order. This procedure was fol 
lowed, making no assumption as to the form of 
equation of state for the fluid, and the following 
results were obtained for the attenuation param 


eter, and the velocity in an infinite tube: 
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is the Laplacian velocity of sound in the 
fluid ; 
are arguments of the Bessel functions for 
unit tube radius; 
J, and J; are the zero™ and first order Bessel funce- 
tions; 


g and h 


o» is the mean Prandtl number of the fluid 
e 


__ “Poo \. 
( Ko ) 
Ky is the mean thermal conductivity of the 
fluid; 
Cro is the mean specific heat at constant pres- 
sure of the fluid; 
vo is the kinematic viscosity of the fluid 


“s 


1+ 








We may regard eq 78 as an extended definition 
of the attenuation parameter yr, and as_ the 
modified velocity that replaces Poiseuilles law. 
it is therefore used without the zero subscript, 
vhich is used to denote the Poiseuille regime. 

If we now bring in the end boundary condition, 
iamely 
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The attenuation parameter in eq 74 is to be 
interpreted as before (see eq 47) as the exponent 


° ty¥r—- 
in the form « L 
Equal 64 and 65 are of doubtful value for 


Vo 

. >] 
CD 
Vow 
a>. 

These restrictions are violated at high vacuum 
or very high frequencies. 

It is instructive to evaluate eq 74 and 75 for 
small values of the Bessel function arguments. 
They become 

327 wry 
C?LP 


rl* Op 
128 Ox J 





which are precisely the results assumed in eq 1 
and 42 under the condition in eq 42 that the 
“polytropic”’ This 
because the value C?/y in eq 77 is the square of 


coefficient is unity. arises 
the Newtonian velocity of sound, which for gases 
iS Po/po. Eq 74 and 75, which take into account 
the heat conduction, thus demonstrate that, when 
the previous results are valid, the equation of 
At higher frequen- 
cies the modifying term in eq 74 may be regarded 
as the “polytropic”’ 


condition is the isothermal. 


coefficient. To bring this 


out explicitly, eq 74 and 75 may be written as 


~} 
hy A(h sy) 
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for a gas, or 
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for a liquid, we obtain 


dé __12840Vjw : 
dx xrD*mpo 





for a gas, or 
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for a liquid. 
Let 


(1) 
128u)Vol . 
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for a gas, or 


(a3) 

(128u9VoLlb oo 
a rl*po oJ, (; > 
A . JTo(h ‘ 








— 


for a liquid. 

We have thus corresponding extensions of our 
definition of x, to cover all frequencies. The 
limiting value of x, for small arguments obviously 
becomes the previous value for xz. 

Since the only modification has been to extend 
the definitions of yr and ¥, without changing the 
form of the equations to be solved (namely eq 
42 and 46), the previous result (eq 49) is strictly 
valid. 
frequencies well into the audio range. 


The results however are now correct for 


It is not possible to use the results of this section 
to extend the range of validity of the calculated 
distortion for finite pressure amplitudes. To do 
this rigorously would require going back to the 
second-order terms neglected in Kirchoff’s equa- 
tions, which is an extremely arduous procedure. 


102 


It must therefore be concluded that the distortio 
calculated in section (3) is valid whenever t! 
Poiseuille regime holds, which also means th: 
the “polytropic” coefficient in the distortion mu 


be taken as unity. The distortion may be valid! 


calculated from eq 73 when 
wl? _ 
dy, = (83 
and when the applied pressure amplitude 
sufficiently small at the applied frequency (sul- 
ficiently small enough Reynolds number) to per 
mit laminar flow. 


5. Theory Corrected for Finite Length-End Effects 


There is one additional factor that must be con- 
for completeness—the end effect. An 
estimate of its magnitude will be made for the 
Poiseuille regime. It arises from the fact that it 
takes an appreciable length of tubing to set up 
the Poiseuille velocity distribution in the trans- 
mission tube. The character of the entrance flow 
is that the axial velocity is flat at the entrance, 
gradually developing an approximately parabolic 
(laminar) boundary with a core of uniform veloc- 
ity, until the approximately parabolic distribution 
fills the tube. It is evident that boundary layer 
theory may be used, and for our purposes an ex- 
tremely crude boundary layer theory. 

We go back to the equations of hydrodynamics 
and make the following assumptions: (1) that the 
entrance flow is incompressible, (2) that the varia- 
tion of pressure in the radial direction is negligible 
in the entrance portion, (3) that quadratic terms 
in velocity are negligible in the boundary layer, 
(4) that the core of the velocity distribution is 
potential. 


sidered 


For our purposes we need only write the equa- 
tions of motion for the potential core as 
ou Ou 1 Op 
P+, P ee exe 7 . 
ot Or po OF 
where u, is the axial velocity in the potential core. 
Let 


(84) 


U,* o¢ 
» Ox 
where ¢ is the velocity potential 
then 
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here f(t) is an arbitrary function of time. 
Eq (86) represents the Bernoulli integral. It 
n be written in the form 


_ 260, 1) 


g7 
dt (84) 


p(z, t) +2 pou? =f (t) 

The boundary conditions are a prescribed pres- 

ire variation at z=0, the entrance, with a flat 
elocity profile (jp=@ where @ is the average 
velocity across the section) and a parabolic dis- 
tribution of velocity at some point z=! down- 


Up 
stream 9° 


sible flow makes @ the same at both sections. 


a): The assumption of incompres- 


these conditions lead to the result that 


0¢(0, f) 


p(0, t) PL) +5 aie +29 (l,t)— ot 


(88) 

If we now refer to the arguments given in Gold- 
stein “Modern Developments in Fluid Mechanics,”’ 
vol. 1, pp 299 to 308 for the static case, we find 
on p. 302 that 


> a 16p,) u? eee, 
P(0)=pO+5 awt+ D - (1—.0575 - R.)) (89) 
ss 5 Rk, 


where R, is the Reynolds number. 

It follows from these two equations that the 
leading term for the entrance loss in the oscillatory 
case is the usual p,i? loss. 

We may therefore adopt the exact static result 
see p. 308 of Goldstein) that 


¥ro 


£ 
s 
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where the last term represents the Poiseuille vis- 
cous resistance. 

Therefore the effect of the entrance is to cause 
a pressure drop given by 


Ap=pu’, (91) 


the exact coefficient 1.2 being unimportant for 
our purposes. We will regard eq 91, not as being 
exact, but as indicating the order of magnitude 
of the entrance correction. 

If we substitute the Poiseuille velocity into eq 
(91) and evaluate the pressure gradient of the 
Poiseuille distribution from eq (47), we arrive at 
the result that the pressure excess just inside the 
tube ' is approximately given by 


¢} t eiwt — 
~ SU 


jP eo ( ro tanh ro Yoo) peti , (92) 
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Equation 92 is the desired result. It shows 
that the approximate effect of the entrance is to 
distort each input harmonic. It can be inter 
preted as meaning that the effect of the entrance is 
the same as if it did not exist, but with the funda- 
mental harmonic generator replaced by a funda- 
mental and a second harmonic generator. The first- 
order terms are thus left unaffected, and the only 
equation requiring modification is the attenuated 
second harmonic. 

If eq 92 is used as the input pressure for a pure 
sinusoidal input in the fluid conduit, eq 73 becomes 


[2]"*V rot pay 
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where 
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N is a dimensionless parameter. 


The second term in eq 93 gives the second 
harmonic distortion. 
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Actually from the condition under which eq 93 
is valid (namely eq 83), the value of N*? must be 
small, so that it is a matter of indifference whether 
it is used in eq 93 or not, and we will therefore 
neglect it. 

6. Summary 


There only remains the task of recapitulating 
the pertinent results and presenting them for 
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computational convenience. To accomplish this 
some minor notational changes will be made. 
In eq 49, it was shown that the complex attenu- 


ation of the fundamental is given by 


vr | 
¥r cosh r+ ¥~ sinh y, 


is the complex amplitude of the frac- 


(95) 


where tor 

; twWrol¥in (cosh 2Wro 
where §, is the complex amplitude of the frac- 
tional pressure excess of the second harmonic at 
the end of the tube; Wry and Wy,» are the values of 


the attenuation parameter computed on the basis 
of the Poiseuille velocity distribution. 


» 
- 
$0 /0 


In eq 76, it was shown that eq 95 is valid if 
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( 2 


In eq 83, it was shown that eq 96 is valid if 


wl)? 
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in <' (98) 


In eq 48, the attenuation parameters were de- 


fined as 


r= Xr) 
vy ix) 
In eq 78, 81, and 82, with slight modifications 


for generality, the attenuation factors were deter- 


mined to be most generally 


tional pressure excess of the fundamental at t! 
end of the tube, and the subscript O refers to t! 
fundamental. 

In eq 93, it was shown that the complex amp! 
tude of the second harmonic distortion due to 
pure sinusoidal pressure input is given in the for 
of its ratio to the input amplitude of the fund; 


mental by 
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while the attenuation factors for Poiseuille flow 


are 
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rom eq 75 
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lere Cy is the velocity of sound appropriate to the 


tube: 


(’, is the velocity of sound appropriate to the 


instrument volume. 


For computational purposes, the attenuation 
factors can be made less complicated by the 
ntroduction of two new functions. 
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where F; and F, are correction functions to the 
Poiseuille attenuation factors. 
The attenuation factors then become 


x= "C2 5) ar) FF (104) 

In order to obtain consistency with our previous 
results, we introduce the following definitions: 
For a gas: 

It can be shown both from kinetic theory and 
from experimental data that the value of the 
Prandtl number for a gas is approximately unity. 
Differences from unity are unimportant for our 
Therefore g and A in eq (102) may be 
regarded as equal. We 


F, and F, as 


purposes 


may therefore define 
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2(y- (hs ) 
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(105) 


These functions therefore depend only on two 
variables, instead of three as in eq (103). 

Equation 78 shows that the velocity of sound 
appropriate to the tube for eq 104 is the Newtonian 
velocity, which for a gas is [po/po| *. Equation 81 
indicates that the velocity of sound appropriate 
to the instrument volume for eq 104 is the ‘‘poly- 
tropic” velocity or [mpo/p|*. Therefore for a gas, 
the attenuation factors may be computed from 


; ; ) 
where F), F,, and hj may be computed from 


eq 102 and 105. 


It is convenient to utilize one more variable, 
the ratio x,/xr, Which from eq 106 has the value 


(107) 


For a liquid: 

It has been stated that y can be satisfactorily 
taken as unity for a liquid. This similarly makes 
the functions F; and F, (see eq. 103) independent 
In that 


case, the definition of F; and F, for a gas (eq. 


of g (or really of the Prandtl number). 


105) holds for a liquid, if y is taken as unity. 
Continuity of definition is thus provided for both 
liquid and gas attenuation. 

Equation 78 shows that the velocity of sound 
appropriate to the tube for eq 104 is the Newtonian 


velocity of sound However, consistent with eq. 
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55 and 58, the Newtonian velocity must be based and 
on the effective compressibility of the liquid and 
From eq 82 it is seen that the velocity of 
appropriate to the instrument 


tube. Xro 32uok ow ( D 


sound volume, 
however, is based on the real compressibility of 
the fluid (the difference 


“polytropic,”” and isothermal compressibilities is 


X19 =32 pk ow (pb ¢ 


between adiabatic, 


assumed negligible). It has been assumed that ,; 
It is now possible to compute the real attenua 


tion and lagging phase angle for the first and 
If the attenuation parameters 
in eq 95 are regarded as the low frequency param 
eters based on the Poiseuille distribution (i. e. th: 


the compressibility of the instrument volume is 
included in the definition of the effective instru- 
ment volume. It follows therefore that the 
attenuation factors for a liquid can be computed 
from 


second harmonics. 


ones with zero subscripts), then the real attenua 


no. (LN. , Po—Pe 1 Dol, D. ) tion and lagging phase angle can be computed 
xr B2ukee(p) (1 PP es Fills ” from 
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The zero subscript means that these are the 
values for the Poiseuille flow regime. Graphs of 
these equations are quite useful for computing 
attenuation. Since x7) is proportional to fre- 
quency (see eq 101), while x/x7r9 is proportional 
to the ratio of instrument volume to line volume 
(see eq 107), a family of curves of attenuation or 
phase angle plotted against y7) for different values 
Of x0/xro are frequency curves for 
different volume ratios. 


response 
These curves are pre- 
sented as figures 2 and 3. 

At higher frequencies, where the functions F; 
and F, take on values appreciably different from 
unity, the expressions become extremely com- 
plicated. It is therefore of utility to examine 
their high frequency behavior. 

At high frequency, we will use the approxi- 
mation 
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If we define a new parameter z ( related toh 4 ); 
which characterizes the fluid regime, as 


TPPw 
4y,’ 


a frequency parameter wL/C, where 
wl Xro< we 
c’ Sy : 
and a volume ratio R defined as 


R= 7%, 
XTo 


it can be shown that at high frequency 
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This is the solution for the undamped acoustic 


R — tan 


resonance of a tube and instrument. 


In computing these quantities in eq 115, it is 
that 
In the solution for 


issumed that x7) is smaller than x7 
is small compared to wL/C. 


Ze, OF 
x 
the undamped case, the phase angle lag is usually 
regarded as zero up to the first resonance. How- 
ever, the given expression permits first-order com- 
putation of the phase angle lag valid for values of 
small compared to wL/C, even though the 


vs 
overshoot is given as undamped. 


means for values of | 


Practically, this 
re] ~ of the order of one 

or less. 
These quantities are presented in figures 6 and 
It can be shown that they are valid for values 


Sor 


c 
$0 


| soak 2¢,+ cos 2¢,+ 2¢; sinh 2¢, 


tan do 


W here 


sin ¢;]' a 2s Fi] : 
-sin ¢;]”2 [ss K ] 2 


[1 


1/2 
(COS Cy [1 . sin ¢5]' 2 sin Cy [1-++sin C5!' _ | 43 F, ] 


ae. 
sin ¢;]'/7—sin ¢.[1—sin ¢s]”?) 33 F,| x 
- x 


(cos ¢,{14 


hitg 
fs T 92) 


fi 
F, 


COS ¢; 


Attenuation of Pressure in Tubes 


Lb 
*) tan +R : 
SY ( 


2ab 
0 


» 


- 


9 


2wl 
0 


R? 
4 


( 


wh 


wl 


‘a 


») 


2c, sin 2¢,+ (c?+¢ 


of z> 100, whereas the low frequency curves (figs. 
2 We will 


2 and 3) are valid for z<1 (see eq 98). 
2, which 


state without proof that the parameter 


characterizes the flow regime, is closely related to 


2 and 3 


a damping coefficient. Figures will 
therefore be referred to as the large damping 
curves, and figures 6 and 7 will be referred to as 
the undamped curves. 

Unfortunately, in many instances, a knowledge 
of the highly damped behavior (figs. 2 and 3) and 
the undamped behavior is not sufficient. Curves 
have therefore been drawn for a value of 2 about 
6.25 


2 «J 


“half-way” between 1 and 100, namely 
In order to preserve a scale 
XTo 


(see figs. 4 and 5). 


proportional to frequency, the quantity is 


used as abscissa. 


These curves were computed from the formulas 
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Equation 116 is valid at all frequencies, and is 
presented without further explanation for com- 
pleteness and the use of those with great compu- 


distortion at low frequencies (neglecting the en 
effect) are presented in figures 8 and 9. The fo 
mula used in their computation was 


tational fortitude. 
The amplitude ratio and leading phase angle 
(angle of lead on the time scale of the fundamental 
where both the fundamental and double frequency 
waves are cosine terms) of the double frequency where 
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